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The  space,  energy,  and  time  or  frequency  dependent 
Boltzmann  equation  in  the  diffusion  approximation  has  been 
used  within  the  framework  of  neutron  thermalization  theory 
to  predict  thermal  neutron  wave  propagation  in  a  moderat- 
ing medium.   This  theory  was  then  compared  with  the 
experiment  conducted  on  a  graphite  pile. 

In  the  theoretical  section  the  energy  dependence 
was  expanded  in  the  space-energy  eigenfunctions  of  the 
Boltzmann  operator,  and  the  transverse-space  dependence 
was  expanded  in  the  orthogonal  solutions  to  the  Helmholtz 
equation  with  homogeneous  boundary  conditions.   It  was 
then  shown  that  each  of  these  space  and  energy  modes  thus 
formed  propagate  with  characteristic  complex  inverse 


Xlll 


relaxation  lengths.   The  perturbation  method  developed  by 
Perez  was  used  to  expand  these  complex  eigenvalues  in  a 
Maclaurin's  series  of  the  frequency  variable.   In  cases 
where  this  series  did  not  rapidly  converge  a  direct 
solution  of  the  complex  secular  determinant  was  obtained. 
Two  scattering  kernels  were  selected  for  implementation 
into  the  theory,  the  MASSMK  kernel  and  the  SUMMIT  kernel 
based  primarily  upon  derivations  by  Parks.   The  numeri- 
cal calculations  of  the  theory  were  performed  on  an 
IBM-709  digital  computer  and  founded  upon  the  previous 
work  of  Ohanian, 

In  the  experimental  section  an  accurate  technique 
for  conducting  neutron  wave  experiments  in  multiplying  or 
diffusing  media  was  presented  along  with  the  results  of 
the  graphite  experiment.   The  excitation  of  neutron  waves, 
both  by  modulated  and  pulsed  sources,  and  the  electronic 
systems  associated  with  their  detection  were  discussed. 
Special  attention  was  given  to  the  analysis  of  data  for 
which  several  computer  codes  were  developed. 

In  the  comparison  of  theory  and  experiment  it  was 
shown  that  the  neutron  wave  technique  allows  for  the 
separation  of  diffusion  and  thermalization  effects,  thereby 
yielding  a  more  exacting  comparison  of  theoretical  and 
experimental  results.   Polynomial  fits  of  functions  formed 
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from  the  experimental  results  yielded  coefficients  to  be 
compared  with  the  theoretically  derived  terms  of  the 
perturbation  expansion.   The  results  of  the  thermal  pulse 
propagation  experiment  were  also  analyzed  by  a  method 
suggested  by  Moore,  and  confirmation  of  his  predictions 
was  established. 

The  diffusion  and  thermalization  parameters 
determined  from  the  data  analysis  were  then  compared  with 
results  of  several  die-away  experiments  and  with  the 
theory.   The  complementary  relationships  between  the 
neutron  wave  and  die-away  experiments  were  also  presented. 

The  diffusion  and  thermalization  parameters 
determined  in  this  work  for  graphite  of  density  1.60  gm/cm 
were  vZ   =  ( 91 . ll . ) sec"! ,  vD  =  ( 2. 16^ . 01 )xlO^  cm^sec'^. 


C  =  (39.12. )xlO^  cm^sec-l,  p  =  (12.12. )xlO 


"^  cm^sec"-'-. 
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CHAPTER  I 
INTRODUCTION 

Placement  of  a  modulated  or  pulsed  source  of  thermal 
neutrons  at  one  boundary  of  a  nuclear  assembly  creates  a 
disturbance  which  can  be  analyzed  in  terms  of  wave  components 
propagating  from  this  source  with  frequency  denendent  phase 
shift  and  attenuation  factors.   The  importance  of  the 
neutron  wave  experiment  is  based  upon  the  fact  that  these 
two  independently  measurable  quantities  are  strongly  dependent 
upon  the  diffusion  and  thermalization  properties  of  the 
assembly.   Thus  accurate  determinations  of  physical  con- 
stants and  severe  tests  of  theoretical  models  are  possible 
through  an  experimental  determination  of  the  frequency 
dependent  dispersive  properties  (dispersion  law)  of  the 
system. 

Several  experimental  methods  for  investigation  of 
moderating  and  multiplicative  assemblies  are  now  in  exten- 
sive use.   They  include,  among  others,  exponential  experiments, 
"poisoning"  experiments,  and  the  die-away  technique.   It  is 
interesting  to  note  that  the  neutron  wave  technique,  which 
is  the  subject  of  this  dissertation,  is  related  to  each  of 
the  three  methods  mentioned  above.   Basically,  the  neutron 


wave  experiment  is  an  extension  of  the  exponential  experiment 
in  that  the  constant  source  of  the  latter  is  replaced  by  a 
modulated  or  pulsed  source.   One  then  measures  the  complex 
inverse  relaxation  length,  p,  of  the  neutron  waves  as  they 

travel  through  the  assembly.   An  extrapolation  to  zero  fre- 

2 
quency  of  the  real  part  of  p  yields  the  same  value,  Bjj^  ,  as 

obtained  by  the  exponential  experiment.   Using  a  modulated 
source,  the  material  buckling  is  obtained  through  an  extrap- 
olation of  a  curve,  which  is  nearly  constant  near  the  zero 
frequency  intercept,  rather  than  by  measuring  a  single 
point  on  the  curve.   Also  the  imaginary  part  of  p  is  still 
available  for  determining  another  parameter  or  comparing 
with  the  theory.   Still  another  possibility  with  the  neutron 
wave  experiment  is  the  appearance  of  resonance  phenomena 
arising  from  the  splitting  of  each  spatial  mode  into  various 
energy  modes  which  propagate  with  different  complex  inverse 
relaxation  lengths. 

Increasing  the  frequency  of  modulation  in  a  neutron 
wave  experiment  can  be  considered  equivalent  to  "poisoning" 
the  medium  with  a  1/v  poison.   Therefore,  determining  the 
bucklings  and  relaxation  lengths  at  several  poison  concen- 
trations in  the  "poisoning"  experiment  and  finding  the  complex 
relaxation  lent^th  as  a  function  of  frequency  in  the  wave 
experiment  are  complementary  measurements.   Moreover,  the 
frequency-dependent,  energy  soectrum  of  the  neutro:.  wave 
experiment  bears  a  close  resemblance  to  the  energy  spectrum 


associated  with  the  "poisoning"  experiment  as  a  function  of 
the  poison  concentration.   In  the  wave  experiment,  however, 
solids  as  well  as  liquids  can  be  homogeneously  "poisoned." 
The  neutron  wave  technique  and  the  die-away  experi- 
ment also  bear  a  close  resemblance.   The  modulus  of  the 
complex  inverse  relaxation  length  can  be  thought  of  as  a 
negative  buckling.   Its  relationship  to  u/v,  the  "poisoning" 
term,  is  quite  similar  to  the  relationship  between  the  geo- 
metric buckling  and  the  decay  constant  as  measured  by  the 
die-away  technique.   In  the  die-away  experiment  one  determines 
X,  the  lowest  time-energy  eigenvalue  of  the  Boltzmann  operator, 
as  a  function  of  the  geometric  buckling.   This  method,  though 
quite  effective  for  the  determination  of  diffusion  parameters, 
has  severe  limitations  when  applied  to  the  determination  of 
thermalization  parameters  unless  the  time-dependent,  energy 
spectrum  is  measured.   At  high  bucklings,  thermalization 
effects  cause  the  curve  of  X  versus  buckling  to  deviate 
slightly  from  the  linear  relationship  predicted  by  diffusion 
theory.   This  second-order  effect  occurs  precisely  at  the 
range  of  buckling  where  the  least  confidence  can  be  given 
to  the  experimental  results;  for  it  is  in  this  range  of 
buckling  that  the  extrapolation  distance  is  most  important, 
transport  corrections  become  significant,  and  the  energy 
dependence  takes  the  longest  time  to  stabilize.   Also  an 
asymptotic  decay  is  never  attained  in  some  crystalline  modera- 
tors because  of  the  so-called  "trapped  neutrons"  which  decay 


in  time  much  slower  than  those  associated  with  the 
fundamental  eigenvalue.   Thermalization  effects  are  of 
first-order  significance  in  the  neutron  wave  experiment. 
Diffusion  theory  predicts  that  the  real  part  of  p^  should 
be  independent  of  frequency,  whereas  it  is  in  fact  nearly 
linearly  related  to  the  square  of  the  frequency.   The 
measurement  of  thermalization  parameters  via  the  neutron 
wave  technique  can  be  accomplished  in  large  systems  where 
angular  dependence  and  extrapolation  distances  are  of 
comparatively  minor  importance. 

The  purpose  of  this  dissertation  is  to  present  both 
theoretical  and  experimental  investigations  of  the  propa- 
gation of  neutron  waves  in  a  moderating  medium.   The 
Boltzmann  equation  in  the  diffusion  approximation  is  used 
as  the  theoretical  model,  and  its  energy  dependence  is 
expanded  in  the  space-energy  eigenfunctions  of  the  Boltzmann 
operator,   A  perturbation  method  is  used  to  expand  the 
square  of  the  complex  inverse  relaxation  length  in  a 
Maclaurin's  series  of  the  joj  variable  and  in  the  x  variable 
which  is  related  to  the  frequency  of  excitation  by  the 
expression  w  =  xe'^.   In  the  experimental  section  an  accurate 
technique  for  conducting  neutron  wave  experiments  in  multi- 
plying or  diffusing  media  is  presented  along  with  the 
experiment  in  graphite.   Special  attention  is  then  given 
to  the  analysis  of  the  data  thereby  obtained  and  to  its 
comparison  with  the  theory. 


CHAPTER  II 

REVIEW  OF  THE  LITERATURE 

The  propagation  of  neutron  waves  caused  by 
oscillating  a  thermal-neutron  absorber  in  a  chain-reacting 
pile  was  investigated  by  Weinberg  and  Schweinler  (1^). 
Weinberg  and  Wigner  (2)  derived  equations  for  the  atten- 
uation distance,  wave  length,  and  propagation  velocity  as 
a  function  of  frequency  for  waves  excited  by  a  localized 
source  in  an  infinite  diffusing  or  multiplying  medium. 
It  was  shown  that  the  diffusion  coefficient  and  the  prompt 
material  buckling  could  be  determined  by  measuring  the 
propagation  characteristics  of  these  neutron  waves. 
Experimental  measurements  of  diffusion  parameters  in 
graphite  and  heavy  water  were  performed  by  Raievski  and 
Horowitz  (3)  by  using  a  mechanical  source  and  a  frequency 
range  which  varied  from  50  cps  to  500  cps .   Their  results 
were  in  good  agreement  with  theoretical  predictions,  but 
no  attempt  was  made  to  extend  the  theory  to  include  therma- 
lization  effects. 


*   Underlined  numbers  in  parentheses  refer  to  the  List 
of  References. 


Preliminary  experiments  in  graphite  and  in  water- 
moderated  subcritical  assemblies  were  done  by  Uhrig  (4,5), 
He  used  one-group  diffusion  theory  to  derive  equations 
relating  the  prompt  material  buckling  and  diffusion  co- 
efficient to  the  neutron  wave's  phase  shift  and  attenuation 
and  compared  his  results  to  experiments  conducted  on  the 
graphite,  subcritical  reactor  at  Iowa  State  University  (6). 
A  mechanical  source  was  used  in  these  experiments;  and  the 
basic  unit  in  the  detection  system  was  a  dual,  coincidence- 
anticoincidence  switching  circuit  which  directed  the  detector 
pulses  to  four  slave  scalers  so  that  each  scaler  counted 
during  a  different  phase  interval  of  the  modulated  source 
(7^).   It  was  determined  that  the  waves  do  follow  the  ten- 
dency predicted  by  theory;  however,  the  wave  velocity  was 
in  all  cases  higher  than  expected,  with  the  deviation  becoming 
greater  as  the  oscillation  frequency  increased. 

The  theoretical  work  of  Uhrig  was  extended  by  Booth 
(.8)    to  a  two-group  diffusion  theory  model.   It  was  shown 
that  when  two  energy  groups  are  considered  two  thermal  waves 
propagate  in  the  fundamental  spatial  mode,  that  the  relative 
magnitude  of  the  Fermi  age  compared  with  the  square  of  the 
diffusion  length  is  of  predominant  importance  in  determining 
the  propagating  characteristics  of  the  two  waves,  that  K„ 
is  a  second-order  effect  with  respect  to  these  same  charac- 
teristics ,  and  that  the  material  buckling  can  be  measured 


from  experiments  conducted  in  the  central  region  of  the 
assembly.  Several  methods  for  the  detection  of  neutron 
waves  were  also  presented  by  Booth. 

Perez  and  Uhrig  (9_)  investigated  the  propagation 
of  neutron  waves  in  a  moderating  media  with  a  one-group 
diffusion  model  and  a  thermalization  model.   In  the 
thermalization  model  the  energy  dependence  was  expanded 
in  the  Associated  Laguerre  Polynomials  of  order  one.   The 
complex  inverse  relaxation  length  was  expressed  in  a  power 
series  of  the  ratio  u/vD,  and  it  was  shown  that  the  co- 
efficients of  this  expansion  are  functions  of  the  thermal- 
ization and  diffusion  properties  of  the  moderator.   The 
authors  demonstrated  that  the  neutron  wave  technique  affords 
a  complementary  method  to  the  "poisoning"  technique  and  that 
the  difficulties  encountered  in  the  die-away  experiment 
because  of  the  trap  effect  should  not  be  present  in  a 
neutron  wave  experiment. 

The  diffusion  and  cooling  coefficients  of  graphite 
were  determined  experimentally  by  Perez,  Booth,  and  Hartley 
(10)  by  the  neutron  wave  technique.   The  value  of  Dq  agreed 
with  previous  measurements,  but  the  value  of  Cg  was  a  factor 
of  two  larger  than  that  obtained  by  Starr  and  Price  (11) 
using  the  die-away  technique.   The  propagation  of  neutron 
waves  in  a  subcritical  assembly  using  multigroup  diffusion 
theory  was  investigated  by  Perez,  Booth,  Denning,  and 


Hartley  (1^).   It  was  shown  that  for  "n"  groups  each  spatial 
mode  splits  into  "n"  damped  waves  of  the  same  frequency  but 
with  different  phase  shift  and  attenuation  factors ,  and  that 
the  interaction  of  these  waves  presents  the  possibility  for 
measurable  resonance  effects  which  are  quite  sensitive  to 
the  lattice  parameters  of  the  system.   An  attempt  was  made 
to  demonstrate  thermalization  effects  in  moderators  by  the 
reflection  of  neutron  waves  by  Denning,  Booth,  and  Perez 
(13^).   In  the  experiment  a  sinusoidally  modulated  source 
of  neutrons  was  placed  at  one  face  of  a  26,7  cm.  long 
graphite  block  with  transverse  buckling  of  H.O  x  10"^  cm~^ 
which  preceded  a  62  cm.  long  heavy-water  region  of  the  same 
transverse  area.   Theoretical  investigations  indicated  that 
if  any  reflection  occurred  at  the  graphite,  hoavy-water  inter- 
face it  would  have  to  be  caused  by  energy  modes  which  pro- 
pagate only  when  the  energy  dependence  of  the  neutrons  differs 
from  a  Maxwellian  distribution.   This  was  qualitatively  con- 
firmed by  the  experiment. 

The  propagation  of  thermal  pulses  was  demonstrated  by 
Hartley  (14)  and  Miley  and  Doshi  (15).   Hartley  produced 
thermal  pulses  of  about  one  millisecond  width  by  use  of  a 
neutron  generator  and  thermalizing  tank,  recorded  the  pulses 
at  various  distances  from  the  source,  and  compared  the 
results  in  the  time  domain  to  a  one-group  diffusion  theory 
model.   An  extension  of  the  theory  to  that  of  Perez  and 


Uhrig  met  with  severe  difficulties  when  the  inversion  to 
the  time  domain  of  the  complicated  Laplace-transformed  flux 
was  encountered.   The  disagreement  between  the  experiment 
and  the  one-group  diffusion  theory  showed  the  sensitivity 
of  the  experiment  to  thermalization  effects.   Miley  and 
Doshi  used  the  TRIGA  Reactor  as  a  source  and  thermal  pulses 
were  produced  of  width  ranging  between  32  and  66  milli- 
seconds.  The  measured  asymptotic  velocity  of  the  pulse 
yielded  a  value  of  D^ ;  but  because  of  the  large  pulse  width 
(i.e.,  lack  of  high  frequency  content),  the  more  interesting 
thermalization  effects  observed  by  Hartley  were  not  detected, 
Diffusion  theory  was  compared  with  the  experiment,  and 
corrections  were  made  for  the  source  distribution. 

A  method  of  comparing  theory  and  experiment  in  the 
neutron  wave  experiment  with  pulsed,  thermal  sources  was 
suggested  by  Moore  (16)  so  that  the  inversion  of  the  theory 
to  the  time  domain  (which  nroved  so  difficult  for  Hartley) 
would  not  be  necessary.   In  essence,  Moore  suggested  that 
comparisons  be  made  in  the  frequency  domain  by  numerically 
determining  Fourier  moments  and  widths  of  the  data  measured 
in  the  time  domain.   Fourier  widths  were  defined  and  shown 
to  be  proportional  to  the  derivatives  with  respect  to 
frequency  of  the  dispersion  law  of  the  medium.   This  theory 
was  shown  to  be  especially  applicable  when  the  fundamental 
solution  of  the  dispersion  law  was  expressed  in  the  per- 
turbation expansion  introduced  by  Perez  and  Uhrig  (9,  17), 
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Kunaish  (1^)  used  time,  energy,  and  space  dependent 
neutron  thermalization  theory  in  the  consistent  P-^   approx- 
imation to  develop  a  theory  for  comparison  with  neutron 
wave  and  die-away  experiments.   The  non-Maxwellian  component 
of  the  energy  dependence  was  expanded  in  a  Taylor  series, 
a  heavy-gas  scattering  kernel  was  used,  and  the  results  pro- 
vided an  accurate  description  of  the  neutron  wave  experi- 
ment for  source  frequencies  up  to  500  cps. 


CHAPTER  III 


THE  BOLTZMANN  EQUATION  IN  THE  DIFFUSION 
APPROXIMATION  WITHIN  THE  FRAMEWORK 
OF  NEUTRON  THERMALIZATION  THEORY 


Presentation  of  the  Theoretical  Model 
and  the  Method  Used  for  Its  Solution 


Introduction 

Placement  of  a  sinusoidally  modulated  or  pulsed 
source  of  thermal  neutrons  at  one  boundary  of  a  moderating 
medium  generates  thermal  neutron  waves  which  propagate  from 
the  source  with  frequency  and  space  dependent  amplitude 
and  phase.   The  positional  variation  of  this  frequency 
dependent  thermal  neutron  flux  is  what  one  measures  experi- 
mentally, and  it  will  receive  the  principal  emphasis  in 
the  theory  which  follows.   Attention  will  be  given  to  the 
space,  energy,  and  time  or  frequency  dependent  thermal 
neutron  flux  in  a  moderating  medium  within  the  framework 
of  neutron  thermalization  theory.   The  main  objective  will 
be  to  calculate  the  complex  inverse  relaxation  length  as  a 
function  of  frequency  for  the  fundamental  space  and  energy 
mode , 

The  energy  distribution  of  thermal  neutrons  in  a 
finite  moderating  medium  deviates  from  a  Maxwellian  distri- 
bution for  several  reasons.   The  preferential  leakage  of 

11 


12 

higher  energy  neutrons  from  the  assembly  produces  a 
"cooling"  effect  upon  the  spectrum,  and  the  preferential 
absorption  of  slower  neutrons  causes  a  "hardening"  effect. 
The  prediction  of  this  spectrum  is  further  complicated  by 
the  thermal  velocity  of  the  scattering  nuclei  and  the 
chemical  binding  and  crystalline  effects  in  the  moderator. 
The  probability  of  a  scattering  interaction  between 
neutrons  and  the  moderator  nuclei  is  expressed  by  the 
differential  scattering  cross  section  Eg(E*,E)  for  energy 
transfer  from  E*  to  E.   It  is  the  essential  physical 
parameter  in  the  thermal  region,  and  the  prediction  of 
transient  and  asymptotic  phenomena  can  be  only  as  accurate 
as  this  scattering  kernel.   For  this  reason,  care  has  been 
taken  to  select  as  accurate  a  scattering  kernel  as  possible. 
The  more  accurate  scattering  kernels  are  much  too  complicated 
to  be  considered  analytically;  and  when  using  them,  one  must 
employ  numerical  methods. 

The  Boltzmann  Equation  in  the  Diffusion  Approximation 

To  investigate  the  effect  of  the  energy  spectrum  upon 
thermal  neutron  wave  propagation,  one  can  use  the  Boltzmann 
equation  in  the  diffusion  approximation.   For  a  finite, 
isotropic,  homogeneous,  and  non-multiplying  medium  it  is 
expressed  as 
1  3<ti(r,E,t) 


at 


D(E)v2^(r^E^t)  -  [EgCE)  +  E.(E)]4.(r,E,t)  + 


(1) 


where 
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ao 

+   /   EgCE' ,E)*(r,E' ,t)dE' 

0 

t  =  the  time  variable 

r  =  the  position  vector 

E  =  the  energy  variable 

V  =  the  speed  variable 

V^  =  the  Laplacian  operator 

(|r(?,E,t)drdE   =   the  total  flux  of  neutrons 

at  time  t  in  the  volume 
element  d'?  about  point  r 
which  have  energies  between 
E  and  E  +  dE  (neutrons/cm^/ 
sec/ev. ) 

-♦• 
E^(E)  (fr  (r ,E ,t )dE   =   the  total  number  of  absorptions 

per  unit  volume  about  r  and  per 
unit  time  at  time  t  due  to 
neutrons  with  energy  between 
E  and  E  +  dE 

Eg(E)  <|i(r  ,E,t  )dE   =   the  total  number  of  scattering 

collisions  per  unit  volume  about 
r  and  per  unit  time  at  time  t  due 
to  neutrons  with  energy  between  E 
and  E  +  dE 

Eg(E' ,E)<^(r,E' ,t)dE'dE   =   the  total  number  of  scattering 

collisions  per  unit  volume  about 
r  and  per  unit  time  at  time  t 
due  to  neutrons  of  initial  energy 
between  E'  and  E'  +  dE'  and  of 
final  energy  between  E  and  E  +  dE 

D(E)   =   the  energy  dependent  diffusion 
coefficient  (cm.) 

Within  the  diffusion  approximation,  D(E)  is  given  by  the 
expression 


D(E) 


3[E^(E)  -  Egj(E)]  (2) 
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where 


and 


^^(E)  =  Es(E)  +  E^(E)  (3) 


r  .(E)   =   the  second  Legendre  moment  of  the 
angular  dependent  scattering  cross 
section 


Eg^(E)  =  /  n.n'j:g(E,f^«n')d(n'n' ) 


-1 


From  the  definitions  of  E^CE)  and  Zg(E' ,E)  it  is  evident 
that 

Eg(E)   =   /   Eg(E,E')dE'  (4) 

0 

In  the  experiment  the  source  is  placed  near  the 
exposed,  z=0  face  of  a  long  rectangular  stack  of  cadmium- 
covered  moderator.   Thus  homogeneous  boundary  conditions  can 
be  applied  at  extrapolated  boundaries  of  the  x  and  y  variables, 
the  source  condition  is  applied  at  z=0,  and  the  flux  can  be 
specified  to  vanish  at  an  infinite  value  for  the  z  variable 
(i.e.,  end  effects  are  considered  negligible  since  measure- 
ments will  not  be  made  at  the  far  end  of  the  pile). 
Specifically,  the  boundary  conditions  that  should  be  applied 
to  equation  (1)  for  the  solution  of  this  problem  are 
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*(;a,y,z,E,t)   =   0  (5) 

^(x,  ;b,z,E,t)   =   0  (6) 

♦(x,y,«,E,t)   =   0  (7) 


and 


32 
-D(E)  - —  ♦(x,y,z,E,t) 
3z2 


=  1/2  S(x,y,E,t)       (8) 
z  =  0 


where  S(x,y,E,t)  is  the  space,  energy,  and  time  dependent 
source  of  thermal  neutrons. 

Under  these  conditions ,  it  will  be  advantageous  to 
expand  the  flux  in  an  orthogonal  set  of  functions  in  the 
X  and  y  variables.   In  view  of  the  linear  character  of  the 
problem  and  since  our  main  interest  is  focused  upon  the 
frequency  dependence  of  the  flux,  it  can  also  be  separated 
into  a  steady  part  (due  to  the  constant  component  of  the 
source)  and  an  oscillatory  part.   Considering  only  the 
oscillatoi^  component  of  the  solution  to  equation  (1),  we 
can  expand  it  in  the  form 

♦  (r.E.t)  =  Re  [  J   *i,m^'^»y^*l,m^^»^^  exp(jait)] 

t  ,m 

(9) 

where  ^^^^jCx.y)  satisfies  the  Helmholtz  equation  with 
homogeneous  boundary  conditions,  which  is  expressed  as 


(v2    ♦  b;-„   )  *    (x,y)   =   0 

x,y     l£,m    e,m  (10) 
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with 


B      r   the  transverse  buckling  of  the  i,m 
spatial  mode  and 

^      =   the  Laplacian  operator  for  the  x  and  y 
variables . 

The  substitution  of  equation  (9)  into  equation  (1) 
with  the  use  of  equation  (10),  multiplication  of  the  result 
by  <<>   )<(x,y),  and  integration  over  the  x  and  y  variables 
results  in 

C-D(E)b2     _  J.  (E)  _  J.  (E)  +  D(2)  -L-  .  lil  ]  ^   v(z,E)  + 


■»■   /   5:„(E'  ,E)*^  ,  (z,E')dE'   =   0  (11) 

o    ^        P»^ 

In  the  experiment,  one  would  like  to  measure  only  the 
fundamental  spatial  mode  so  it  is  of  major  concern  in  the 
theory.   Since  the  equations  for  the  various  spatial  modes 
differ  only  by  the  magnitude  of  the  transverse  buckling, 
one  can,  without  serious  loss  in  generality,  remove  the 
subscripts.   The  symbols  B^^  and  (>(z,E)  are  now  understood 
to  represent  the  transverse  buckling  and  the  flux  associated 
with  the  fundamental  spatial  mode,  respectively.   It  is  also 
convenient  at  this  point  to  express  equation  (11)  in  a  more 
compact,  operator  formalism  by  use  of  the  substitution* 


*   The  author  wishes  to  express  his  thanks  to  Dr.  Frank  J. 
Munno  for  helpful  discussions  in  connection  with  this 
substitution. 
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X(E,Z)  VvM(E) 

*(z,E)   =   


Vd(7)  (12) 

where  MCE),  the  Maxwellian  distribution,  is  expressed  as 

MCE)   =   C(E/kT)     expC-E/kT)  (13) 

By  substituting  equation  (12)  into  equation  (11)  and 
dividing  the  result  by  VvMCE)DCE),  one  obtains 

[.B| +  .  ]xCz,E)   + 

DCE)     DCE)    3z2    vDCE) 

I    CE'  ,E)xCz,E')  Vv'MCE') 

OD      5 

+   /     dE'  =  0 

0   VvMCE)  VDCE')  VDCE)  (14) 


or 


3^      ju) 


^-'^  '77^-  ^jFnr^^^^'"^  -^  px(z.E)  =  0        (15) 


where 


E^CE' ,E) Vv'MCE') 

PxCz,E)  =  /  dE'{  ,      ,   —  . ■ 

0      VvMCE)  VDCE')  VDCE) 


l^iE)         r„CE) 


_  6CE'-E)[ ♦  ]}  x(z,E')  (16) 

D(E)     DCE) 
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Solution  of  the  Eigenvalue  Problem  Associated  with  the 
Space  and  hJnergy  Dependent  Boltzmann  Operator 

It  is  instructive  to  consider  equation  (15)  for  the 

case  when  the  transverse  leakage  and  frequency  of  excitation 

are  both  zero  (i.e.,  a  semi-infinite  plane  with  a  constant 

source).   The  incentive  for  this  deviation  comes  from  the 

work  of  Ohanian  (19)  and  Ohanian  and  Baitch  (20)  in  their 

theoretical  investigation  of  the  die-away  experiment.   We 

shall  assume  that  the  solution  to  the  equation  thus  formed 

may  be  written  as 

x(z,E)   =    I      a^>|-j^(E)  exp[-p^^z]  (17) 

n  =  l 

with  the  understanding  that  the  above  expansion  is  convergent 
only  if  the  4i„(E)  form  a  complete  set.   Substitution  of 
the  above  expansion  into  equation  (15)  with  both  B?  and 
0)  equal  to  zero  results  in 

m  m 

-      I      ^n^no^^^^^^P^-^no^^  =  I      P^'*'n^E)exp(-p^^z) 

n=l  n=l  (18) 

If  the  initial  assumption  of  completeness  is  correct,   the 
i(pj^(E)  will  be  linearly  independent.   Therefore,  equation  (18) 
must  hold  separately  for  each  value  of  "n"  and  reduces  to 

-^no  *n^E)   =   P*^(E)  (19) 

Thus  p^^  and  ip  (E)  are  the  eigenvalues  and  corresponding 
eigenfunctions  of  the  P  operator.   An  inspection  of  this 
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operator  with  reference  to  the  principle  of  detailed  balance, 


v'EgCE' ,E)M(E')dE'dE  =  vI^CE  ,E» )M(E)dEdE' 


(20) 


reveals  that  it  is  both  real  and  symmetric.   Thus  its  eigen- 
functions  must  be  real,  satisfy  the  orthogonality  condition, 

m 

I      ii>^iZ)^^(,E)dE   =  0         for  n?<m,  (21) 

o 

and  be  linearly  independent  (2_1).   This  latter  characteristic 

assures  that  the  expansion  expressed  by  equation  (17)  and 

the  subsequent  manipulations  were  valid. 

The  facility  of  the  substitution  expressed  by 

equation  (12)  and  of  the  investigation  of  the  semi-infinite, 

zero- frequency  problem  now  becomes  evident.   If  one  expands 

the  energy  dependence  of  equation  (15)  in  terms  of  the 

eigenf unctions  obtained  above,  equation  (19)  can  be  used 

to  exchange  for  the  complicated  P<<'n^^^  expression  its  simple 

equivalent  -p^  ij,  (E).   Thus  the  solution  of  equation  (15)  by 
no  n  ^  J 

this  method  is  not  directly  complicated  by  the  consideration 
of  realistic  and  therefore  involved  energy  exchange  kernels 
since  more  difficult  kernels  simply  result  in  different 
eigenfunctions  for  use  in  the  expansion.   Or,  expressed 
differently,  the  complicated  characteristics  of  the  energy 
exchange  kernel  are  carried  to  the  solution  of  equation  (15) 
by  the  eigenfunctions  and  eigenvalues  obtained  from  the  P 
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operator.    Moreover,  the  homogeneous  equations  and  secular 
determinant  resulting  from  the  above  mentioned  expansion 
will  be  much  simpler  than  that  obtained  by  using  other 
orthogonal  functions  which  do  not  bear  as  close  a  relation- 
ship to  the  actual  solution. 

The  eigenfunctions  and  their  corresponding  eigen- 
values also  have  physical  significance.   They  are  the 
thermalization  theory  solution  to  the  exponential  experi- 
ment in  a  semi-infinite  plane.   The  ij/,  (E)  function  is  the 
fundamental  energy  mode  which  is  exponentially  attenuated 
with  distance  from  the  source,  and  its  corresponding  eigen- 
value is  the  square  of  the  asymptotic  inverse  relaxation 
length.   The  one-group  diffusion  theory  equivalent  to  this 
eigenvalue  is  1/L^.   Thus  ijj^CE)  and  p^   represent,  respec- 
tively, the  eigenfunctions  and  eigenvalues  of  the  space  and 
energy  dependent  Boltzmann  operator  and  are  what  one  en- 
counters when  performing  the  exponential  experiment  in  a 
very  large  moderating  medium.   More  important  physical 
significance  can  be  assigned  to  these  eigenvalues  in  the 
context  of  this  work  when  one  recalls  that  the  exponential 
experiment  can  be  considered  a  neutron  wave  experiment 


This  argument  is  based  upon  the  assumption  that  the 
difficulty  involved  in  the  solution  to  the  eigenvalue 
problem  is  not  dependent  upon  the  complexity  of  the 
scattering  kernel.   This  is,  in  general,  the  case  since, 
as  explained  in  the  second  part  of  this  chapter,  a 
discrete  energy  representation  of  the  eigenvalue  problem 
is  solved  on  a  computer.   Thus  more  involved  scattering 
kernels  simply  result  in  different  numbers  being  read 
as  input  to  the  code. 
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performed  at  zero  frequency.   From  this  viewpoint  it  becomes 
apparent  that  the  p2^  are  the  zero  intercepts  of  the  squared 
complex  inverse  relaxation  lengths  plotted  as  a  function  of 
frequency.   A  knowledge  of  these  zero  intercepts  is  useful 
because  their  closeness  to  the  fundamental  eigenvalue  is 
quite  indicative  of  the  amount  of  higher-energy-mode  con- 
tamination predicted  by  the  particular  scattering  kernel 
under  investigation. 

Further  insight  into  the  nature  of  these  eigen- 
functions  and  eigenvalues  can  be  obtained  by  a  more  detailed 
investigation  of  the  eigenvalue  problem  represented  by 
equation  (19).   with  the  P  operator  explicitly  expressed 
and  some  terms  rearranged,  equation  (19)  becomes 

^^^^  0   VvM(E)Vd(F)  A/D(E') 

(22) 

For  the  special  case  in  which  E^(E)  =  0,  the  lowest  eigen- 
value  is  zero  and  its  corresponding  eigenfunction  is 

VvM(E)D(E);  which  is  easily  verified  by  direct  substitution 
and  by  use  of  the  principle  of  detailed  balance.   This  is 
consistent  with  the  fact  that  in  an  infinite  medium  with 
no  absorption  the  flux  will  attain  an  equilibrium  Maxwellian 
distribution  in  energy  and  a  constant  value  in  space.   From 
the  arrangement  of  equation  (22)  it  is  evident  that  it  will 
be  satisfied  by  two  types  of  solutions: 
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1.  p^  i    I^(E)/D(E)  for  any  value  of  E.   This 

no    '- 

results  in  discrete  eigenvalues  and  the  more  conunon 
continuous  eigenfunctions  as  solutions, 

2.  p^^  =  j:-t;(E)/D(E) .   In  this  case  the  expression 

enclosed  in  brackets  is  zero.   The  right-hand  side  of 

equation  (22)  is  not  necessarily  zero,  however,  so  that  if 

the  eigenfunctions  are  to  satisfy  the  equation  in  this 

circumstance  they  must  exhibit  in  part  a  singular  behavior. 

Moreover,  condition  (2)  can  be  satisfied  for  any  value  of 

the  continuous  variable  E  under  the  condition  p^   > 

no 

[E^(E) /D(E)  ]j^jjj.   Thus  above  this  limit  point  a  continuous 
eigenvalue  spectrum  is  obtained. 


Expansion  of  the  Frequency  Dependent  Equation  in 
^igenfunctioi    r-  . ,  '  „ 
Boltzmann  Opt 


Eigenfunctions  of  the  Space  and  Energy  Dependent 
ioitzmann  Operator 


The  energy  dependence  of  the  flux  can  now  be  developed 
in  terms  of  the  eigenfunctions  mentioned  above.   An  exponential 
dependence  upon  the  z  variable  will  be  assumed  so  that  the 
solution  to  equation  (15)  can  be  expressed  in  the  form 

CO 

x(z,E)  =   J   An<J'n(E)exp(-pz)  *  (23) 

nal 


*  A  more  rigorous  expression  for  this  expansion  would  be 

N 
X(Z,E)  =  I      {A^d^^d>(E)  *    I      A^^\p2)/c)^j,^^2)dp2}^^pf.p,3 

n=l 


Ce^(E)/D(E)]^jj^ 
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where  the  i>    (E)  satisfy  the  eigenvalue  equation,  and  they 


have  been  normalized  so  that 


(2^) 


Substitution  of  equation  (23)  into  equation  (15)  results 
in 

I      [-b2_  +  p2 1!1-  ]  Aj^(exp[-pz]),|;^(E)   + 

+  I      A^(expC-.pz])Pt);^(E)  =  0  (25) 

n=l 

which,  when  multiplied  by  'i^jjj(E)  and  integrated  over  all 
energy,  yields 


n  =  l 


0 


m  =  1,2,3,...  (26) 


where  the  superscripts  "d"  and  "c"  represent,  respectively, 
the  discrete  and  continuum  sets.   However,  in  the  actual 
solution  a  discrete  energy  formalism  is  used  and  calcula- 
tions are  performed  only  at  a  finite  number  of  discrete 
energy  values.   Eigenvalues  in  the  continuum,  which 
actually  are  a  function  of  energy,  and  their  corresponding 
eigenfunctions  are  obtained  only  at  particular  energy 
points.   Thus  equation  (23)  correctly  expresses  the  actual 
expansion  used,  which  is  the  discrete  approximation  to  the 
more  rigorous  and  less  convenient  representation. 
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where  use  has  been  made  of  equations  (19)  and  (24)  and 

00 

^mn   =   /   ^m^E)  [l/vD(E)]  4'^(E)dE 

0  (27) 

Equation  (26)  expresses  in  a  compact  form  a  set  of  linear, 
homogeneous,  algebraic  equations  for  the  modal  amplitudes 
Aj^  and  the  eigenvalues  p^.   One  is  able  to  solve  the  set 
only  for  the  relative  magnitudes  of  the  A^^,  and  the  number 
of  eigenvalues  which  allow  a  solution  will  be  equal  to  the 
number  of  equations  considered.   Also  there  will  be  a 
different  set  of  ratios  for  each  eigenvalue.   The  compatibility 
condition  for  the  above  set  of  homogeneous  equations  yields 
the  secular  determinant  or  dispersion  law 

I  ^B^  -  Pno^^n  -  >L    I   =   0  (28) 


no  mn    -^   mn 


with 


b2   =   p2  _  B?  (29) 


which  can  be  used  to  determine  the  complex  eigenvalues  that 
will  allow  solutions  for  equation  (26).   These  eigenvalues 
are  the  squared  complex  inverse  relaxation  lengths  which 
are  to  be  compared  with  the  experiment.   It  is  satisfying 
to  note  that  for  u)=0  the  secular  determinant,  equation  (28), 
reduces  to  a  diagonal  determinant  yielding  the  p^   for  the 
B^  eigenvalues. 
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Determination  of  the  Fundamental  Complex  Inverse 
Relaxation  Length  from  the  Complex  Homogeneous  Equations 

The  perturbation  method  developed  by  Perez  and 
Uhrig  (9_)  and  systematized  by  Kunaish  (18_)  will  be  used 
to  develop  each  of  the  eigenvalues  and  modal  amplitudes 
in  a  power  series  expansion  of  the  frequency  variable. 
This  power  series  representation  is  used  because  it  allows 
one  to  compare  very  easily  various  theoretical  models 
with  one  another  and  with  the  experimental  results.   Once 
the  coefficients  are  calculated,  the  eigenvalue  can  be 
determined  for  any  frequency  within  the  radius  of  con- 
vergence of  the  series.   To  initiate  this  expansion,  the 
first  equation  and  the  first  term  of  each  of  the  other 
equations  in  the  homogeneous  set  represented  by  equation 
(26)  are  explicitly  expressed.   Considering  the  first  "N" 

equations  of  this  set,  one  obtains 

N 
(b2  -  pJq)  -  L^^j^  -  I      4„ju,S^  =  0 

n=2 

for  m  =  1  (30) 

N 
-I^lico  -  I      {(b2  .  o2^)6^  -  L^j.  }  S^  =  0 
ns2 

for  m  =  2,3,4,...       (31) 

where  each  equation  has  been  divided  by  A-,    so  that  the  set 
could  be  solved  for  the  ratios 
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n 


\ 


/A. 


n  =  2  ,3,4, 


(32) 


and  for  the  eigenvalue  B^,   In  this  method,  B^  and  each 

of  its  corresponding  ratios,  S  ,  are  expanded  in  a  Maclaurin's 

series.   Thus  we  let 


B' 


I 

I 

K=0 


B^^h].)^ 


K! 


(33) 


and 


n 


I 

K  =  0 


K! 


(34) 


where 


B 


(K) 


.K 


3^ju) 


jojsO 


(35) 


and 


,(K) 


n 


a^ju 


n 


jojrO 


(36) 


The  K"*^^  coefficient  of  B^  and  those  of  each  of  the  S^^  are 
obtained  from  equations  (30)  and  (31),  respectively,  after 
differentiating  them  "K"  times  with  respect  to  jw  and 
setting  jti)  equal  to  zero. 
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It  should  be  emphasized  at  this  point  that  we  are 
presently  expanding  only  the  fundamental  eigenvalue  along 
with  its  corresponding  ratios.   However,  it  will  be  shown 
that  the  other  N-1  eigenvalues  and  their  ratios  are  obtained 
in  exactly  the  same  way  after  rearranging  the  initial 
equations  and  redefining  some  terms. 

By  letting  jiorO  in  equation  (30)  and  equation  (31), 
it  can  easily  be  determined  that 


B^°^   =   p2 


10 


(37) 


and 


^n     =   0  n  =  2,3,M,...  (38) 


Using  the  identity 


~  [A(x)B(x)]   =    I      (J)  aK(x)b(I-K)(x) 


3X 


where 


KrO  (39) 


(J) 


If 


K!(I-K)!  (40) 


one  can  differentiate  equation  (30)  and  equation  (31)  "I" 
times  to  obtain  at  ju)  =  0 


K=0  n=2 

(41) 
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n=2  K=0 

m  =  2,3,4,...        (42) 


where  u(I-2)  is  the  unit-step  function.  By  evaluating  the 
summation  over  "K"  in  equation  (41)  and  rearranging  terras, 
one  obtains 


B 


N 

10  10     11  II  ^   ''l^   In  n 

"*2  (43) 


which  can  be  used  to  determine  each  B^    if  the  S^'^'"'" 

n 

are  known.   Eqi  ation  (42)  can  be  simplified  by  separating 
out  from  the  rest  of  the  series  in  the  "K"  summation  the 


terms  for  which  K=0  and  K=l.   By  performing  this  manipulation, 

hicl 
(I) 


extracting  the  S    coefficient  from  the  term  for  which 


K=0,  and  rearranging  terms,  an  expression  for  the  S 
ratios  is  obtained 

N 


L   6    -  u(I-2)  I      (i)(B^^^    _  L   )S^^ 
mi  li  -   ■*■       mn    mn  n 

,(I)  n=^ 

'm    =   


(p2   -  p2  ) 

mo    10 
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"(1-3)   I    I   (J)B«'s<^-'^' 

ns2  K=2 


(p2   -  p2  ) 

mo    10 


m    -    2,3,4, 


(41+) 


The  summation  over  "K"  in  the  above  expression  is  extended 

only  to  "I-l"  since  S^   is  zero,  and  the  step  functions 

are  introduced  because  some  terms  were  extracted  from  the 

summations . 

The  procedure  used  for  evaluating  the  B^^^  and  S^^^ 

m 

is  initiated  by  first  determining  B^^^  and  S^"^^  from 

m    ^ 

equations  (43)  and  (44),  respectively.   Knowing  these  terms, 

(  2 ) 

B    can  be  calculated  from  equation  (43)  and  placed  in 

equation  (44)  in  order  to  find  the  S^2)^   ^^.^  ^^^^^^^   ^f 

alternately  using  the  two  expressions  to  first  calculate 

B  ^   then  S^^^  when  all  the  B<^^  and  S^J^  are  known  for 

m 

J<I  can  be  repeated  to  determine  as  many  terms  as  one  needs. 
Then  using  the  B^    in  equation  (33),  the  frequency  dependent 
complex  inverse  relaxation  length  for  the  fundamental  mode 
can  be  evaluated. 

It  may  happen,  however,  that  the  series  does  not 
converge  very  rapidly.   In  fact,  this  is  the  case  for  one 
of  the  scattering  kernels  selected.   For  this  situation,  a 
more  convenient  series  can  be  obtained  by  a  change  of 
variable  of  the  form 
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jw   =   Jxe^  (45) 


One  can  then  derive  a  power  series  expansion  for  B^  in 
terms  of  x*   Recognizing  that 


d^ 


dx 


Cixe'^) 


X=0  (46) 


and using  the  procedure  as  explained  with  the  w  expansion, 


(^^^  and  the  S^^^ 
m 

the  power  series  expansions  in  the  x  variable.   As  in  the 


one  can  obtain  expressions  for  the  B^-'-''  and  the  S^^^   for 

m 


u  expansion 


B^O^   =   p2 


10  (47) 


and  all  of  the  S    are  zero.   For  1=1,  one  finds 


■^  11  (48) 


^m^   =   ^^ml/^^mo  "  'lo^  ^   =    2,3.4....       (49) 


and ,  in  general ,  for  I  >  2 


n=2  K=l  (50) 
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and  I-l  N   I-l 

K=l  n=2  K=l 


g(I) 

ra 


'"^o   -»Jo> 


m  =  2,3,4,..-       (51) 

In  the  determination  of  p^  as  a  function  of  frequency,  both 
of  these  power  series  expansions  will  be  used. 

An  alternative  to  the  series  representation  is  to 
solve  the  complex  secular  determinant  directly  for  several 
values  of  w.   For  reasons  previously  discussed,  this  method 
is  not  as  appealing  as  the  ones  derived  above;  and  it  will 
be  used  only  as  a  last  resort,  that  is,  for  cases  in  which 
both  of  the  series  converge  slowly  or  not  at  all. 

Determination  of  the  Higher  Order  Complex 

Inverse  Relaxation  Lengths 

It  was  stated  that  the  equations  obtained  for  the 
B    and  the  S^   were  derived  to  evaluate  only  the  funda- 
mental p  .   To  calculate  the  other  eigenvalues ,  one  can 
rearrange  the  equations  expressed  by  equation  (26)  so  that 
the  same  expressions  can  be  used  to  calculate  coefficients 
associated  with  all  of  the  eigenvalues.   To  explain  this 
procedure,  consider  the  case  where  m=3  in  equation  (26). 
Allowing  n  to  range  from  one  to  three,  one  obtains 
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C(p2  .  B|  -  p2^)  -  L^iJ'»']A^  -  L^2^u,A2  -  ^g^A.   =  0 


3 

(52) 


-  L2;l^''^1  *    ^^^^    '  ^i  "  ^20^^  ■  L22^'^-'^2  "  ^23^"^3  ""  ° 


(53) 


-  hl^^\    -    ^2^*^^   -  ^^''^  -  ^1-  ^30^  -  hz^'^^^    -    ° 


(54) 


Let  U3  specify  that  the  eigenvalue  associated  with  the 
second  equation  is  desired  (i.e.,  the  eigenvalue  whose 
zero  intercept  is  0^^+  Bj),      If  we  exchange  equations  (53) 
and  (52)  and  rearrange  terms  in  each  of  the  equations,  the 
result  is 


[(p2  .  b2  .  p^2^)  .  L^^j.lA^  .  L^^j.A^  .  L23J.a,A3  =  0 

(55) 

-  42^'^^2  '    ^''    -  ^i  -  'lo^    '    4l^«^^  -  43^^^3  =  ° 

(56) 


-  ^2^  ^2  -  Sl^-^  -  ^^<>"    -  ^i  -  ^30^  -  S3^-^A3  =  ' 


(57) 


An  inspection  of  these  last  three  equations  and  the  ones 
from  which  they  were  derived  indicates  that  both  sets  are 
of  exactly  the  same  form  since  they  differ  only  in  sub- 
scripts.  Thus  the  same  expressions  used  for  the  solution 
of  the  fundamental  eigenvalue  are  also  valid  for  the  second 
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if  care  is  taken  to  correctly  rearrange  the  L„„  and  p^ 

-^        °       nrn      no 

matrix  elements  and  to  assign  subscripts  to  them  correspond- 
ing to  their  new  positions  before  applying  the  equations. 
By  analogy,  the  arrangement  of  terms  for  calculating  the 
n   eigenvalue  and  its  ratios  is  obtained  by  exchanging 
the  first  and  the  n"*^^  equations  and  then  exchanging  the 
first  and  n   terms  in  all  of  the  equations. 

Numerical  Calculations  Associated 
with  the  Theory 

Introduction 

The  derivations  thus  far  presented  are  founded  upon 
the  theoretical  formalizations  of  Perez  an  i  'J!-irif  (9)  and 
Ohanian  and  Daitch  (2_0),   The  work  of  the  former  authors 
indicated  the  method  by  which,  beginning  with  the  Boltzmann 
equation  in  the  diffusion  approximation,  one  could  formulate 
a  power  series  expansion  for  the  complex  inverse  relaxation 
length.   They  used  Associated  Laguerre  Polynomials  in  ex- 
panding the  energy  dependence,  whereas  here  the  eigenfunctions 
of  the  space  and  energy  dependent  Boltzmann  operator  form 
the  orthogonal  functions  for  this  expansion.   This  more 
elegant  method  has  practical  value  since  the  secular  deter- 
minant used  to  determine  the  eigenvalues,  B^   has  a  much 

n  * 

simpler  form  and  one  need  not  calculate  the  complicated 
matrix  elements  of  the  scattering  operator.   As  a  by- 
product of  this  method,  one  obtains  all  of  the  zero  inter- 
cepts of  the  B^  (i.e    o^  ) 
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For  the  determination  of  the  B^   and  their 

n 

corresponding  eigenfunctions ,  use  was  made  of  the  numerical 
techniques  developed  by  Ohanian  for  his  calculation  of  the 
eigenfunctions  of  the  time  and  energy  dependent  Boltzmann 
operator.   In  fact,  after  minor  modification,  several  of 
the  computer  codes  developed  by  Ohanian  were  used  in  this 
work. 

Evaluation  of  the  Eigenvalues  and  Eigenfunctions  of 
the  Space  and  Energy  Dependent  Boltzmann  Operator 
in  the  Diffusion  Approximation 

The  Discrete  Energy  Formulation  of  the  Eigenvalue 

Problem  —   In  the  discrete  energy  representation,  equation 

(15)  becomes  a  set  of  "J"  equations  of  the  form 

-  P^r.'^J^A^    =   I   K(E.,E.).j;  (E.)AE.  -  [  ^  J^n^E-) 


i  =  1,2,3,. ..,J 


(58) 


where 


K(Ej,Ei) 


Eg(Ej,E^)Vv.M(E. 


) 


Vv^M(E^)  Vd(E.)Vd(E^) 


(59) 


J 
t 


^^(Ei)   =    I         Eg(E.,E,)AE, 

t=l  (60) 
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AE.  depends  upon  the  numerical  integration  scheme  selected, 
"J"  is  equal  to  the  number  of  energy  points,  and  the  prime 
on  the  summation  indicates  that  the  term  involving 
j:  (E^.E^)  has  been  eliminated  from  the  summation  in  the 
interest  of  greater  numerical  accuracy  (1_9_)  ,   An  inspection 
of  the  above  equations  with  all  terms  included  reveals  that 

for  i  =  i  the  term  E  (E.  ,E.)  in  the  j  summation  in  equation 

•'  s   1   1 

(58)  should  be  canceled  by  an  identical  term  in  the 

J:<,(E--) 
expression   ?.  .-^   .   It  was  determined  by  Ohanian  that 
D(Ei) 

greater  accuracy  could  be  achieved  by  elimination  of  these 
terms  altogether  before  introducing  equation  (58)  into  the 
computer.   This  is  a  very  fortunate  development  since,  as 
explained  later  in  this  section,  it  was  necessary  to  adjust 
the  diagonal  elements  of  the  scattering  kernel  calculated 
by  the  SUMMIT  code.   Since  the  diagonal  elements  are  not 
needed  in  this  numerical  solution  to  the  eigenvalue  equation, 
this  modification  should  not  have  any  direct  influence  upon 
that  solution. 

The  set  of  equations  represented  by  equation  (58) 
can  be  more  simply  expressed  as 

-  p2  (*  )   =   (P)(AE)(.{/  )  (61) 

no  n       —       n 


where 


( i|;  )  is  a  column  vector  composed  of  the 
values  of  the  i|)n  eigenfunctions  at 
each  of  the  j  energy  points, 

(aE)   is  a  diagonal  matrix  containing  the 

weights  for  the  numerical  integration. 
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and 

(P)   is  a  square  matrix  representing  the 
""   P  operator. 

The  elements  of  the  (P)  matrix  are  given  by  the  expression 

J' 

5ij     I   ^s^^i'^i^^^i    J^a^Ei) 

Pij   =   K(Ej,E.)(l  -  6ij)  -  —  [  i^i .  ] 

^^j         D(E^)         D(E^) 

(62) 

We  have  now  arrived  at  the  equation  which  is  to  be 
solved  by  the  EIGENVALUE  code.   Note  that  the  synunetry 
of  the  P_  operator  is  preserved  in  the  discrete  representa- 
tion and  that  since  (aE)  is  simply  a  diagonal  matrix  to 
provide  weighting  for  the  integration  the  numerical  eigen- 
value problem  also  contains  a  symmetric  operator. 

Evaluation  of  the  Scattering  Kernel  and  Formulation 
of  the  Symmetrized  P  Operator  —   All  of  the  theoretical 
calculations  were  performed  on  an  IBM  709  digital  computer. 
In  essence,  the  implementation  of  the  theory  derived  in  the 
preceding  section  was  accomplished  through  the  use  of  an 
echelon  of  computer  codes,  beginning  with  one  which  calcu- 
lated the  scattering  kernel  and  ending  with  the  EIGENVALUE 
code.   Each  code  made  use  of  the  results  of  one  or  more  of 
its  predecessors  and,  in  turn,  contributed  to  the  codes 
farther  on  in  the  chain.   The  format  for  the  discussion  in 
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this  section  can  therefore  be  logically  based  upon  the 
flow  diagram  shown  in  Figure  1,"   In  the  discussion  which 
follows,  pertinent  information  pertaining  to  the  computa- 
tional methods  will  be  presented.   For  a  complete  discussion 
of  the  workings  and  input-output  information  for  the  com- 
puter codes  used  the  reader  is  referred  to  the  reports 
submitted  by  the  author  to  his  committee  supervisor. 

One  of  the  more  important  considerations  in  the 
implementation  of  this  numerical  method  was  the  determina- 
tion of  what  variable  to  use  in  order  to  most  advantageously 
express  the  energy  dependence.   This  problem  was  solved  by 
Ohanian  Cl9_)  by  the  use  of  a  variable  defined  as 

n   =   (E/KT)^^^  /  (1  +  (E/KT)^/2)  (63) 

The  33  energy  points,  corresponding  to  32  equally  spaced 
increments  of  the  n  variable,  that  he  used  were  adopted  for 
this  work  and  are  listed  in  Table  1,   Also  shown  in  the 
table  are  the  D(E^)    as  obtained  by  the  use  of  SUMMIT, 
PINNACLE,  and  DIAGONAL  CORRECTOR  codes.   For  the  numerical 
integrations  encountered,  except  those  in  SUMMIT  and 


*   This  computation  was  done  under  the  direct  supervision  of 
Dr.  M.  J.  Ohanian,  with  initial  inspiration  from  Dr.  R.  B, 
Perez.   The  author  is  grateful  to  Robert  B.  Razminas  for 
the  use  of  his  DIAGONAL  CORRECTOR  CODE  and  for  his  con- 
tribution to  the  matching  of  several  of  the  codes. 


38 


0)  ^^ 
13    (t 


t4 

a 

^  e 

O    <D 

CO  v^ 
CO 

s:  y; 


c 

o 

•H 
fJ 

JD 
■H 

U 
4-) 

OT 
•H 

n 

> 

o 

S  w 
cr  c 

5^    O 


c 
o 

•H 

c 
a 

X 

w 
to 
o 
t^ 

> 

O 


c 
o 


en 
<u 

4-' 
•H    0) 

■•-'  e 
a  ft) 
c  t- 

Cu 


c 
o 
c 
o 


Cu  u 


o 

•H 

o  o 
o  c 

W    O 

M  O 


o 


•H 

c- 

•P 

< 

o 

■M 

•> 

a 

CO 

E 

4-> 

>> 

c 

CO 

•H 

< 

O 

^^ 

CL 

!W  r:^ 

t-j 

bT 

tl 

^ 

(U 

F^ 

C 

y: 

[-J 

c 


3 

< 

> 


o 

c 


— ?r 

•H     ^ 
^   -M 

gl< 

W 
N 

< 

U 
O 


e 


o 

CO 


4-' 


W     w      U       -H 

CO     ■•-'        W^ 

t-J     t-0      1-;    O 


W 
•H 


< 
O 


c 

M 


<     C< 

M     O 


to 


2: 

CO 


a 
•n 
W 

«k 

•H 

W 

CO 


W 


u 


0) 
4-' 

•r4 

x:  CD 

(1)    fO   c 

+J 

rO  4-. 

P4      O 

a; 

cu  en 
e  to  (u 

(U   TJ    c 

E-  £  W 


O 

CU 

bC 


C/3 
C/5 
< 


CO 

■p 

C    c 

•H   < 

o    * 


O  bC 
>       P^    D 

v^        (1)        •< 


U 

w 
or: 


w 


w 


w 


(U 
■M 
•H 

C.T  x: 

C    It 
•H    f^ 

o  u 

>^    O 

bO 
U  CO 

(U    CO 
C   fO 


< 
2 


>, 

(4 

o 

0) 

^ 

Eh 

0) 

,c 

+J 

c 

•H 

■r 

• 

(U 

d) 

en 

•p 

^ 

0) 

r^ 

CO 

a 

<D 

E 

■V 

o 

o 

o 

(J 

p 

{-. 

o 

CD 

c 

P 

3 

0) 

c 

c^ 

E 

q; 

O 

5 

C' 

H 

(U 

M 

^ 

r 

p 

r 

3 

U- 

cn 

o 

E 

rH 

0 

fC 

t. 

Jh 

U-. 

0) 

> 

tn 

CD 

p 

m 

c 

0) 

U-: 

C 

o 

CD 

rH 

E 

(U 

m 

u 

iH 

be- 

fO 

lt 

c 

•H 

o 

c 

bO 

(0 

3 

•tH 

0 

a 

rH 

39 
lABLE  1 


THE  ENtRGY  POINTS  USED  IN  FHt  DISCRETt  ENERGY  REPRESENTATION 
ANU  THE  CALCULATEIj  VALUES  Of     THE  OIEFUSION 
CUfcFFICIENI  AT  THESE  POINTS 


ENERGY  POINTS 

DIFFUSION  CUEH 

(tV.) 

(CM.  ) 

U. 

.  19A8989^>E- 

-0-!. 

0. 

.90540472E 

00 

0. 

.824/n09)E- 

-04 

0, 

.  16644209E 

01 

0, 

.  l966b200E- 

-03 

0. 

.28357907E 

01 

0. 

.17ill399E- 

-03 

0. 

.  3a6'i6l64E 

01 

u. 

.61667299E- 

-03 

0. 

.50007^96E 

01 

0. 

.94622^0GE- 

-03 

0. 

.61531631E 

01 

0, 

.l37b2lOJE- 

-02 

0, 

.  7J0)2Z7tit 

01 

0. 

.  I9221999t- 

-02 

0. 

.51939103E 

OU 

0. 

.26096600E- 

-02 

0, 

,  70397882E 

00 

0. 

,i<»648699E- 

-02 

0. 

.94295391b 

00 

0. 

.45211899C- 

-02 

0. 

.  12733748t 

01 

0, 

.b8l96499E- 

-02 

0. 

.960b0945b 

00 

0, 

.7^il029  vt- 

-02 

0. 

.95736114b 

00 

0. 

.93586800E- 

-02 

0, 

.91lU5106b 

00 

u. 

.11742199E- 

-01 

0. 

.11272739b 

01 

0. 

.i4662799E- 

-01 

0. 

.107ia94lt 

01 

0, 

.  1824950JE- 

-01 

0. 

,8l394893b 

00 

0, 

.22669999E- 

-01 

0. 

.79333499E 

OU 

0. 

.28I43399E- 

-01 

0. 

.84320298L 

00 

0. 

.34960400E- 

-01 

0. 

,82992154b 

00 

0. 

.43bl2300b- 

-01 

0. 

,82778580t 

00 

0. 

.t>4334600E- 

-01 

0. 

.809?7708b 

00 

0. 

.6817309:>E- 

-0  1 

0. 

.77l86675t 

00 

0. 

.860a9300C- 

-01 

0. 

,786706L0b 

00 

0. 

.  10963000E- 

-00 

0. 

,79197136E 

00 

0. 

.14111499E- 

-00 

0. 

,  77184059E 

00 

0, 

.  1841360UE- 

-00 

0. 

.78138231b 

00 

0. 

,24448000E- 

-00 

0. 

,  76lb2067t 

00 

0. 

,33I9160UE- 

-00 

0. 

,78947944b 

00 

0. 

,46393500E- 

-00 

0. 

,91681817E 

00 

0. 

.67426600E 

00 

0. 

,  96415816b 

00 

0. 

,  10345899E 

01 

0. 

,  I48b92  7/b 

02 

0. 

,1719159 >t 

01 

0. 

,lb7?G6l0G 

01 
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PINNACLE,  Simpson's  rule  was  used  since  most  of  the  functions 
to  be  integrated  were  generally  smooth.   Comparatively  few 
quantities  were  needed  as  input  parameters ,  and  the  more 
important  ones  are  listed  in  Table  2. 

Figure  1  shows  the  sequence  of  the  numerical  cal- 
culations and  the  input  and  output  of  each  code  in  the  chain. 
For  the  first  task,  that  of  calculating  the  scattering  kernel, 
two  different  models  were  used.   The  first,  and  more  realistic, 
kernel  was  formulated  by  Parks  (22,23)  and  incorporated  by 
Bell  (2^)  into  a  computer  code  called  SUMMIT.   The  P,  kernel 
was  calculated  from  the  output  of  SUMMIT  by  the  PINNACLE 
code.   The  code  multiplied  the  angular  dependent  scattering 
kernel  from  SUMMIT  by  the  cosine  of  the  scattering  angle, 
performed  a  Gaussian  integration  over  the  angle,  and  multi- 
plied the  result  by  am  appropriate  energy  dependent  factor. 
The  second  scattering  kernel  considered  was  a  free-gas  energy- 
exchange  kernel  obtained  by  the  use  of  the  MASSMK  code  ( 2_5)  . 
Note  that  the  calculations  involving  the  MASSMK  kernel 
differed  from  those  of  the  SUMMIT  code  only  by  the  scattering 
kernel  Es (Hi ,Ej ) .  All  other  inputs  to  the  OMEGA-ZERO  SYMBOL 
via  DIAGONAL  CORRECTOR  (26_)  were  the  same  as  those  obtained 
from  the  SUMMIT  kernel. 

The  scattering  kernel  calculated  by  SUMMIT  does  not 
include  the  elastic  contribution.   For  this  reason,  it  was 
necessary  to  correct  the  diagonal  elements  of  this  kernel 
so  that  the  scattering  cross  sections  obtained  from  it  were 
consistent  with  values  taken  from  the  "barn  book"  (BNL-325), 
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TABLE  2 


PARAMETERS  USED  Ifl  THE  NUMERICAL 
CALCULATIONS 


Parameter 


Mass  of  Graphite 

Oo 

U    (number  density) 

kT 

An 

°s   FREE 


Value 


12.011 

90. OU 
2.198x105 

0.0803 
0.02526 
0.027027 
4.6949 


Units 


amu. 

sec"  1 
cm/sec 

n/cm^ . 

ev. 


barns 
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This  was  accomplished  in  the  DIAGONAL  CORRECTOR  code. 
As  shown  in  Figure  1,  several  other  quantities  were  also 
obtained  from  this  code.   The  (P)  matrix  was  calculated 
using  the  OMEGA-ZERO  SYMBOL  code.   This  code  was  obtained 
by  modifying  the  SYMBOL  code  written  by  Ohanian  to  account 
for  the  differences  between  the  (P)  matrix  of  equation  (52) 
and  the  one  obtained  by  him  for  his  eigenvalue  problem. 

The  eigenvalues  and  eigenfunctions  of  the  space  and 
energy  dependent  Boltzmann  operator  were  computed  by  the 
EIGENVALUE  code  using  the  Jacobi  method  of  matrix  diagonali- 
zation  (Appendix  A  of  reference  (19)).   This  was  not  the 
same  code  as  used  by  Ohanian  but  it  was  similar  to  it. 
The  code  used  in  this  work  was  modified  to  accept  the 
diagonalization  when  the  smallest  diagonal  element  became 
more  than  10-^   greater  than  the  largest  off-diagonal  element. 

Presentation  of  the  Eigenfunctions  and  Eigenvalues 
Several  of  the  eigenfunctions  and  all  of  the  eigenvalues 
computed  by  the  method  outlined  above  are  presented, 
respectively,  in  Figure  2  and  Table  3  for  the  SUMMIT 
kernel  and  Figure  3  and  Table  3  for  the  MASSMK  kernel. 
Some  problems  were  encountered  in  drawing  the  eigenfunctions 


*  The  code  was  obtained  by  Dr.  Ohanian  and  converted 
for  use  on  the  University  of  Florida  IBM  709  computer 
by  the  author. 
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(.2) 


(33)xl/2 


Log  of  Enerj^y 
(uev. ) 


(31) 


(11) 


(21) 


2.   Several  Fif^enfunctions  of  the  Space  and  Energy  Dependent 
Boltzmann  Operator  Using  the  SUMMIT  Kernel. 


4U 


Fig.  3.   Several  Ei;Tenfunctions  of  the  Space  and  Ener^^y  Dependent 
Boltzmann  Operator  Using  the  MASSMK  Kernel. 
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TABLt  3 


CIGENVALUES  UF  FHt  SPACE  AND  ENERGY  DEPGNDGNT 
BULrZMANN  OPERATOK  USING  THE  SUMMIT  KERNEL 
AND  T^^H    MASSMK  KERNcL 


SUMMIT  MASSMK 

0.^284^2^lE-03  0. 429942796-03 

U. i087b46yE-02  0 . 30b94949E-02 

0.  170 7  865 OC-02  0 .  14  3^0  7 I2C-01 

0.44850244E-0?  0 . 6302 51 52E-0 I 

0.6595399dE-02  0. 8041^B6 /G-01 

0.  10662339E-01  0 .  107349 I6E-00 

0.  1494969)E-01  0 . 984748 lOE-Ql 

0.17229941E-01  0.  1230726bE-00 

0.2OO7674yE-0l  0. 1444 r561 E-00 

u.i9b29l6yE-01  0.17164934E-00 

0.20320657E-01  0.227286 I3t-00 

U.2l37806dE-0l  0  .  244'i2724E-00 

0.2230108i£-01  0. 268 J1973C-00 

0.2451l52iE-01  0.  1626110UE-00 

0.27752671E-01  0. 302  l4746t-00 

0.29222453E-01  0. 344t.9a67E-0G 

0.35146919E-01  0. 32484750E-00 

0.4047440bE-01  0  .  3'>9'i3601E-00 

0.366l90i)7t-0i  0.39137711t-0G 

0.43478405E-01  0. 4092651 8E-00 

0.50642820E-01  0  .  383c!d654E-00 

0.49833021E-01  0 .43537266E-00 

0.55298758fc-01  0 .44520015Q-00 

0.6501758 3E-01  0. 439/5324E-00 

0. 77065729E-01  0.4478059  IE-GO 

0.81092229E-01  0 .4498 3647E-00 

0.9l62a064E-0l  0  .476  I0663b-00 

0.  11717452E-0U  0.4824065iE-00 

0.  12339253t-00  0.50220024b  00 

0. I0396006C-00  0.71235961E  00 

0.  13176729E-00  0.10231198C  01 

0.  14078926C-00  0.12931190E  OL 

0.19005135E-00  0.48972284E  01 
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because  they  can  have  more  than  one  singularity  and  in  the 
discrete  representation  it  is  difficult  to  distinguish 
between  multiple  discontinuities  and  oscillations.   The 
possibility  for  multiple  discontinuities  is  attributed  to 
the  many  fluctuations  in  E-,- ( E) /D(  E)  .   This  function  can  be 
identical  to  a  particular  eitrenvalue  for  several  energy 
values  thereby  causing  the  corresponding  eigenf unction  to 
have  several  discontinuities.   From  the  characteristics  of 
the  eigenfunctions ,  it  appears  that  only  the  fundamental 
eigenvalue  is  discrete;  and  the  remaining  ones  represent 
particular  eigenvalues  in  the  continuum.   As  shown  in  Table 
3,    beyond  p2g  the  larger  eigenvalues  associated  with  the 
SUMMIT  kernel  are  closer  to  the  fundamental  eigenvalues  than 
are  the  ones  from  the  MASSMK  kernel. 

The  Energy  Expansion  of  the  Frequency  Dependent 
Boltzmann  Equation  in  the  l3if fusion  Approximation 

The  expansion  of  the  energy  dependence  of  equation 
(15)  in  terms  of  the  eigenfunctions  of  the  P  operator  re- 
quires the  evaluation  of  only  one  matrix  element,  Ljjy^, 
expressed  by  equation  (27).   In  the  discrete  energy  repre- 
sentation this  integration  becomes 


J 

^mn   =    I   V^i^  ^ ' 3  *ntEi)AEi 

i=l  v(Ei)D(Ei)  (64) 
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where  t^L^   is  the  weighting  factor  for  the  numerical  integration 

A  L-MATRIX  ELEMENTS  code  was  written  which  evaluated  these 

matrix  elements  using  both  Simpson's  rule  and  Weddle's  rule 

(27)  and  which  also  checked  the  orthonormality  condition. 

The  numerical  integrations  were  evaluated  using  the  n  variable 

since  the  eigenfunctions  obtained  from  the  EIGENVALUE  code 

were  expressed  in  this  variable  and  the  discrete  energy 

representation  was  based  upon  the  selection  of  equal  n 

increments  for  the  evaluation  of  the  energy  points.   It  was 

felt  that  since  each  of  the  continuum  eigenfunctions  exhibited 

at  least  one  discontinuity,  and  since  DCE^)  is  not  a  smooth 

function,  it  would  be  necessary  to  apply  very  elaborate 

integration  schemes  in  order  to  obtain  accurate  results.   If 

any  integration  other  than  Simpson's  rule  was  used,  it  was 

found  that  for  some  of  the  eigenfunctions  the  orthonormality 

check  was  not  satisfied  (i.e.,  the  integration  over  energy  of 

all  i(/_(E.)  times  iJ'„(E-)  was  not  necessarily  much  less  than 
m   1         n   1  -' 

the  integration  of  all  l^    (.L-)']^)  .      This  fact  is  not  surprising 
since  the  eigenfunctions  obtained  from  the  EIGENVALUE  code 
were  obtained  with  a  Simpson's  rule  weighting  function. 
For  this  reason,  it  was  felt  that  continuity  of  the  numerical 
methods,  that  is,  the  use  of  Simpson's  rule  throughout,  would 
result  in  the  more  consistent  set  of  L^^   elements.   Some 
justification  for  this  reasoning  is  found  in  the  work  of 
Ohanian  (1_9^)  since  a  similar  integration  was  required  in  his 
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calculation  of  the  cooling  coefficient  for  water  and  accurate 

it 

results  were  obtained, 

A  pragmatic  consideration  of  this  calculation  provides 
some  satisfaction  since  its  use  resulted  in  theoretical 
curves  which  were  in  good  agreement  with  the  experiment. 
Also  the  diagonal  elements  of  the  L   matrix  were,  in 
general,  much  larger  than  the  off-diagonal  elements  since 
in  the  evaluation  of  the  former,  the  large  discontinuity 
of  both  eigenfunctions  occurred  at  the  same  energy  point. 
Thus  after  a  careful  investigation  of  the  integrands  arising 
in  the  evaluation  of  the  L^jjj^  matrix  elements,  some  justifi- 
cation can  be  obtained  for  their  accuracy.   But,  in  the 
general  case,  the  numerical  integration  of  such  "wild" 
functions  by  Simpson's  rule  cannot  be  defended. 

Evaluation  and  Presentation  of  the  Complex 
Inverse  Relaxation  Lengths 

With  the  Lj^   matrix  elements  determined,  one  can 

apply  equations  (43)  and  (44)  to  evaluate  the  coefficients 

of  the  expansion  of  p^  in  powers  of  (ju)  and  aquations  (50) 

and  (51)  to  evaluate  the  same  expansion  in  powers  of  x» 

In  cases  where  neither  series  converges  rapidly,  the 


Also  in  c.ilcalat  ion-  oer  formed  by  Razrr'.ir..i-;  d.id  the  author 
under  the  direction  of  Ohanian,  the  cooling  coefficient 
for  graphite  was  calculated  by  Ohanian 's  method.   Tt 
agreed  quite  well  with  the  one  obtained  by  Honeck  using 
the  same  kernel. 
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eigenvalues  were  evaluated  directly  from  the  secular 
determinant  at  particular  values  of  frequency.   Computer 
codes  were  written  to  calculate  the  coefficients  for  both 
the  Jul  expansion  and  the  x  expansion  of  any  of  the  33  p^ 
along  with  their  corresponding  ratios,  S^.   If  more  than 
a  few  terms  were  used  in  the  ju>  expansion,  however,  it 
was  necessary  to  divide  the  equations  by  L]^]^  and  derive 
expressions  for  the  expansion  in  powers  of  L2^]_jci).   This 
was  because  at  high  frequencies,  u^   became  larger  than 
10"^^,  the  largest  number  acceptable  to  the  computer.   This 
problem  did  not  arise  in  the  x  expansion  since  for  a 
frequency  of  about  1000  ops  x  need  be  only  about  7. 
Problems  were  also  encountered  in  the  coefficients.   For 
the  ju)  and  L]_]_ju)  expansions,  the  series  could  not  be 
carried  past  the  fifteenth  terra  since, in  the  former  case, 
underflow  occurred  in  the  computer  aind  in  the  latter  case, 
overflow  occurred.   In  the  x  expansion  the  coefficients  did 
not  change  as  rapidly  and  no  underflow  or  overflow  problem 
occurred;  but  since  factorial  33  is  nearly  10   ,  the  series 
was  cut  off  at  the  33rd  term.   This,  however,  was  not  the 
limiting  factor  in  the  x  expansion  for  it  appeared  that 
after  the  20th  term  round  off  errors  built  up  rapidly. 
Shown  in  Figure  *+  and  Figure  5  are  the  frequency 
dependent  real  and  imaginary  components  of  the  complex 
inverse  relaxation  length  for  the  fundamental  energy  mode 
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and  the  first  four  energy  harmonics  as  generated  by  using 
the  MASSMK  kernel.   Presented  in  Figure  6  and  Figure  7 
are  these  same  eigenvalues  up  to  the  third  harmonic  ob- 
tained by  use  of  the  SUMMIT  code.   Figure  8  and  Figure  9 
show  the  fundamental  and  first  two  space  harmonics  for 
the  MASSMK  calculation.   The  transverse  bucklings  used  in 
this  last  set  of  curves  are  indicated  on  the  figures ,  and 
for  all  other  curves  a  value  of  5.95  x  10~^  cm"^  was  used 
since  this  was  the  calculated  value  for  the  graphite  pile 
investigated.   The  omega  expansion  was  used  in  the  genera- 
tion of  all  the  MASSMK  curves,  and  the  method  presented 
previously  for  the  generation  of  higher  energy  harmonics 
was  successfully  applied  with  this  scattering  kernel.   Six 
terms  were  sufficient  for  convergence  in  this  case.   For 
the  SUMMIT  kernel,  however,  the  omega  expansion  did  not 
converge  when  as  many  as  15  terms  were  used;  and  for  this 
scattering  kernel  the  x  expansion  and  the  direct  numerical 

solution  of  the  complex  secular  determinant  yielded  the 

J* 

eigenvalues."   As  shown  in  Figures  5  and  7,  the  18-term 

X  expansion  appears  to  have  converged  since  agreement  exists 
between  it  and  the  direct  solution  for  the  fundamental 


*      The  code  used  for  this  calculation  was  obtained  by 

Victor  Cain  from  the  3HARE  Program  Library.   The  code, 
called  riGEMVALUr.S  0^  COMPLEX  MATRICES,  has  as  its 
identification  number  NU  EIG4. 
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eigenvalue.   The  higher  energy  mode  eigenvalues  for  the 
SUMMIT  kernel  were  calculated  by  the  direct  method. 

Some  concern  can  be  caused  by  a  careful  inspection 
of  the  higher  energy  modes.   For  higher  frequencies,  these 
modes  appear  to  be  less  attenuated  than  the  fundamental 
mode;  and,  therefore,  modal  contamination  is  predicted. 
The  less  accurate  MASSMK  kernel  does  not  predict  as  serious 
energy  mode  contamination.   One  must  recall,  however,  that 
in  order  for  these  higher  energy  modes  to  propagate,  they 
must  be  excited  by  the  source  or  the  medium;  and  even  in 
the  finite,  absorbing  medium  to  be  investigated  they  will 
account  for  only  a  small  percentage  of  the  energy  spectrum 
compared  with  the  fundamental  Maxwellian  distribution. 
Also  Kunaish  has  shown  that  for  the  heavy  mass  scattering 
kernel  the  total  energy  spectrum,  calculated  from  as  many 
as  nine  of  his  energy  harmonics,  will  propagate  with  a 
complex  inverse  relaxation  length  which  is  within  1 
per  cent  of  that  of  the  fundamental. 

From  the  curves  showing  the  higher  spatial  mode 
eigenvalues,  it  can  be  concluded  that  for  this  experiment 
they  will  be  of  secondary  importance  in  comparison  with  the 
higher  energy  modes.   Calculations  by  Perez  (9)  and  Kunaish 
(18)  predict  that  higher  spatial  contamination  should  be 
reduced  to  a  negligible  amount  at  points  farther  than  20 
cm.  from  the  source.   The  calculations  presented  here 
confirm  these  predictions. 


CHAPTER  IV 

THE  EXPERIMENT  IN  GRAPHITE 

Introduction 
The  experimental  techniques  currently  used  at  the 
University  of  Florida  for  the  production  and  detection  of 
neutron  pulses  and  waves  were  developed  over  a  period  of 
nearly  four  years  through  the  efforts  of  many  persons. 
As  a  half-time  graduate  assistant  on  the  neutron  wave 
project,  the  author  has  been  intimately  involved  in  all 
phases  of  this  developmental  work.   The  methods  for 
analyzing  the  data  have  also  evolved  over  a  period  of 
years.   The  principal  contribution  of  the  author  has  been 
in  this  area;  for  he  has,  under  the  guidance  of  his  chair- 
man and  with  suggestions  from  co-workers ,  written  every 
computer  code  required  for  the  data  analysis  and  analyzed 
or  co-analyzed  all  the  data  taken  thus  far  in  the  neutron 
wave  project.   Also  this  work  represents  by  no  means  the 
first,  nor  the  last,  experimental  results  to  be  derived 
from  the  neutron  wave  project,  as  is  evident  from  an  in- 
spection of  the   Review  of  the  Literature. 
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The  medium  investigated  in  this  experiment  was  a 
rectangular  pile  of  AGOT  graphite  UU.7  cm.  by  71  cm.  by 
125  cm.  covered,  except  on  the  source  face,  with  30  mil. 
thick  cadmium.   The  discussion  of  the  experiment  which 
follows  is  divided  into  three  general  areas.   Consideration 
is  given  to  the  methods  used  to  excite  the  thermal  neutron 
waves.   Attention  is  then  focused  upon  the  systems  employed 
in  their  detection.   The  third  section  deals  with  the 
analysis  of  data  and  will  be  given  special  emphasis. 

The  Source 
The  source  of  thermal  neutron  waves  used  in  this 
experiment  was  a  Cockcroft-Walton  type  neutron  generator 
placed  in  a  thermalizing  tank  as  shown  in  Figure  10.   The 
neutron  generator,  manufactured  by  Texas  Nuclear  Corporation 
of  Austin,  Texas,  produced  either  a  sinusoidally  modulated 
or  a  pulsed  intensity  of  14  Mev  neutrons  which  were  then 
slowed  down  to  thermal  energy  and  directed  to  the  graphite 
pile  by  the  thermalizing  tank.   The  sinusoidally  modulated 
intensity  of  neutrons  was  obtained  by  varying  the  extraction 
voltage  of  the  ion  source  bottle,  thereby  controlling  the 
intensity  of  the  deuterium  beam  which  impinged  upon  the 
tritium  target.   The  pre-acceleration  pulsing  unit  consists 
of  deflection  plates,  Einzel  lens,  and  focus  lens.   This 
unit  deflects  the  deuterium  beam  to  the  base  of  the  Einzel 
lens  during  the  wait  time  between  pulses  ajid  allows  it  to 
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Fig.  10.   Arrangement  of  the  Experiment  and  Details 

of  the  Thermalizing  Tank  (from  Reference  14) 
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pass  on  to  the  target  for  the  duration  of  the  pulse.   For 
the  data  obtained  with  the  modulated  source,  the  frequency 
varied  between  50  and  1226  cps ,  the  ratio  of  maximum  to 
minimum  neutron  count  rate  was  roughly  9  to  1,  and  harmonic 
and  60  cycle  contamination  were  quite  small.   For  data  taken 
in  the  pulsing  mode,  the  maximum  to  minimum  thermal  neutron 
yield  was  about  1000  to  1;  and  the  total  pulse  width  was 
about  1  millisecond  at  the  z  position  in  the  graphite  which 
was  closest  to  the  source.'   More  detailed  discussions  of 
the  sinusoidal  and  pulsed  mode  are  given,  respectively,  in 
References  (8_)  and  (14). 

As  previously  mentioned,  the  therraalizing  tank  was 
used  to  therraalize  the  14  Mev  neutron  waves  produced  by  the 
neutron  generator  and  then  transport  them  to  the  external 
port  of  the  tank  where  the  graphite  medium  was  placed.   This 
was  accomplished  by  placing  light  water,  steel  bricks,  lead, 
and  graphite  in  the  tank  and  facing  the  tank  with  a  poured- 
paraffin  section.   By  systematically  varying  the  dimensions 
of  the  various  materials  and  carving  out  a  8  cm.  by  8  cm. 
square  hole  in  the  center  of  the  paraffin  section,  the 


*  It  should  be  mentioned  that  the  attainment  of  such  a 
high-quality  source  was  due  primarily  to  the  efforts 
of  Robert  H.  Hartley  and  George  W.  Fogle. 
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uncovered-to-covered  count  rate  of  a  1  atm.  He  detector 
was  optimized  at  50  to  1  at  the  source  interface  of  the 
graphite.   The  initial  thinking  associated  with  this 
thermalizing  tank  is  presented  in  Reference  (8),   The 
unwanted,  but  ever  present,  fast  background  emitted  from 
the  tank  was  eliminated  experimentally  as  discussed  below 
and  in  the  section  dealing  with  the  data  analysis. 


The  Accumulation  of  Data 
It  is  fitting  at  this  point  to  indicate  precisely 
what  data  are  to  be  presented  and  how  it  was  recorded.   The 
main  objective  of  the  neutron  wave  experiment  is  to  measure 
the  frequency  dependent  complex  inverse  relaxation  length, 
p ,  of  the  fundamental  space  and  energy  mode.   The  most 
obvious  method  to  accomplish  this  objective  is  to  modulate 
sinusoidally  the  neutron  current  which  enters  one  boundary, 
z=0,  of  the  medium  under  investigation  and  to  measure  the 
amplitude  and  phase  of  this  neutron  wave  as  it  propagates 
from  this  source.   The  slope  of  the  semilog  plot  of  amplitude 
versus  z  is  the  real  part  of  p,  or  the  attenuation  factor, 
a;  while  a  plot  of  the  phase  of  the  wave  versus  z  is  linear 
with  a  slope  equal  to  the  imaginary  part  of  p ,  or  the  phase 
shift  factor,  5.   The  source  is  modulated  for  as  many  fre- 
quencies as  desired  by  the  experimenter.   The  process  out- 
lined above  was  used  to  obtain  the  wave  data  presented  in 
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this  work.   Still  another  way  of  obtaining  the  same  information 
is  to  record  the  dispersion  of  a  pulse  of  thermal  neutrons  as 
it  propagates  from  the  source  boundary.   The  amplitude  and 
phase  of  any  frequency  component  of  this  pulse  at  a  particu- 
lar z  position  is  obtained  by  numerically  performing  a  Fourier 
transform  of  the  measured,  time  response  at  that  position. 
These  methods  are  not  new.   The  analogy  in  electrical  theory 
to  the  discussion  above  is  that  the  frequency  characteristics 
of  a  "black  box"  can  be  determined  by  either  directly  measur- 
ing its  response  at  various  frequencies  or  by  performing  a 
Fourier  transform  of  its  unit  impulse  response.   One  method 
is  not  necessarily  superior  to  the  other  for  all  applications, 
but  in  this  measurement  it  was  found  convenient  to  obtain 
most  of  the  data  by  exciting  the  neutron  waves  v;ith  thermal 
pulses , 

The  pulse  propagation  data  were  collected  by  Hartley, 
and    first  appeared  in  his  master's  thesis  (14),   The 
purpose  of  his  investigation  was  accomplished  by  considering 
data  in  the  time  domain  only.   In  this  work  the  analysis  of 
Hartley's  data  is  extended  beyond  that  which  he  completed 
since  Fourier  transforms  are  performed  under  the  guidance 
of  Reference  (15)  and  the  complex  inverse  relaxation  length 
is  obtained.   Actually,  in  the  final  comparison  with  the 
theory,  the  wave  data  and  pulse  propagation  data  are  combined 
as  discussed  in  the  next  section. 


6U 

From  the  above  discussion  it  is  evident  that  it  will 
be  necessary  to  discuss  two  data  accumulation  systems,  that 
associated  with  a  sine  wave  input  and  the  one  used  by 
Hartley  in  his  pulse  propagation  experiment.   Before  con- 
sidering these  systems  individually,  however,  it  is  con- 
venient to  present  those  characteristics  which  were  common 
to  both  experimental  set-ups. 

In  each  of  the  experiments  use  was  made  of  a  movable 
detector  which  could  be  located  at  several  z  positions,  and 
a  stationary,  standard  detector  placed  in  the  thermalizing 
tank.   To  facilitate  accurate  placement  of  the  movable  detector 
and  cause  as  little  perturbation  as  possible,  length-wise 
holes  were  bored  in  several  of  the  graphite  stringers  to 
accomodate  the  detector  and  its  preamplifier.   For  measure- 
ments at  a  particular  z  position,  a  detector-filled  stringer 
was  inserted  into  the  graphite  pile  so  that  its  center  was 
aligned  with  the  transverse  center  of  the  pile  and  its  length 
was  in  the  transverse  plane.   By  using  several  of  these 
stringers  with  holes  bored  at  different  distances  from  their 
front  face,  the  detector  could  be  moved  from  the  source  in 
about  1  cm.  increments  from  1  cm.  to  72  cm.   Two  measurements 
were  made  at  each  position  of  the  movable  detector;  one 
where  both  epicadmium  and  thermal  neutrons  entered  the 
moderating  medium,  the  uncovered  run,  and  a  second  where 
thermal  neutron  penetration  was  blocked  with  a  cadmium  shutter, 
the  covered  run.   The  standard  detector  was  used  to  correct 
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for  source  variations  and  to  facilitate  subtraction  of  the 
epicadmium  background  from  the  movable  detector  data  of 
the  uncovered  run,  thereby  obtaining  the  thermal  wave.   The 
same  source  (i.e.,  neutron  generator  and  thermalizing  tank) 
was  used  in  both  experiments ,  and  changing  from  the  wave 
experiment  to  the  pulse  experiment  was  accomplished  by 
changing  accordingly  the  mode  of  operation  of  the  neutron 

generator.   The  same  detector  systems  were  also  common  to 

3 
both  measurements.   They  consisted  of  He   detectors  and  Nuclear- 

Chicap:o  Corporation  oreampli  f  iers  and  scalers.   The  system 
dead  time  was  about  3  microseconds,  which  was  found  to  be 
adequate  for  our  purposes.   The  arrangement  of  electronics 
beyond  the  scalers  was  unique  for  each  of  the  two  experi- 
ments and  they  will  now  be  discussed  separately. 

Electronic  Equipment  Associated  with  the  'Vave  Experiment 

A  block  diagram  of  the  electronics  used  in  the  wave 
experiment  is  shown  in  Figure  11,   The  main  component  of  this 
system  is  the  multi-channel  analyzer  manufactured  by  Technical 
Measurement  Corporation.   A  modified  212  plug-in-unit  is 
used  in  the  multi-channel  analyzer  so  that  it  stores  the 
input  pulses  in  a  channel  corresponding  to  their  time  of 
arrival  after  an  initial  trigger.   The  analysis  time  of 
each  of  the  channels  is  stepwise  variable  from  10  to  2,560 
microseconds  so  that  the  accumulation  may  cover  a  wide 
range  in  frequency.   The  remaining  circuitry  is  arranged  so 
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Fig.  11.   A  Block  Diagram  of  the  Klectronic  Equipment  Used 
in  the  '/ave  Experiment, 
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that  for  two  periods  of  the  input  sine  wave  the  output 
of  the  standard  detector  and  that  of  the  movable  detector 
are  recorded,  respectively,  in  the  first  half-memory  and 
second  half-memory  of  the  multi-channel  analyzer.   This  is 
accomplished  by  alternately  allowing  the  pulses  from  each 
detector  to  reach  the  212  plug-in-unit ,  controlling  their 
storage  in  the  multi-channel  analyzer,  and  adjusting  its 
channel  width  so  that  the  total  sweep  time  is  slightly  less 
than  four  periods  of  the  input  frequency. 

A  time-sequence  diagram  for  the  system  is  shown  in 
Figure  12.   The  HP202A  function  generator  provides  a  pulse 
at  the  positive  peak  of  its  sine  wave  output  to  the  neutron 
generator.   This  pulse  train  is  fed  to  a  Tektronix  type 
161  pulse  generator  which,  when  triggered  by  a  sync,  pulse, 
produces  another  pulse  whose  width  is  almost  as 
long  as  two  periods  of  the  input  sine  wave.   This  output 
pulse  from  the  Tektronix  type  161  pulse  generator,  which 
is  triggered  by  every  other  sync,  pulse,  serves  two  purposes 
It  is  used  as  a  trigger  for  a  Tektronix  type  16  3  pulse  gen- 
erator whose  output,  connected  to  the  trigger  input  of  the 
212  plug-in-unit,  initiates  the  sweep  of  the  multi-channel 
analyzer.   Thus  the  analyzer  sweep  always  begins  at  the 
same  point  in  the  cycle  of  the  input  wave.   The  pulse  train 
is  also  connected  to  the  external  trigger  input  of  the 
Tektronix  type  5  3/5UC  plug-in-unit,  and  it  controls  the 
switching  of  this  gate  circuit.   Other  inputs  to  the  Tek- 
tronix type  53/5UC  plug-in-unit  are  the  pulses  from  the 
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Fig.  12.  Time  Sequence  of  Electrical  Signals  Associated 
with  the  Electronic  Equipment  Used  in  the  Wave 
Experiment. 
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standard  detector  and  those  from  the  movable  detector. 
The  Tektronix  53/5HC  plug-in-unit  contains  two  identical 
amplifier  channels;  and  it  is  usually  used  with  a  dual 
trace  oscilloscope  so  that  two  signals,  one  connected  to 
input  A  and  the  other  to  input  B,  appear  on  alternate 
sweeps  of  the  oscilloscope.   In  this  application  the  pulses 
from  the  two  detectors  are  allowed  to  pass  alternately  to 
the  input  of  the  212  plug-in-unit.   An  alternate  output 
from  the  53/54C  plug-in-unit  is  connected  to  the  EXT  input 
of  the  multi-channel  analyzer,  and  it  is  used  to  control 
the  half-memory  in  which  the  detector  pulses  are  being 
stored.   Thus  the  53/54C  plug-in-unit  controls  which 
detector  pulses  are  being  recorded  by  the  multi-channel 
analyzer  and  also  controls  the  half-memory  in  which  the 
pulses  are  being  stored.   This  eliminates  the  possibility 
of  the  storage  locations  for  the  two  detectors  being  inter- 
changed during  the  accumulation  of  data. 

Referring  again  to  Figure  12,  the  functioning  of  the 
detection  system  during  one  analysis  sweep  of  the  multi- 
channel analyzer  is  as  follows : 

1.   A  trigger  from  the  Tektronix  type  15  3  pulse 
generator  starts  the  multi-channel  analyzer  on  its  analysis 
sweep.   The  pulses  from  the  standard  detector  will  be  stored 
in  its  first  half-memory,  and  then  it  will  stop  and  wait 
for  another  trigger. 
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2.  The  next  trigger  will  start  the  multi-channel 
analyzer  again,  but  now  the  Tektronix  type  53/54C  plug-in- 
unit  is  passing  only  the  movable  detector  signal  and  holding 
the  multi-channel  analyzer  in  the  second  half  of  its  memory. 

3.  The  next  trigger  will  again  start  the  cycle, 
which  is  repeated  until  enough  statistics  have  been  obtained. 

Electronic   Equipment  Associated  with  the  Thermal 
f-^eutron  Pulse  Propagation  Experiment 

A  block  diagram  of  the  circuitry  used  by  Hartley 

is  shown  in  Figure  13.   For  his  experiment,  only  the  integral 

counts  of  the  standard  detector  were  needed.   The  output  of 

the  movable  detector  was  recorded  in  the  multi-channel  analyzer 

after  shaping  in  a  Tektronix  type  153  pulse  generator.   The 

sweep  of  the  multi-channel  analyzer  was  initiated  by  a  signal 

synchronized  with  the  start  of  the  neutron  burst  from  the 

neutron  generator. 

The  Data  Analysis 

Introduction 

The  focal  point  in  the  analysis  of  neutron  wave 
data  is  the  counts  per  channel  recorded  for  the  movable 
detector  by  the  multi-channel  analyzer.   To  indicate  the 
typical  appearance  of  these  data,  the  200  cps  thermal  wave 
is  shown  in  Figure  14  for  three  of  the  positions  used  in 
the  wave  experiment;  and  the  thermal  pulse  at  several 
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positions  is  shown  in  Figure  15.   In  the  wave  experiment 
these  data  are  composed  of  a  fundamental  sine  wave  con- 
taminated by  harmonic  distortion,  and  a  Fourier  series 
expansion  yields  the  amplitude  and  phase  of  the  fundamental 
thermal  wave.   The  thermal  pulse,  on  the  other  hand,  was 
considered  to  be  composed  of  neutron  waves  in  a  continuous 
spectrum  from  z;ero  frequency  up  to  a  high  frequency  cutoff. 
Thus  the  amplitude  and  phase  of  a  particular  frequency  was 
obtained  by  numerically  Fourier  transforming  the  pulse  at 
the  frequency  of  interest.   After  the  attainment  of  the 
amplitude  and  phase  points  versus  z  position  for  several 
frequencies,  the  data  analysis  and  comparison  with  theory 
can  proceed  in  an  identical  manner  for  both  the  wave  and 
pulse  propagation  experiments.   In  fact,  in  this  work  the 
amplitude  and  phase  curves  obtained  by  the  pulse  propagation 
experiment  and  those  from  the  wave  experiment  are  normalized 
together  before  determining  their  slope.   In  this  way  the 
complex  inverse  relaxation  length,  p,  arises  from  measure- 
ments in  both  the  time  domain  and  the  frequency  domain.   The 
fact  that  these  two  sets  of  data  can  be  normalized  together 
effectively  before  proceeding  with  the  evaluation  of  p 
is  significant  since  it  indicates  that  the  basic  assumptions 
associated  with  their  analysis  were  correct  and  also  points 
out  the  fundamental  equivalence  in  these  two  methods  of 
excitation  of  neutron  waves. 
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Fig.  15.   Thermal  Pulse  for  the  Pulse  Propagation  Experiment 
for  Several  z  Positions  (from  Reference  14 ) . 


75 

In  the  discussion  which  follows  in  this  section, 
the  data  evaluation  will  be  carried  through  the  calculation 
of  the  real  and  imaginary  components  of  the  complex  inverse 
relaxation  length.  The  data  acquired  from  wave  experiments 
will  be  considered  separately  from  those  of  the  pulse  pro- 
pagation experiment  in  obtaining  the  amplitude  and  phase  as 
a  function  of  position  and  frequency.  Then  these  data  will 
be  normalized  together  for  the  evaluation  of  p. 

The  data  presented  here  has  been  taken  over  a  period 
of  two  years  in  two  separate  wave  experiments  and  in  Hartley's 
pulse  propagation  experiment.   The  curves  of  amplitude  and 
phase  versus  z  position  have  been  determined  for  55  fre- 
quencies ranging  froit  zero  to  1450  cps.   For  each  of  these 
curves  a  minimum  of  14  z  positons  have  been  taken;  and  in 
many,  where  all  three  experiments  have  been  normalized  to- 
gether, up  to  40  z  positions  are  considered.   Moreover,  to 
obtain  each  of  the  14  to  40  points  on  the  130  amplitude  and 
phase  curves  either  a  Fourier  transform  was  performed  over 
the  256  channel  output  of  the  multi-channel  analyzer  or  a 
Fourier  series  expansion  was  made  with  between  50  and  120 
channels  per  cycle. 

The  complex  inverse  relaxation  lengths ,  obtained 
from   least-squares  fits  of  the  amplitude  and  phase  curves, 
were  derived  from  a  careful  examination  of  from  2  to  8  fits 
of  each  of  these  curves  in  which  points  with  largest  error 
were  systematically  eliminated.   The  presentation  of  all 
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of  these  data  would  be  impractical.   One  can,  however, 

present  some  of  the  data  which  are  indicative  of  the 

results  obtained  and  show  sample  calculations  to  indicate 
the  method  of  data  analysis. 

Determina.tion  of  the  Amplitude  and  Phase  of  the 
Fundajnental  xhermal  Wave  of  the  Wave  Data 

For  the  initial  analysis  of  the  wave  data  use  was 
made  of  the  University  of  Florida  Non-Linear  Least  Squares 
(UF-NLLS)  code  C2£).   Several  new  subroutines  were  written 
for  this  code ,  and  it  was  incorporated  into  another  main 
program  to  facilitate  a  Fourier  series  expansion  up  to 
ten  harmonics  of  the  thermal  neutron  wave.   The  bulk  of 
the  numerical  calculations  were  performed  in  the  WAVE- 
ESTPAR  subroutines  which  provided  estimates  of  the  para- 
meters to  the  UF-NLLS  code.   The  WAVE-ESTPAR  subroutine  will 
now  be  discussed. 

Information  was  first  read  into  the  code  which 
allowed  it  to  calculate  the  period  and  number  of  channels 
per  cycle  of  the  fundamental  sine  wave  variation.   The  code 
then  read  in  the  counts  per  channel  recorded  in  the  multi- 
channel analyzer  for  the  standard  detector  with  no  cadmium 
placed  between  the  source  and  the  graphite  pile  (bare 
standard  data) ,  for  the  standard  detector  with  the  cadmium 
sheet  in  place  so  that  only  neutrons  of  epicadmium  energy 
entered  the  graphite  (covered  standard  data),  and  for  the 
movable-detector  counts  per  channel  for  the  latter  situation 
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(covered  movable  data).   The  data  of  interest,  that  arising 
only  from  the  thermal  content  of  the  source,  are  contained 
in  the  movable  detector  counts  per  channel  for  the  bare 
run  (bare  movable  data) ;  but  this  information  is  contaminated 
by  counts  arising  from  the  epicadmium  content  of  the  source. 
The  fundamental  amplitude,  phase,  and  harmonic  content  up 
to  the  tenth  harmonic  of  the  former  three  sets  of  data 
were  then  calculated.   From  the  two  sets  of  standard  data 
the  source  variation  from  the  covered  to  the  uncovered 
runs  for  each  of  the  frequencies  analyzed  could  be  determined. 
This  information  was  then  used  to  generate  correction 
factors  to  be  applied  to  each  of  the  frequency  components 
of  the  covered  movable  data.   A  set  of  data  was  generated 
which  represented  the  contribution  of  epicadmium  neutrons 
from  the  source  to  the  bare  movable  data.   This  was 
accomplished  by  individually  multiplying  the  fundamental 
and  each  of  the  harmonics  of  the  covered  movable  data  by 
the  ratio  of  the  standards  content  of  the  same  frequency  for 
the  covered  and  uncovered  runs  and  by  using  a  subroutine 
NEWCOV  which  synthesized  the  generated  waves.   This  generated 
epicadmium  content  of  the  bare  movable  data  was  then  sub- 
tracted from  it  to  obtain  the  thermal  neutron  wave  arising 
only  from  the  thermal  neutron  content  of  the  source.   Note 
that  this  background  subtraction  has  been  made  individually 
for  each  frequency;  therefore,  changes  in  the  harmonic 
content  of  the  source  from  one  run  to  the  next  should  not 
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affect  the  result.  This  elaborate  method  of  background 
subtraction  was  necessary  because  of  the  instability  of 
the  reutron  generator. 

The  fundamental  component  of  the  thermal  wave  and 
as  many  of  its  first  ten  harmonics  as  desired  were  then 
determined,  and  the  results  were  normalized  to  the  bare 
standard  data.   Corrections  were  also  made  for  changes  in 
the  constant  component  of  the  source  as  determined  from 
the  total  counts  of  the  detectors  for  both  runs.   The  data 
were  normalized  so  that  the  channel  width  of  the  multi- 
channel analyzer,  its  total  number  of  sweeps,  counting 
time,  dead  time  between  channels,  and  channels  per  cycle 
could  be  varied  without  affecting  the  results.   The  print- 
out of  WAVE-ESTPAR  included  the  normalized  amplitude  and 
phase  of  the  frequencies  analyzed  for  the  two  detectors  for 
both  runs  so  that  any  variation  in  the  source  was  recorded. 

The  numerical  methods  used  in  WAVE-ESTPAR  were  rather 
crude  compared  to  those  which  could  have  been  applied.   The 
correction  for  dead  time  between  channels  was  made  by  cal- 
culating the  average  of  the  counts  accumulated  in  two  con- 
secutive channels ,  multiplying  this  by  the  dead  time  between 
channels ,  and  adding  the  result  to  the  counts  of  the  first 
of  these  two  channels.    The  numerical  integrations  associated 


The  212  plug-in-unit  used  for  this  experiment  requires 
10  microseconds  to  store  the  counts  accumulated  in  a 
channel  before  moving  to  the  next  channel.   During  this 
storage  time  no  input  pulses  are  counted. 
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with  the  calculation  of  the  Fourier  coefficients  were 
approximated  by  a  summation.   The  method  can,  however, 
be  justified  for  these  data.   The  operation  of  the  multi- 
channel analyzer  is  such  that  it  will  store  in  each  channel 
all  of  the  counts  accumulated  during  the  count-time  of  that 
channel.   In  effect,  the  counts  recorded  in  a  channel  are 
the  integral  of  the  detector  count  rate  over  the  time  interval 
allotted  to  the  channel.   If  exact  adjustment  could  be  made 
for  the  dead  time  between  channels,  the  summation  of  the 
counts  accumulated  in  all  the  channels  would  be  identical  to 
the  integral  of  the  count  rate  over  the  total  sweep  of  the 
multi-channel  analyzer.   The  weighted  summation  performed  by 
the  code  approximates  the  exact  integral  over  a  period  of 
the  fundamental  sinusoidal  variation  by  a  weighted  summation 
of  integrals  over  each  channel.   In  each  of  these  "channel 
integrals"  the  sinusoidal  function  has  been  effectively  moved 
outside  the  integral,  anc  the  integrated  count  rate  for  the 
time  interval  of  the  channel  is  multiplied  by  the  value  of 
the  sinusoidal  function  at  the  center  of  the  time  interval. 
This  approximation  is  mathematically  expressed  by 

T  N  ^.^^ 

/   sin  nwt  y(t)dt  «   J   sin  nuCt.  +  t/2)   /   y(t)dt 

j=l         ^        tj 

(55) 

where  T  is  the  period  of  the  fundamental  sinusoidal 
variation  of  the  count  rate,  y(t),  t  is  the  channel  width, 
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tj  is  the  time  when  the  j    channel  is  opened;  and  it  has 
been  assumed  that  the  dead  time  was  zero  and  that  there  were 
exactly  N  channels  per  cycle.   The  accuracy  of  the  above 
expression  is  based  primarily  upon  the  value  of  N,  which 
varied  between  50  and  120  for  this  work. 

To  prevent  inaccuracies  in  the  data,  the  harmonic 
distortion,  fast  neutron  content,  and  D.C.  component  of  the 
source  were  all  held  at  as  low  a  value  as  possible;  and 
every  attempt  was  made  to  hold  the  neutron  generator  in 
stable  operation  during  all  of  the  runs  of  a  particular 
frequency.   These  objectives  were  well  realized,  especially 
for  the  second  wave  experiment;  and  for  these  data  the 
elaborate  corrections  seemed  somewhat  extreme.   The  original 
purpose  of  WAVE-ESTPAR  was  to  extract  the  thermal  wave  from 
the  bare  movable  data,  keep  track  of  variations  in  the  source, 
and  provide  estimates  to  the  UF-NLLS  code  of  the  amplitude 
and  phase  of  the  fundamental  thermal  wave.   The  UF-NLLS  code 
was  then  to  be  used  to  modify  these  initial  estimates  to 
values  which  yielded  the  best  fit  to  the  experimental  results. 
Actually,  WAVE-ESTPAR  did  such  a  good  job  of  estimating  the 
parameters  that  changes  in  them  by  the  UF-NLLS  code  were 
made  only  in  the  third  and  fourth  decimal  places.   In  the 
interest  of  saving  computer  time,  the  parameters  calculated 
by  WAVE-ESTPAR  were  accepted  for  the  majority  of  the  data 
considered. 
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To  indicate  the  accuracy  achieved  by  this  technique, 
Figures  16  and  17  show  the  amplitude  and  phase  versus  z 
position  of  the  400  cycle  data  obtained  from  runs  made 
when  the  fundamental  sinusoidal  variation  was  200  cps  and 
then  UOO  cps.   For  the  200  cps  runs,  the  magnitude  of  the 
UOO  cps  distortion  was  about  3  per  cent.   As  shown  in  the 
Figures,  the  UOO  cps  results  obtained  from  this  small  per- 
centage of  the  thermal  source  agree  quite  well  with  that 
for  the  400  cps  data.   In  fact,  plotting  the  first  harmonic 
of  a  particular  frequency  became  a  very  useful  way  to 
determine  if  any  mistakes  had  been  made  in  the  data  input 
since  this  small  component  of  the  measured  signal  was  very 
sensitive  to  inaccuracies  or  inconsistencies. 

Analysis  of  the  Pulse  Data  by  Moore's  Method 

Some  enlightening  derivations  by  Moore  (16,17)  provide 
the  guidelines  for  the  analysis  of  data  obtained  from  the 
pulse  propagation  experiment.   Moore  expressed  Fourier 
moments,  t|)^(q,v)  ,   of  the  measured  data  as 

-*■  •         -2irivt    ^ 

♦j^(q,v)   =   /   dt  t"  e       *(q,t) 

(66) 


*  The  author  is  indebted  to  Dr.  M.  N.  Moore  for  several 
stimulating  discussions  concerning  the  analysis  of 
pulse  propagation  data. 
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Fig.  16.   Amplitude  of  the  400  cps  Wave  as  Obtained  from  the 
First  Harmonic  of  the  200  cps  Wave  Data  and  the 
Fundamental  of  the  400  cps  Wave  Data. 
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Fig.  17.   Phase  of  the  400  cps  Wave  as  Obtained  from 
the  First  Harmonic  of  the  200  cps  Wave  Data 
and  the  Fundamental  of  the  400  cps  Wave  Data, 
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where  ^(q,t)  is  the  experimentally  measured  pulse  at  space 
point  q  as  a  function  of  time,  t,  and  v  is  the  frequency  in 
cycles  per  second.  He   then  defined  Fourier  widths,  W  (a.v). 
by  the  expression 

W^Cq.v)   =   C  /"  dt(t-t)^  e        if»(q,t)3/4'Q(q.v) 

(67) 

where 

t   =   *i(q,\')/^Q(q,v)  (58) 

and  proved  that,  neglecting  higher  space  and  energy  mode 
contamination,  these  widths  are  related  to  the  frequency 
derivatives  of  the  system  dispersion  law  by  the  expression 


(69) 


+  [-a,(v,B2)  +  iC.(v,B2)]z  }    n  >  1 

dv^      1*1      1   ♦  i.    ^ 

where 

a" 

^n  =   T;;^  l°g  A;L0Q23^q2»q3»Bi)  (70) 

and 
<|'Q(q,v)  =  A]_Q(v)Q23(q2.q3»B_[>exp  (  [-03^(  v,  B^^)  +  i^^^  (  v  ,  BJ_)  ]z  } 
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In  the  above  expression,  'I'-Cqiv)  is  the  fundamental  space 
and  energy  mode  which  propagates  through  the  system.  The 
source  condition  is  specified  through  A,q,  and  the  trans- 
verse space  dependence  is  expressed  by  Qot  Wo  »*^3  »^J_^ '  A-^^° 
ai  and  C]_  are,  respectively,  the  real  and  imaginary  parts 
of  the  fundamental  complex  inverse  relaxation  length. 

The  significance  of  these  expressions  is  that  they 
form  the  link  between  the  pulse  propagation  experiment  per- 
formed in  the  time  domain  and  the  theory  derived  in  the  fre- 
quency domain.   If,  for  example,  a  unit  impulse  were  intro- 
duced at  z  =  0,  ijj(q,t)  would  be  the  unit  impulse  response  whose 
Fourier  transform,  "iioWj^)*  is  the  transfer  function  of  the 
system.   Herein  lies  the  key  to  the  normalization  of  data 
from  the  pulse  propagation  experiment  and  from  the  wave 
experiment.   These  expressions  indicate  that  the  amplitude 
and  phase  curves  measured  directly  by  the  wave  experiment 
can  also  be  obtained  by  numerically  Fourier  transforming 
the  data  obtained  from  the  pulse  propagation  experiment. 
However,  much  more  than  this  is  evident  from  a  careful 
inspection  of  the  equations.   Not  only  are  a  and  C  readily 
obtainable  from  the  raw  data  of  the  pulse  propagation 
experiment,  but  their  derivatives  with  respect  to  frequency 
can  also  be  calculated  by  the  use  of  t  and  the  Fourier 
widths.   Expressing  t   in  the  form 
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t(v,B]^)  =  ~  ^  Gi  +  f-  [-a.(v.B2)  +  iCi(v.B2)]z  } 

(71) 

indicates  that  plots  of  the  imaginary  and  real  parts  of 
t  versus  z,  numerically  evaluated  for  a  particular  frequency 
by  the  use  of  equations  (66)  and  (58),  should  yield  linear 
curves  with  slopes  equal  to  a'  and  C',  respectively.   The 
same  calculation  performed  at  several  frequencies  yields 
the  frequency  dependence  of  the  first  derivative  of  the 
complex  inverse  relaxation  length.   Similarly,  an  inspec- 
tion of  equation  (69)  shows  that  the  Fourier  width,  VJ 
calculated  numerically  from  equation  (67)  can  be  used  to 
determine  experimentally  the  n   derivative  of  a  and  C  as 
a  function  of  frequency.   Thus  one  cannot  only  compare  the 
theoretically  derived  a  and  Z    in  the  frequency  domain  to 
those  obtained  in  the  experiment;  but  one  can  also  perform 
the  more  demanding  comparison  of  the  theoretically  and 
experimentally  determined  derivatives  of  a  and  E, . 

Moore  has  also  shown  that  this  comparison  is 
especially  applicable  when  the  perturbation  expansion  method 
is  used  to  determine  the  complex  eigenvalue  p^.   In  fact, 
a  direct  determination  of  the  coefficients  of  this  expansion 
is  possible  through  application  of  his  suggested  procedure 
for  data  analysis.   This  is  easily  illustrated  by  considering 
the  p^  expansion,  which  for  convenience  is  expressed  in  the 
form 
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N    B^"^c."         ^*1   B^"^«" 


p2   (^)   =      [      ; +   i        I  ; 

T        '^     n!  ^     n! 

n=0  n=l 

(even)  (odd)  (72) 

where  B  is,  as  before,  the  n^"  coefficient  of  the  power 
series  expansion;  and  the  T  subscript  of  p^  indicates  that 
this  is  a  theoretical  prediction.   Since 

p   s   a   ■»•   iC  (73) 

and  a  and  ^  are  what  one  obtains  experimentally,  the 
contrast  between  theory  and  experiment  is  expressed  by 

nCn)  n  o(n)  n 

2         2     ..2.'>-  r^U..     rDO) 

n=0  n=l 

(even)  (odd) 


(74) 

The  double  arrow  above  indicates  reciprocal 
correspondence,  and  the  E  subscript  stands  for  experimental 
quantities.   Thus  a  power  series  fit  in  the  frequency 
variables   a^   -  ^^  and  20^?^  yields  coefficients  which 
should  be  equal  to  the  theoretically  determined  coefficients 
of  the  perturbation  expansion.   If  o'  and  5*  were  also  cal- 
culated,  one  could  contrast  theory  and  experiment  through 
the  relationship 
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i;   'E       =   2(a.a^  -  ^^^^)    .    2i(a'5^  *    C'a^)   - 


n=2      ^^  n=l       ^- 

(even)  (odd) 

(75) 

Allowing  0)  fo  approach  zero  in  the  above  relationship 
indicates  that  the  zero  intercept  of  the  experimental 
curve  2(0^5^;  +  ^E**!^  yields  a  value  to  be  compared  with 
the  theoretically  determined  B^-^  .   Likewise,  if  the  second 
derivatives  of  o  and  5  could  be  determined,  the  experimen- 
tal curve  arising  from  the  derivative  of  the  real  part 

of  equation  (75)  would,  when  extrapolated  to  w=0,  have 

(2) 
a  value  to  be  compared  to  B    of  the  theory.   Moreover, 

one  could  also  fit  equation  (75),  and  others  derived  from  it 
using  higher  derivatives  of  a  and  £; ,  to  power  series  expan- 
sions in  still  another  attempt  to  determine  these  coefficients. 
The  above  discussion  illustrates  that  the  analysis  of  pulse 
propagation  data  by  the  method  suggested  by  Moore  can  be 
visualized  as  the  experimental  equivalent  to  the  perturbation 
expansion  developed  by  Perez, 

Moore  has  also  derived  a  very  elegant  method  by  which 
the  derivatives  of  a  and  C  can  be  calculated.   As  indicated 
previously,  one  calculates  the  Fourier  width,  W  ,  of  the 
thermal  pulse  at  several  z  positions  and  fits  their  real 
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and  imaginary  components  to  a  linear  function  of  z.   From 
equation  (59)  it  is  evident  that  the  slopes  of  these  curves 
yield  a   and  E,^ ,      One  does  not  have  to  take  pulse  propagation 
data  to  be  able  to  determine  these  derivatives ,  however. 
If  only  wave  data  were  available,  one  could  arrive  at  a 
set  of  points  on  the  curves  of  a  and  C  versus  frequency  and 
then  apply  numerical  differentiation  techniques  for  the 
determination  of  the  derivatives.   This  latter  method, 
though  mathematically  equivalent  to  that  suggested  by  Moore, 
is  quite  different   experimentally.   Using  Fourier  widths, 
the  derivatives  are  obtained  almost  directly  from  the  raw 
data,  that  is,  from  iiiq^t).      For  the  other  method,  a'  and 
i\   for  example,  are  obtained  through  the  intermediate  step 
of  calculating  a  and  ^;  and  the  errors  accumulated  in  the 
latter  two  will  surely  propagate  to  their  derivatives. 

One  further  comment  concerning  the  Fourier  widths 
will  be  made  before  proceeding  to  the  discussion  of  the 
computer  code  used  to  analyze  the  pulse  propagation  data. 
An  inspection  of  equation  (67)  reveals  that  W   =  1  and  W,  =  0. 
These  relationships  can  be  used  to  check  the  correctness 
of  the  method  of  data  analysis,  the  accuracy  of  the 
numerical  methods,  and  the  consistency  of  the  data. 

The  MOORE  MOMENTS  code,  though  complicated  by  the 
fact  that  the  integrands  contain  complex  numbers,  is  simple 
in  terms  of  its  overall  sequence  of  operations.   One  reason 
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for  this  was  that  the  subtraction  of  the  signal  caused  by 
source  neutrons  with  energies  above  cadmium  cutoff  was  able 
to  be  based  only  upon  the  total  counts  accumulated  in  the 
standard  detector  during  both  the  covered  and  the  uncovered 
runs.   As  shown  by  Hartley  (14).  the  fast  pulse  was  localized 
in  the  first  few  channels  of  the  multi-channel  analyzer  while 
the  thermal  pulse  extended  over  most  of  the  255  channels. 
Thus  extraction  of  the  thermal  pulse  from  the  bare  movable 
data  by  taking  the  source  corrected  difference  between  it 
and  the  covered  movable  data  was  accepted  as  adequate  for 


these  data. 


For  the  numerical  approximation  to  equations  (66)  and 
(67),  the  same  method  was  used  as  has  been  explained  in 
connection  with  the  WAVE-ESTPAR  subroutine.   Some  difficulty 
arose  in  justifying  the  fact  that  theory  called  for  an 
integral  over  an  infinite  range  and  the  data  extended  over 
a  finite  interval  of  time.   However.  Hartley's  exacting 
experimental  technique  reduced  this  problem  to  an  academic 
one.   He  accumulated  the  data  on  one  long  run  of  about  30 
hours  during  which  time  the  neutron  generator  was  held  within 
narrow  operating  limits.   The  ratio  of  maximum  to  minimum 
value  of  the  thermal  pulse  was  maintained  well  above  1000  to  1 
and  showed  only  a  slight  drift  over  the  course  of  the  experi- 
ment.  Enough  channels  were  allotted  to  times  after  the 
complete  decay  of  the  pulse  so  that  the  D.C.  contribution. 
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Which  was  due  to  the  constant  drift  current,  could  be 
determined.   Thus  the  contributions  to  the  integral  beyond 
the  counting  range  on  the  multi-channel  analyzer  was  held 
to  a  negligible  amount  and  could  be  subtracted  from  the 
thermal  pulse  before  performing  the  numerical  integrations. 

The  MOORE  MOMENTS  code  was  set  up  to  evaluate 
numerically   for  any  desired  frequency   the  Fourier  moments, 
equation  (66) ■  t,  equation  (68);  and  Fourier  widths  for  any 
value  of  n,  equation  (67).   Corrections  were  made  for  the 
dead  time  between  channels  and  the  drift  current  contribu-  • 
tion.   The  numerical  integration  scheme  could  be  either  a 
direct  summation  or  Simpson's  rule.   The  results  were  normal- 
ized, printed  out,  and  plotted.   Little  difference  (less  than 
1  Der  cent)  was  obtained  when  Simpson's  rule  was  applied 
rather  than  the  direct  summation.   This  was  because  256 
points  were  available  in  the  numerical  integration.   This 
difference  increased  with  frequency  but  was  never  more  than 
1  per  cent  for  the  range  of  frequency  considered.   The 
value  of  W^,  which  should  be  zero,  was  quite  dependent  upon 
the  value  of  the  argument  selected  for  the  sinusoidal  functions 
used  in  the  formation  of  the  integrand.   It  was  found  that 
the  smallest  W^  values  were  obtained  when  these  arguments 
corresponded  to  times  at  the  center  of  the  counting  interval 
of  the  channels.   Values  of  about  10-12  ^^^^  obtained  at 
the  low  frequencies,  but  they  increased  to  about  10"^  at  the 
high  frequencies. 
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The  Fit  of  the  Combined  Amplitude  and  Phase  Curves 
to  Obtain  the  Real  and  Imaginary  Components  of  the 
Complex  Inverse  Relaxation  LengxTT 

As  previously  mentioned,  the  curves  of  amplitude 
and  phase  versus  z  position  for  a  particular  frequency  are 
obtained  from  the  Fourier  series  analysis  of  the  wave  data 
and  from  the  zero  moment,  li'oCqjv)  of  equation  (56),  of  the 
pulse  propagation  data.   The  slope  of  the  semilog  plot  of 
the  amplitude  curves  is  a,  and  the  slope  of  the  linear 
phase  curve  is  C.   To  obtain  these  components  of  the 
complex  inverse  relaxation  length,  the  amplitude  and  phase 
curves  from  both  wave  experiment  and  from  the  pulse  propa- 
gation experiment  were  combined;  and  the  UF-NLLS  code  was 
used  to  obtain  the  slope  of  the  result.   It  would  be 
impractical  to  present  plots  of  each  of  the  130  amplitude 
and  phase  curves.   However,  all  of  the  data  points  of  these 
curves  have  been  tabulated  in  the  Appendix,  and  several  of 
them  are  shown  in  Figure  13  and  Figure  19  to  indicate  the 
agreement  obtained  between  the  wave  data  and  pulse  propaga- 
tion data.   One  point  which  is  clear  from  an  inspection  of 
these  curves  is  that  data  were  able  to  be  obtained  at 
larger  z  values  with  the  pulse  propagation  experiment  than 
with  the  wave  experiment.   This  is  because  in  the  wave 
experiment  the  D.C.  component  of  the  source  must  be  at 
least  equal  in  magnitude  to  the  sinusoidal  component.   The 
D.C.  component  suffers  less  attenuation  with  distance  from 
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Fig.  19.   Combined  Phase  Curves  for  Several  of  the  Frequencies. 
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the  source.   After  a  decade  or  so  of  decay  of  the  wave,  its 
contribution  to  the  movable  detector  count  rate  is  swamped 
by  the  D.C.  component.   In  this  circumstance,  cross  correla- 
tion techniques  could  be  applied  to  extract  the  sinusoidal 
component  of  the  signal  but  this  was  not  attempted.   For 
the  pulse  propagation  experiment,  on  the  other  hand,  the 
D.C.  component  of  the  source  is  determined  only  by  the 
on-off  ratio  of  the  neutron  generator,  which  for  this  data 
was  about  1000  to  1. 

The  amplitude  curves  were  normalized  to  a  value  of 
unity  at  the  first  z  position  at  which  data  were  obtained, 
and  their  logarithms  were  least-squares  fit  to  a  linear 
function  of  z  with  a  slope  of  a  and  a  zero  intercept.   The 
phase  curves  were  normalized  to  a  value  of  zero  at  the 
closest  z  position,  and  they  were  least-squares  fit  to  a 
straight  line  with  slope  5  and  zero  intercept.   Only  one 
parameter,  the  slope,  was  varied  through  the  iterative 
process  of  the  code.   It  is  important  to  note  that  for 
the  calculation  of  a  and  C  only  the  relative  amplitude  and 
phase  of  the  neutron  wave  at  various  z  positions  were  re- 
quired.  Thus  the  frequency  spectrum  of  the  source  did  not 
need  to  be  determined.   The  interest  here  was  how  the  wave 
of  a  particular  frequency  was  attenuated  with  distance  from 
the  source,  not  with  how  its  amplitude  and  phase  compared 
with  that  of  another  frequency  at  a  certain  point  in  the 
medium.   The  ability  to  adequately  analyze  the  data  without 
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knowledge  of  the  frequency  dependency  of  the  source  was 
found  to  be  very  fortunate,  especially  when  an  attempt  was 
made  to  obtain  this  knowledge  in  connection  with  the  re- 
flection experiment  (13). 

The  final  determination  of  a  and  C  resulted  from 
several  fits  of  each  of  the  amplitude  and  phase  curves  in 
which  points  were  systematically  eliminated  from  both  their 
high  z  values  and  low  z  values.   The  points  for  large  z 
were  poor  because  of  statistical  errors.   The  thermal  pulse 
was  detected  over  an  attenuation  of  nearly  3  decades ,  but 
for  the  same  range  of  z  values  the  pulse  resulting  from  the 
epicadmium  energy  source  neutrons  was  only  slightly  attenuated, 
Thus  the  ratio  of  useful  signal  to  background  decreased 
rapidly  with  increasing  z  because  of  the  rapid  attenuation 
of  the  useful  signal.   The  points  obtained  at  small  z  values 
were  contaminated  by  space  and  energy  modes  and  were,  there- 
fore, not  as  useful.   The  first  fit  of  the  data  included  all 
of  the  points.   The  channel-chop  option  of  the  UF-NLLS  code 
was  then  used  to  eliminate  the  smaller  z  data  in  sets  of  3 
to  8  points  at  a  time,  depending  on  the  total  number  of 
points  and  their  z  values.   For  this  first  sequence  of 
channel  chops ,  the  points  fit  covered  the  approximate  ranges 
of  4  through  69  centimeters,  10  through  69  centimeters,  15 
through  59  centimeters,  and  20  through  59  centimeters.   The 
highest  z  points,  those  with  largest  error,  were  then 
eliminated  from  each  curve.   The  number  of  points  eliminated 
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increased  slowlv  with  ^requencv  from  zero  to  3  or  'i  at 
around  1000  cps ,  and  then  more  rapidly  to  6  to  8  at  the 
high  frequency  limit  of  IHSO  cps.   The  channel-chop 
sequence  was  then  applied  to  these  depleted  curves  in 
the  same  fashion  as  before. 

The  results  obtained  frorr.  this  procedure  were  verv 
gratifying.   At  the  lower  frequencies  (10  to  300  cps),  the 
variation  of  the  slopes  with  the  number  of  points  chopped 
clearly  approached  an  asymptotic  value  in  a  fashion  which 
indicated  the  progressive  elimination  of  contributions  from 
higher  space  and  energy  modes.   For  this  frequency  range, 
values  of  the  slope  versus  the  number  of  points  chopped 
were  plotted  and  a  smooth  curve  drawn  through  them.   All 
of  these  curves  tended  towards  an  eisymptotic  value  which 
was  taken  to  be  that  of  the  fundamental  space  and  energy 
mode.   The  change  in  slope  observed  over  these  curves  was 
about  5  per  cent.   For  the  rest  of  the  frequency  range, 
no  such  trend  was  observed;  and  up  to  about  1000  cps,  the 
final  slopes  were  obtained  from  a  simple  averaging  of  the 
various  values  obtained.   The  percentage  variation  of  these 
slopes  was  slightly  less  than  for  the  lower  frequencies 
except  between  the  range  of  about  700  cps  to  900  cps  where 
the  points  appeared  to  exhibit  more  scatter.   Only  one  or 
two  chops  were  made  and  fewer  points  were  available  above 
1000  cps  so  the  final  values  of  these  slopes  are  less  accurate, 
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With  the  attainment  of  a   and  c  from  the  amplitude 
and  phase  curves,  the  first  stage  of  the  data  analysis  is 
complete.   It  is  hoped  that  these  experimentally  determined 
values  are,  in  fact,  the  real  and  imaginary  parts  of  the 
fundamental  complex  inverse  relaxation  length.   The 
presentation  of  a  and  C  versus  frequency  is  left  to  the 
next  chapter. 


CHAPTER  V 

COMPARISON  OF  THE  THEORY  WITH  THE  EXPERIMENT 

Presentation  of  the  Measured 

Complex  Inverse  Relaxation 

Length 

It  has  been  proposed  by  Beckurts  (29)  that  it  is 
more  useful  to  contrast  the  experimentally  obtained  eigen- 
values with  the  ones  computed  theoretically  than  to  compare 
parameters  extracted  from  the  data  to  their  values  as  pre- 
dicted by  the  theory.   Consistent  with  this  idea,  Fipures 
20  and  21  show  the  real  and  imaginary  components  of  the 
experimentally  determined  complex  inverse  relaxation 
length  versus  frequency,  along  with  those  curves  from  the 
theory  using  both  the  SUMMIT  kernel  and  the  MASSMK  kernel. 
From  an  inspection  of  these  figures,  it  appears  that  both 
theoretical  models  predict  the  experiment  quite  well  up  to 
about  300  cps .   For  higher  frequencies,  the  theory  and 
experiment  diverge  slightly,  with  the  discrepancy  of  the 
MASSMK  kernel  exceeding  that  of  the  SUMMIT  kernel  by  a 
noticeable  amount.   In  the  frequency  range  between  400 
cps  and  1000  cps,  the  deviation  between  the  SUMMIT  kernel 
and  the  experiment  rises  from  4  per  cent  to  about  10  per  cent, 
while  in  this  same  frequency  range,  the  MASSMK  deviation 
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ranges  between  8  and  15  per  cent.   Both  theories  tend 
to  underestimate  the  damping  coefficient,  a,  and  to  over- 
estimate the  phase  shift  factor,  ^,  especially  at  higher 
frequencies . 

One  could  conclude  from  this  direct  comparison  that 
the  disagreement  between  theory  and  experiment  exists  only 
for  frequencies  higher  than  300  cps .   However,  an  ambiguity 
exists  in  that  this  disagreement  can  be  attributed  to  one 
or  more  of  the  following  facts:   (1)  higher  energy  mode 
contamination,   (2)  higher  spatial  mode  contamination, 
(3)  transport  effects,   (4)  inadequacy  of  the  particular 
kernels  which  were  used  to  describe  the  thermalization 
effects.   Consideration  will  now  be  given  to  each  of 
these  possibilities  in  the  same  order  as  they  were  listed. 

Both  scattering  kernels  were  used  to  generate  the 
complex  inverse  relaxation  lengths  associated  with  higher 
energy  modes.   It  was  found  (see  Figure  6)  that  the  damping 
coefficients  for  these  higher  energy  modes  are  significantly 
larger  for  low  frequencies  but  become,  in  general,  smaller 
than  the  fundamental  a  beyond  the  first  crossing  point  of 
300  cps.   This  result  agrees  with  the  observation  by  Starr 
and  Price  (11_)  that  for  large  bucklings  waiting  times  on 
the  order  of  two  milliseconds  were  needed  to  obtain  the 
asymptotic  behavior  in  die-away  experiments  in  graphite. 
In  terms  of  the  neutron  wave  experiment,  this  is  equivalent 
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to  the  prediction  that  higher  energy  modes  will  become 
prevalent  for  frequencies  larger  than  500  cps. 

As  frequency  increases,  all  of  the  spatial  modes 
associated  with  the  fundamental  energy  mode  approach  the 
same  values  for  their  damping  and  phase  shift  coefficients. 
This  is  because  they  differ  only  by  the  transverse  leakage 

term,  B     ^  which  becomes  progressively  less  significant 
-»-,  I  ,m 

compared  to  the  frequency  dependent  terms  of  p^  as  the 
frequency  increases.   At  low  frequencies  these  higher 
spatial  modes  can  be  easily  discriminated  against  by 
obtaining  data  only  at  points  far  from  the  source  boundary 
and  by  applying  the  channel-chop  technique  previously 
discussed.   At  higher  frequencies,  however,  one  cannot  be 
sure  that  higher  spatial  modes  are  not  contributing  to  the 
data  since  their  attenuation  and  phase  shift  factors  will 
be  nearly  equal  to  those  of  the  fundamental  spatial  mode. 
Perez  and  Uhrig  (9)  estimated  that  for  a  point  source  in 
the  assembly  whose  transverse  buckling  is  1.65  x  10    cra~^ 
higher  spatial  mode  contamination  could  become  as  large 
as  5  per  cent  in  the  upper  frequency  range. 

Another  source  of  difficulty  arises  from  the  fact 
that  transport  theory  should  be  introduced  when  the  time 
between  two  collisions  becomes  of  the  same  order  of  mag- 
nitude as  the  inverse  of  the  frequency.   This  effect  is  less 
important  than  the  thermalization  effects,  and  it  has  been 
estimated  by  Kunaish  (1_8)  to  be  on  the  order  of  2  per  cent. 


104 

Hence  the  direct  comparison  test  without  additional 
information  can  both  preclude  any  firm  conclusion  about  the 
disagreement  between  theory  and  experiment  at  the  higher 
frequencies  and  lead  one  to  assume  that  since  there  was 
agreement  over  a  reasonable  range  of  frequencies  the 
diffusion  and  thermalization  parameters  derived  from  the 
theory  are  adequate.   From  this  fact  the  question  arises 
as  to  how  sensitive  the  direct  comparison  test  is  to 
thermalization  phenomena.   Still  another  example  of  the 
difficulty  in  deriving  firm  conclusions  concerning  thermali- 
zation effects  by  the  direct  comparison  of  eigenvalue  curves 
can  be  obtained  by  applying  the  reasoning  presented  thus 
far  to  typical  die-away  experiments.   By  doing  this,  one 
could  conclude  that  diffusion  theory  gives  adequate  theore- 
tical predictions  except  in  the  range  of  bucklings  where 
the  experimental  data  can  be  affected  by  higher  energy 
and  spatial  modes  as  well  as  by  transport  corrections.   The 
inadequacy  of  the  direct  comparison  test  is  confirmed  by  the 
fact  that  the  MASSMK  kernel,  which  is  much  less  sophisticated 
from  a  thermalization  viewpoint  than  the  SUMMIT  kernel, 
gives  results  which  appear  very  similar  to  those  of  the 
SUMMIT  kernel. 

The  difficulties  encountered  in  the  direct  comparison 
test  for  both  the  die-away  and  the  neutron  wave  experiments 
arise  because  the  thermalization  parauneters  enter  as  second- 
order  effects  in  the  frequency  range  in  which  the  asymptotic 
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mode  is  more  easily  measured.   In  order  to  perform  a 
meaningful  comparison  between  the  theory  and  the  experi- 
mental results,  one  should  be  able  to  isolate  thermalization 
effects  from  diffusion  effects. 

In  the  conventional  die-away  experiment  this  is  not 
possible  because  only  one  quantity,  the  decay  constant,  is 
measured.   In  the  neutron  wave  experiment,  however,  two 
independent  parameters,  a   and  ^,  are  measured;  and  they 
can  be  combined  to  yield  two  frequency  dependent  functions, 
2a C  and  a^  -  C^.   The  former  is  primarily  dependent  upon 
diffusion  parameters;  and  the  latter,  except  for  the  zero- 
frequency  intercept,  is  exclusively  a  function  of  thermali- 
zation parameters. 

Figure  22  shows  a  plot  of  2at  as  a  function  of  fre- 
quency.  Since  the  fundamental  dependence  of  this  function 
is  upon  diffusion  parameters,  the  agreement  between  the 
SUMMIT  theory  and  the  experimental  results  appears  within 
the  experimental  error  over  the  full  frequency  range  con- 
sidered.  The  MASSMK  kernel  also  represents  the  data  well 
in  a  reasonably  large  range  of  frequencies.   Note  that  the 
above  is  especially  true  for  frequencies  less  than  300  cps , 
where  no  problems  associated  with  modal  contamination  and 
transport  corrections  are  expected  to  occur. 

The  a^  -  C ^  versus  u^  function,  shown  in  Figure  23  for 
the  frequency  range  from  10  cps  to  1400  cps  and  in  Figure  24 
for  the  frequency  range  below  UOO  cps,  indicates  that  sizeable 
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differences  exist  among  the  MA5SMK  kernel,  SUMMIT  kernel, 
and  the  experimental  results  over  the  full  frequency  range. 
Note  that  the  one-velocity  diffusion  theory  model  predicts 
a  constant  behavior  of  a^  -  C^.   This  is  the  only  curve 
thus  far  presented  in  which  noticeable  differences  exist 
between  the  two  theories  and  the  experimental  results  in 
the  zero  to  300  cps  frequency  range.   Therefore,  this  is 
the  first  instance  in  which  the  various  thermalization 
models  can  be  compared  with  the  experimental  results  in 
the  frequency  interval  where  the  experimental  data  are 
most  accurate  and  reliable  and  with  no  interference  arising 
from  diffusion  phenomena.   In  other  words,  the  small 
differences  in  the  theoretically  and  experimentally 
determined  a    and  C,  which  are  mainly  due  to  thermalization 
effects,  tend  to  compensate  in  the  2a5  curve  and  to  be 
amplified  in  the  a^  -  ^^  curve. 

An  inspection  of  Figures  23  and  2U  clearly  indicates 
that  although  the  trend  of  the  SUMMIT  theory  correctly 
describes  that  of  the  experimental  results  the  coefficient 
associated  with  the  linear  rise  at  low  frequency  will  be 
about  twice  as  large  in  the  experimental  curve  as  the  one 
for  the  theory.   The  evaluation  of  this  coefficient  and 
the  others  to  be  compared  with  the  theory  are  discussed 
in  the  next  section. 
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Least-Squares  Fits  of  the  Squared  Complex 
Inverse  Relaxation  Length 

The  real  and  imaginary  components  of  p^,  a^-^^ ^ 
and  2aC,  respectively,  were  least-squares  fit  to  poly- 
nomials of  the  form 

N 

n  =  0 
(EVEN)  (76) 

U-1 

n=l 
(ODD)  (77) 

SO  that  the  P    thereby  determined  could  be  compared  with 
the  coefficients,  B^'^Vn!,  obtained  for  the  perturbation 
expansion  of  the  theory.   A  STEP-WISE  MULTIPLE  REGRESSION 
coae   was  used  for  these  least-squares  fits.   Each  data 
point  was  weighted  by  the  inverse  of  its  variance  as 
determined  by  the  error  analysis.   The  code  added  one  term 
of  the  above  expressions  at  a  time  in  the  order  of  their 
importance  to  the  "goodness  of  fit"  of  the  power  series  to 
the  experimental  data.   Terms  were  added  as  long  as,  and 
only  if,  their  inclusion  increased  the  "goodness  of  fit" 
by  a  statistically  significant  amount  as  determined  by 
both  the  F-test  and  the  specific  level  of  significance 
specified  by  the  input.   By  analyzing  the  data  in  this  way. 


*   This  code  was  obtained  from  the  University  of  Florida 
Computer  Center. 


Ill 

the  answer  to  the  difficult  question  of  how  many  terms 
to  include  in  the  expansions  was  based  upon  dependable, 
statistical  tests  and  not  left  to  intuition. 

It  was  found  that  for  the  2aC  curve  experimentally 
determined  in  this  work  only  the  p^^^  ^^^  p(3)  coefficients 
were  statistically  significant  enough  to  be  evaluated.   In 
the  a^  -  C^  curve  the  P^'^^  P^2)^  ^nd  P^^^  coefficients 
were  important  enough  to  be  determined.   As  a  test  of  the 
accuracy  of  this  method  for  the  extraction  of  the  P^'^\ 
the  highest  frequency  points  were  progressively  chopped 
from  the  data,  5  at  a  time;  and  the  depleted  data  re- 
analyzed.  The  results  of  this  procedure  are  shown  in 
Figure  25  along  with  the  theoretical  predictions  for 
their  values  as  derived  from  the  perturbation  expansion. 
This  very  small  variation  of  these  P^"^  can  be  interpreted 
in  two  ways: 

(1)  No  space  or  energy  mode  contamination  was 
present  in  the  experimentally  determined  a  and  C 

(2)  Space  and  energy  mode  contamination  were 
present  but  combined  in  such  a  way  that  cancellation  of 
their  contamination  occurred. 

Although  one  would  like  for  the  actual  situation 
to  be  correctly  described  by  possibility  (1),  above,  there 
is  strong  theoretical  evidence,  previously  presented,  to 
indicate  the  contrary.   More  realistically,  one  should 
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expect  that  some  modal  contamination  will  exist  in  any 
neutron  wave  experiment.   The  main  point  of  concern, 
therefore,  is  focused  upon  whether  this  contamination 
affected  the  final  results. 

Kunaish  (1^)  has  shown  that  for  the  heavy-gas 
kernel  theoretical  calculations  confirm  the  fact  that  the 
higher  spatial  and  energy  modes  can  combine  in  such  a  way 
that  the  total  flux  propagates  with  a  complex  inverse 
relaxation   length  which  is  nearly  equal  to  the  one 
associated  with  the  fundamental  complex  eigenvalue.   Also, 
if  higher  mode  contamination  were  present,  the  p^^^  would 
be  frequency  dependent  unless  cancellation  occurred.   The 
above  considerations  have  led  the  author  to  postulate  that 
modal  contamination  was  present  in  the  measured  data;  but 
enough  cancellation  occurred  so  that  it  did  not  affect, 
or  would  not  be  detected  by,  the  data  analysis  procedure. 

From  the  information  presented  in  Figure  25  and 
similar  considerations  associated  with  P^^s  the  best 
estimates  for  the  values  of  the  several  p'"^  determined 


are 


P^°^  =  (63.86  T  .09)  X  lO-"^  cm-2 

P^l^  =  (4.64  +  .01)  X  10-6  gee  cm-2 

P^^^  =  (3.92  ;  .09)  X  10-10  sec2cm-2 

P^^^  =  -(1.02  T  .07)  X  lO-l**  sec3cm-2 

P^"^^  =  -(11.   T  1.   )  X  10-18  sec'+cm-2 
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Application  of  Moore's  Method  to  the 
Analysis  of  the  Pulse 
Pr'opagation  Data 


The  MOORE  MOMENTS  code  was  used  to  determine  the 
Fourier  moments  and  widths  of  the  pulse  propagation  data. 
Shown  in  Figures  2  6  and  2  7  are  the  real  and  imaginary 
components  of  t"  versus  z  for  several  of  the  frequencies 
analyzed.   Consistent  with  Moore's  prediction,  these 
curves  vary  linearly  with  distance  from  the  source.   The 
slopes  of  these  straight  Ixnes ,  C    and  a',  were  determined 
and  are  shown  in  Figure  28.   Also  shown  are  values  of  the 
derivatives  as  determined  from  numerical  differentiation, 
using  six  order  differences,  of  the  a  and  E,    curves.   It 
appears  that  at  these  lower  frequencies  the  Fourier  moments 
allow  for  a  more  accurate  determination  of  the  derivatives. 
From  the  four  quantities  a,  C,  a',  and  C',  one  can  construct 
the  first  frequency  derivatives  of  a^-^^    and  2aC.   The 
former  should  increase  linearly  with  frequency  with  a 
slope  of  2P^2)  until  the  P^^^  term  becomes  significant,  and 
the  latter  should  be  a  constant  at  low  frequencies  with  the 
value  P^-^^.   These  quantities  are  shown  in  Figure  29,  where 
it  is  evident  that  success  in  the  determination  of  the  first 
derivative  of  p^  was  limited  to  its  imaginary  component. 
The  zero  intercept  of  this  curve  yields  a  P^^^    value  of 


*   The  prime  as  used  here  indicates  derivative  with 
respect  to  oj  . 
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Fig.  26.   The  Experimentally  Determined  Real  Part  of  t   versus 
Distance  from  the  Source. 
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Fig.  27.   The  Experimentally  Determined  Imaginary  Part  of  t 
versus  Distance  from  the  Source. 
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Fig.  28.   The  Experimentally  Determined  a'  amd  C' 
as  a  Function  of  Frequency. 
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Fig,  29.   First  Frequency  Derivatives  of  2aC  and  a^  -  C^ 
versus  u.' 
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(4.65  1  ,02)  X  10"^  sec  cm~^,  which  is  in  good  agreement 
with  the  value  obtained  from  the  power  series  fit  of  2ciC. 
As  in  the  case  of  the  power  series  analysis,  the  curve 
from  which  P^ -^ ''  is  calculated,  2(aa'   -C5'),  results 
from  the  difference  of  two  nearly  equal  quantities.   The 
error  propagation  associated  with  this  calculation  was 
more  severe  than  the  ones  encountered  in  the  determination 
of  P^l^  from  the  curve  2(a'C  +  5'a).   An  inspection  of 
the  relative  magnitude  of  the  parameters  indicated  that 
at  these  low  frequencies  a^  >>  Cq  and  C'  >>  a'.   Thus 
the  (2aC)'  curve  is  dominated  by  the  termaC,  both  of 
whose  components  are  nearly  constant  in  the  frequency 
range  considered.   On  the  other  hand,  both  the  terms  aa' 
and  CC'  from  which  the  (a^-C^)'  curve  is  formed  are 
significant;  and  they  contain  terms  which  vary  comparatively 
rapidly  with  frequency.   It  appears  that  difficulties  are 
self-compensating  for  the  p^-'-^  calculation  and  amplified 
when  p^2)  is  to  be  determined.   Although  this  curve  of 
(a^-C^)'  shows  the  correct  trend  and  agrees  with  reasonable 
values  for  p*- '^ -' ^  it  cannot  be  used  for  an  accurate  deter- 
mination of  this  parameter.   Similar  results  were  obtained 
when  the  second  derivative  of  p^  (i.e.,  (a^-^^)")  was 
calculated  from  a,  a',  a",  ^,  ^'  and  C".   These  second 
derivatives,  a"  and  C",  were  determined  from  the  linear 
variation  of  W2  with  frequency.   The  zero  intercept  of  the 
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experimentally  determined  pj^  yielded  a  value  only  within 
a  factor  of  2  of  the  expected  value,  with  an  error  roughly 
equal  to  its  magnitude. 

The  results  of  the  analysis,  however,  appear  en- 
couraging.  The  experimental  results  confirm  the  validity 
of  Moore's  method  of  data  analysis.   It  was  found  that  the 
accuracy  of  the  higher  moments  is  strongly  dependent  upon 
the  number  of  time  channels  used  to  describe  the  pulse  and 
its  foreground  to  background  ratio.   In  this  experiment, 
channel  widths  of  30  microseconds  yielding  2  56  channels 
per  pulse  and  ratios  of  about  1000  to  1  were  used. 
Presently  available  techniques  allow  for  significant  im- 
provements in  both  of  these  quantities,  which  would  hope- 
fully result  in  more  accurate  points  for  curves  of  p'  and 

R 

p"   from  which  P^    could  be  calculated. 

Coefficients  of  the  Power  Series  Expansion 
Through  Zero  Frequency  Extrapolation 

Though  not  as  sophisticated  as  the  methods  previously 
presented,  the  P^^)^  P^^\    and  P^^^  coefficients  can  be 
calculated,  respectively,  from  the  zero  frequency  inter- 
cept of  the  curves  a  versus  w  (intercept  [P^^]-'-'^), 
a2-^2  (intercept  P^°^),  2aC/a)  (intercept  P^^^),  (a^-c^-pO ) /u)2 
(intercept  P^^))^  ^^^    [;(g2  +  ^2)2  .  (p(0^)2]/(^2  (intercept 

(P    )2  +  2P^'^^p^^^).   These  curves  were  determined  and 
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extrapolated  to  zero  in  an  attempt  to  extract  the  coefficients 
Use  was  first  made  of  the  a,  a^-c^,  and  2aC  curves  to  de- 
termine P^°^  and  p(l)    Tv,o««    -1 

.   These  values  were  then  used  in  the 

determination  of  that  value  of  p(2)  which  yielded  the  best 
intercepts  for  the  latter  two  curves.   p(3)  ,^^    calculated 
from  the  slope  of  2aC/u,.   Shown  below  are  the  pCn)  ^3  ^^_ 
termined  by  this  method. 


p(0) 


(63.87  i    .09)  X  10-^  cm-2 


(4.64  t    .05)  X  10-6  ge^  ^^.j 


p(2) 


P 


(3) 


*(3.8)  i  .3)   X  10-10  3ec2c„,.2 
-(1.0   ±  .3)   X  10-14  3ec3cm-2 


Summary  of  the  Kxperimental  Results 
And  uiscussion  ()[:  the  ^.rror  AnalTiT-, 

Shown  in  Table  U  are  theoretical  values  for  the 
first  four  p(">  coefficients,  along  with  the  experimental 
values  obtained  by  the  various  methods  of  data  analysis 
presented  in  this  chapter.   It  is  clear  from  the  table 
that  consistent  experimental  values  for  the  n(n)  ^^^^ 
obtained  in  every  case  and  that  theory  and  experiment  began 
to  diverge  after  pd^. 

The  determination  of  the  errors  associated  with  the 
a  and  ,    values  was  based  upon  the  reproducibility  of  points 
on  the  amplitude  and  phase  curves   and  the  variation  of  a 
and  C  as  the  low  z  points  were  successively  eliminated. 
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These  errors  were  found  to  be  on  the  order  of  1  per  cent 
for  the  lower  frequencies  and  about  5  per  cent  at  the 
higher  frequencies.   A  standard  propagation  of  errors 
technique  resulted  in  estimates  of  the  errors  associated 
with  the  experimental  curves  derived  from  a  and  C.   The 
same  procedure  was  used  in  the  curves  which  were  based 
upon  the  derivatives  of  a  and  C.   Errors  in  the  p^'^^  as 
listed  are  consistent  with  the  ones  calculated  by  the 
STEP-WISE  MULTIPLE  REGRESSION  code. 


CHAPTER  VI 


RELATIONSHIPS  BETWEEN  THE  NEUTRON  WAVE  EXPERIMENT 
AND  DIE-AWAY  EXPERIMENT  IN  THE  DETERMINATION 
OF  DIFFUSION  AND  THERMALIZATION  PARAMETERS 


The  main  objective  of  this  dissertation  has  been 
to  demonstrate  the  effectiveness  of  the  neutron  wave 
technique  for  the  determination  of  diffusion  and  therma- 
lization  parameters.   In  order  to  evaluate  the  fulfillment 
of  this  objective,  the  results  obtained  here  for  the 
experiment  in  graphite  will  be  compared  to  published 
results  of  die-away  experiments  where  the  same  parameters 
were  measured.   The  various  complementary  relationships 
between  the  two  experiments  will  also  be  investigated. 

In  the  die-away  experiment  one  measures  the  decay 
constant  of  the  fundamental  time  and  energy  dependent 
eigenfunction   of  the  Boltzmann  operator  associated  with 
the  fundamental  spatial  mode.   The  real  eigenvalue,  X, 
is  determined  as  a  function  of  the  system's  lowest  geo- 
metric buckling;  and  this  A  versus  B^  curve  yields  the 
parameters  of  interest.   Consistent  with  current  termi- 
nology, this  method  will  be  referred  to  as  a  dynamic 
experiment.   The  neutron  wave  experiment,  in  contrast,  is 
a  static  experiment  (i.e.,  in  the  same  class  as  the 


124 


125 

poisoning  experiment);  and  one  measures  the  complex, 
fundamental  eigenvalue,  p,  of  the  space  and  energy 
dependent  Boltzmann  operator  as  a  function  of  frequency. 

The  complementary  nature  of  the  two  experiments 
is  expressed  pictorially  in  Figure  30.    Shown  in  the 
first  quadrant  of  the  figure  are  results  of  a  typical 
die-away  experiment,  and  in  the  third  quadrant  the 
modulus  of  p^  as  determined  in  this  work  is  plotted  as 
a  function  of  frequency  up  to  1000  cps.   It  is  interesting 
to  note  that  B^,  the  independent  variable  of  the  die- 
away  experiment,  and  |p^|,  the  dependent  variable  of 
the  neutron  wave  experiment,  can  be  mathematically  equated. 
This  result  is  not  surprising  when  one  recalls  that  p  is 
the  complex  eigenvalue  associated  with  the  space  decay  of 
the  neutron  wave  and  B^  is  the  eigenvalue  of  the  funda- 
mental spatial  dependence  of  the  die-away  experiment.   One 
can  also  devise  corr^alations  between  o)  and  time  after  the 
pulse  in  the  die-away  experiment.   This  is  most  easily 
explained  if  one  considers  the  expansion  of  the  energy 
dependence  of  the  neutron  wave  experiment  in  the  time- 
energy  eigenvalues  of  the  Boltzmann  operator  (i.e.,  the 


*   The  die-away  points  were  obtained  from  Reference  (29). 
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.  ft 

energy  modes  of  the  die-away  experiment).    If  this 

expansion  is  used,  very  similar  results  to  those  pre- 
viously presented  are  obtained.   As  frequency  increases, 
the  higher  energy  modes  begin  to  contribute  proportionally 
more  to  the  total  energy  dependence  which  is  dominated  at 
low  f  by  the  fundamental  energy  mode.   This  is  exactly 
the  situation  that  exists  at  short  times  after  the  pulse 
in  the  die-away  experiment.   Thus  in  terms  of  energy  mode 
contamination,  frequency  in  the  neutron  wave  experiment 
and  time  after  the  pulse  in  the  die-away  experiment 
are  inversely  related. 

One  can  also  represent  the  poisoning  experiment  in 
Figure  30.   A  curve  very  similar  to  the  one  associated 
with  the  wave  experiment  would  be  obtained  if  the 
poisoning  experiment's  squared  inverse  relaxation  length 
is  plotted  on  the  same  axis  as  |p^|  and  ^  v^  ,  the  added 
poison  times  v,  plotted  with  w.   It  is  for  this  reason 
that  the  jw/v  term  of  the  neutron  wave  experiment  is 
sometimes  equated  to  a  complex  1/v  poison.   Note  also 
that  the  point  on  the  figure  which  represents  the  con- 
ventional exponential  experiment  is  common  to  both  the 

*   This  was  actually  done  by  the  author  for  the  MASSMK 

kernel  using  45  energy  points.   The  fundamental  complex 
inverse  relaxation  length  thereby  obtained  was  in 
excellent  agreement  with  the  one  presented  in  this 
work.   This  result  indicates  that  33  energy  points 
are  adequate  for  theoretical  results  to  be  compared 
with  the  experiment  and  that  the  energy  expansion 
as  presented  was  valid. 
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poisoning  and  neutron  wave  experiment  and  that  it  is 
represented  as  the  a)  =  0  or  ^^v  =  0    intercept.   Herein  lies 
the  justification  for  the  reference  of  the  space-energy 
eigenfunct ions  of  the  Boltzmann  operator,  used  in  the 
energy  expansion  of  this  dissertation,  as  eigenfunctions 
of  the  exponential  exneriment. 

Perez  and  Uhrig  (9)  have  shown  that  the  eigenvalues, 
^   and  p^,  are  mathematically  related  when  p^  is  expressed 
in  the  form 

n=l  ^o  (78) 

and  use  is  made  of  the  relationships 


jw   =   -X"  (79) 

where  the  w  superscript  indicates  that  X"  is  the  eigen- 
value that  is  to  be  obtained  from  the  0^"^^  in  a  neutron 
wave  experiment  and 

-p   =   jBg  (80) 

Substitution  of  the  above  transformation  into  equation 
(78)  yields 
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N 


.B2   =  -(B2  +  b2)     I  Q^"^  (  1°_I )''        (81) 

Z  J-  \7 


n  =  l  ^o 


°        ^■^-^^r,Q■,•r,n2 


Bv  inverting  the  series  in  ( )  to  one  in  B"^  , 

V. 


o 


one  obtains  for  N  =  3 

XW   =   cx^  +  DW  b2  -  CWB'*  +  F^  B^  (82) 


where 


Vo^^ao  (83) 


D^   =   V  /Q 
o       o 


(1) 


Co   =  -v,Q^2V(Q^l))3 

F^   --   V  /(q(1))S[2(0(2))2  _  qCDqO)-, 


(84) 


(85) 


(86) 


The  coefficients  as  defined  above  differ  slightly  from 
the  ones  measured  in  the  dynamic  die-away  experiment 

because  they  were  derived  from  the  static,  wave  experiment 

ft 

expressions. 


*      In  this  context  the  difference  between  static  and  dynamic 
experiments  is  related  to  the  effect  upon  the  energy 
dependence  of  a  1/v  absorber  in  the  medium.   In  the 
die-away  experiment  a  1/v  absorber  does  not  affect  the 
energy  dependence  of  the  fundamental  eigenfunction  (19) 
which,  for  zero  buckling,  is  a  Maxwellian.   The  funda- 
mental eigenvalue  of  the  exponential  experiment  is  a 
Maxwellian  distribution  only  in  the  zero  absorption 
case.   Thus  the  static  and  dynamic  experiments  are 
equivalent  only  when  aQ=0  in  the  expressions  for  the 

static  experiment. 
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In  the  analysis  of  data  presented  in  the  previous 
chapter,  o^    was  expressed  in  the  form 

N-1  N 

n=0  n=l 

(EVEN)  (ODD)  (87) 

From  equations  (78)  and  (87)  the  following  connection 
between  the  0^"^^  and  P^"^  coefficients  is  obtained  with 

N=3 

N 


n  =  l 


^o  (88) 


P^l^  =  -  {Q^^^2(^)Q^^^3(^)2q(3)j 


(89) 


Vo  °  (90) 


P^3)  _    J_  ^(3) 

^^3  (91) 


The  Q^"^  and  —  obtained  from  the  P^"^  by  the  iterative 

^o 

use  of  these  expresions  can  be  used  to  determine  D  ,  C"  , 

^  o »   o 

and  F^.   For  the  determination  of  these  same  parameters 
o 

for  comparison  to  die-away  experiments  by  use  of  ex- 
pressions (84),  (85),  and  (86),  the  above  expressions 
with  0^=0  should  be  used  to  determine  the  Q    . 
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Table  5(a)  shows  the  values  of  a^ ,  Q^'^\    Q^^\ 
and  Q^"^^  obtained  from  the  experimental  data  as  well  as 
theoretical  results  obtained  by  use  of  the  SUMMIT  kernel. 
In  Table  5(b)  these  same  oarameters  derived  by 
the  use  of  equations  (89),  (90),  and  (91)  with  aQ=0  are 
presented.   From  these  two  tables  and  equations  (83),  (8U), 
(85),  and  (86),  the  values  of  both  the  static  and  dynamic 
Qq,  Dq ,  C^,  and  F^  parameters  were  calculated  and  are 
shown  in  Table  6.   In  Table  7  the  values  for  the  various 
diffusion  and  thermalization  parameters  that  are  to  be 
compared  with  the  die-away  results  are  listed  along  with 
the  results  from  several  die-away  experiments.    From  an 
inspection  of  these  tables  one  can  conclude  the  following: 

(1)  The  static  and  dynamic  parameters  differ 
by  less  than  5  per  cent  with  most  differences  on  the 
order  of  1  per  cent. 

(2)  The  experimentally  obtained  diffusion 


*   Dq  and  Cq  were  also  calculated  using  the  SUMMIT  kernel 
by  Ohanian's  method  with  both  33  points  and  4  5  points 
in  the  discrete  energy  representation.   The  values  ob- 
tained were:   Dq  =  (2.092)  x  10^)  cm  sec"^,   Cq  = 
(20.9014  X  10^)  cm2sec-l,  33  points;   D^  =  (2.1435  x  10^) 
cm  sec-l,  Co  =  (21.128  x  10^)  cm^sec-i,  ^^  points. 
These  values  agree  very  v;ell  with  the  ones  in  Table  6 
especially  when  one  recalls  that  the  time-energy  eigen- 
functions  were  used  in  the  Ohanian  method  and  the 
space-energy  eigenfunct ions  were  used  in  this  work. 


132 


TABLE  5 

THE  EXPERIMENTAL  AND  THEORETICAL  VALUES 
FOR  THE  Q(n)  COEFFICIENTS 


Oq/Vo 


s(l) 


EXPERIMENT 

(p  =  1.60  gm/cm^) 

THEORY 

(o  =  1.60  gm/cm-^) 


U.14 
i  .03 


U  .  0  9  6  5 


1.036 
i  .005 


1.0505 


,(2) 


,(3),,_2 


Q   xlQ- 


-19.1 

i  .9 


-10.9455 


1.09 
+  .05 


4.7582 


(a)   a„  5<  0 


EXPERIMENT 

(p  =  1.60  gm/cm^) 


THEORY 


1.020 
±  .005 


1.0415 


-19.0 
i   .9 


-10.3595 


1.09 
+  .05 


4.7682 


(b)  a^    Set  Equal  to  Zero  to 
Calculate  Q^"^ 
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TABLE    6 


THE    STATIC   AND   DYIJAMIC    PARAMETERS 


a      sec~l 


DqX10-5 


C    xlO-5  F„xl0-7 

o  o 


CTr.2sec~l  cra'^sec"!         cm^sec-l 


EXPERIMENT 
(p  =  1.60 
gm/cm^ ) 


91.  t    1       2.12  +  .01     38.  t    2  12.    t    2 


Made  zero  to 

obtain  Q(n)    2.16  ±  .01     39.  ±2  12.  ±  2 


90.04 


2.0838 


20.9110     -4.4891 


THEORY 
(p  =  1.60 
gm/cm^ ) 


Made  zero  to 

obtain  Q^n)    2.1104 


20.1S53     -5.0574 


13H 


TABLE  7 


PARAMETERS  AS  DETERMINED  FROM  VARIOUS 
DIE-AWAV  EXPERIMENTS  COMPARED  WITH 
THE  RESULTS  OF  THIS  WORK 


Year        Authors      ao(sec-i)  ^o^^^-^        CoxlO-5   FqxIO-^ 

cm2sec-l   cm'*sec-l  cm'^sec-1 


1961  Lalande          79.5*3'"  2. 07*. 04  25.3±2  

Starr  S  Price 

1962  type  GBF         7U.5±.06  2.14±.01  34±3  

type  AA          68.1*. 08  2. 19*. 03  41*4  

38±5  

88.3*1.2  2. 13*. 02  25*5  

88.6*1.6  2. 11*. 02  16*5  -20*10 

1963  Sagot  et  a]_.  (30)  7U.1*3  2.19*. 03  37.9*H  

1955      This  work        91.  *1  2. 16*. 01  39.  *2  12. *2 


1962      S^^^r  ^ 

de  Villiers 


1962     Klose  et  al . ' 


"   All  values  listed  in  the  table  have  been  corrected  to 
a  density  of  1.50  gm/cm^ . 

t   Both  the  three  and  four  parameter  fits  are  included. 
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coefficient  agrees  quite  well  with  those  of  pulsed 
experiments  and  with  the  theory. 

(3)  The  value  obtained  in  this  work  for  the 
cooling  coefficient,  C^ ,  is  about  a  factor  of  two 
higher  than  the  one  predicted  by  use  of  the  SUMMIT 
kernel  but  is  in  very  good  agreement  with  the  more 
recent  die-away  results.   However,  this  result  is 
still  somewhat  lower  than  the  lower  bound  of  45  x 
10  5  cm'*sec-i  predicted  for  Cq  by  de  Saussure  (31). 

(4)  The  value  of  F^  obtained  experimentally 
differs  both  in  sign  and  magnitude  from  the  one 
predicted  here  by  the  use  of  the  SUMMIT  kernel. 

It  also  disagrees  both  in  sign  and  magnitude  with 
the  value  obtained  by  Klose,  et  a}^.    (3_2).   Honeck 
(  33)  used  the  SUMMIT  kernel  and  obtained  a  value  of 
of  -7.4U  X  10^  cm^sec-^  for  Fq.   However,  Melkin 

( 3U )  predicts  a  positive  coefficient  of  the  order 

2 

of  Cq/Dq  which  is  consistent  with  equation  (86)  if 

one  neglects  the  contribution  of  0^^'.      An  inspection 
of  equation  (36)  indicates  that  the  sign  of  F^  is 
determined  by  the  magnitude  of  2(Q^2^)^  compared  to 
Q   Q    .   The  solution  to  this  contradiction  awaits 
more  exacting  theoretical  and  experimental  work. 


CHAPTER  VII 
CONCLUSIONS  AND  RECOMMENDATIONS 

The  space,  energy,  and  time  or  frequency  dependent 
Boltzmann  equation  in  the  diffusion  approximation  has  been 
used  within  the  framework  of  neutron  thermalizat ion  theory 
to  predict  thermal  neutron  wave  propagation  in  a  moderating 
medium.  This  theory  was  then  compared  with  the  experiment 
conducted  on  a  graphite  pile. 

In  the  theoretical  section  the  energy  dependence 
was  expanded  in  the  space-energy  eigenfunctions  of  the 
Boltzmann  operator,  and  the  transverse-space  dependence 
was  expanded  in  the  orthogonal  solutions  to  the  Helm- 
holtz  equation  with  homogeneous  boundary  conditions. 
It  was  then  shown  that  each  of  these  space  and  energy 
modes  thus  formed  propagates  with  characteristic   complex 
inverse  relaxation  lengths.   The  perturbation  method 
developed  by  Perez  was  used  to  expand  these  complex  eigen- 
values in  a  Maclaurin's  series  of  the  frequency  variable. 
In  cases  where  this  series  did  not  rapidly  converge,  a 
direct  solution  of  the  complex  secular  determinant  was 
obtained.   Two  scattering  kernels  were  selected  for 
implementation  into  the  theory,  the  MASSMK  kernel  and  the 
SUMMIT  kernel,  which  is  based  primarily  upon  derivations  by 
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Parks  (2_2_).   The  numerical  calculations  of  the  theory 
were  performed  on  an  IBM  709  digital  computer  and 
founded  upon  the  previous  work  of  Ohanian  (19). 

In  the  experimental  section  an  accurate  technique 
for  conducting  neutron  wave  experiments  in  multiplying 
or  diffusing  media  was  presented  along  with  the  results 
of  the  graphite  experiment.   The  excitation  of  neutron 
waves,  both  by  modulated  and  pulsed  sources,  and  the 
electronic  systems  associated  with  their  detection  were 
discussed.   Special  attention  was  given  to  the  analysis 
of  data  for  which  several  computer  codes  were  developed. 

In  the  comparison  of  theory  and  experiment  it 
was  shown  that  the  neutron  wave  technique  allows  for  the 
separation  of  diffusion  and  thermalization  effects,  there- 
by yielding  a  more  exacting  comparison  of  theoretical  and 
experimental  results.   Polynomial  fits  of  functions 
formed  from  the  experimental  results  yielded  coefficients 
to  be  compared  with  the  theoretically  derived  terms  of 
the  perturbation  expansion.   The  results  of  the  thermal 
pulse  propagation  experiment  were  also  analyzed  by  a 
method  suggested  by  Moore,  and  confirmation  of  his  pre- 
dictions was  established. 

The  diffusion  and  thermalization  parameters  determined 
from  the  data  analysis  were  then  compared  with  results  of 
several  die-away  experiments  and  with  the  theory.   The 
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complementary  relationships  between  the  neutron  wave 
and  die-away  experiments  were  also  presented. 

The  more  significant  results  of  this  dissertation 
are  listed  below: 

A.  It  was  demonstrated  that  the  neutron  wave 
technique  is  a  very  powerful  method  for  the  determination 
of  thermalization  and  diffusion  parameters. 

B.  The  numerical  methods  developed  by  Ohanian 
for  the  determination  of  the  time-energy  eigenfunctions 
of  the  Boltzmann  operator  and  investigation  of  the  die- 
awav  experiment  were  successfully  applied  to  the  deter- 
mination of  the  space-energy  eigenfunctions  of  the 
Boltzmann  operator.   It  was  shown  that  these  latter 
eigenfunctions  are  especially  well  suited  for  theoretical 
studies  of  the  neutron  wave  experiment. 

C.  Pligher  energy  mode  contamination  of  the 
fundamental  complex  inverse  relaxation  length  is  predicted 
by  the  SUMMIT  scattering  kernel  for  frequencies  above  300 
cps  when  the  transverse  buckling  is  on  the  order  of 

6  X  lO""*  cm~2.   The  frequency  at  which  this  contamination 
becomes  significant  is  dependent  upon  the  transverse 
buckling,  and  its  severity  is  not  predicted  by  less 
accurate  scattering  kernels  than  the  one  incorporated 
into  the  SUMMIT  code. 

D.  The  experimental  equivalence  of  excitation 
of  neutron  waves  by  either  modulated  or  pulsed  sources 
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has  been  confirmed.  It  was  found  that,  in  general,  the 
thermal  pulse  propagation  technique  is  more  useful;  but 
it  is  recommended  that  the  results  thereby  obtained  be 
always  "spot  checked"  by  the  use  of  a  modulated  source. 
The  advantages  inherent  in  these  two  methods  are  listed 
below : 

(1)  Advantages  of  the  modulated  source. 

(a)  It  is  much  simpler  to  check  the 
accuracy  of  the  experimental  results 
and  the  accuracy  of  the  method  of 
data  analysis. 

(b)  Hand  calculations  and  integral  count 
methods  can  be  used  to  analyze  the 
data. 

(c)  Resonance  phenomena  can  be  more  easily 
detected.   One  could,  for  example,  set 
the  input  to  the  resonance  point  through 
hand  calculations  of  the  output  at 
several  frequencies  or,  perhaps,  by 
analog  circuits  connected  to  a  con- 
tinuous representation  of  the  sinusoidally 
varying  count  rate. 

(d)  A  sinusoidal  source  is  very  useful  when 
crosscorrelat ion  techniques  are  to  be 
applied.   In  some  measurements  one  cannot, 
or  does  not  want  to,  perturb  the  flux  to 
any  significant  degree.   For  this  situation 
a  sinusoidal  input  together  with  cross- 
correlation  methods  of  data  analysis  may 
prove  to  be  useful. 

(e)  The  data  analysis  is  not  comnlicated  by 
the  approximation  of  a  Fourier  transform 
with  a  finite  numerical  integration. 

(2)  Advantages  of  the  thermal  pulsed  source. 

(a)   One  is  able  to  extract  many  frequency 

points  on  the  system's  transfer  function 
from  a  single  measurement  in  the  system. 
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(b)  The  system's  response  to  thermal 
pulses  can  be  used  as  suggested  by 
Moore  to  determine  derivatives  of 
the  complex  inverse  relaxation  length 
directly  from  the  experimental  results. 

(c)  The  epicadmium  background  and  room 
return  are  present  only  at  the  beginning 
of  the  thermal  pulse.   Therefore,  less 
elaborate  methods  are  needed  for  the  sub- 
traction of  these  unwanted  effects. 

(d)  The  D.C.  background  can  be  made  as 
small  as  the  neutron  generator  will 
allow.   Ratios  of  1000  to  1  have  been 
obtained  here  while  with  a  modulated 
source  the  maximum  ratio  is  1  to  1. 

For  this  reason  one  can  generally  detect 
more  decades  of  decay  of  frequencies 
excited  by  a  thermal  pulse  source  than 
those  from  a  modulated  source. 

(e)  Reflection  of  thermal  pulses  at  an 
interface  may  prove  to  be  a  very 
useful  way  for  testing  various  inter- 
face boundary  conditions. 


E.  The  predictions  by  Moore  in  his  suggested  method 
of  analysis  of  thermal  pulse  propagation  data  have  been 
experimentally  verified. 

F.  The  relationships  between  the  die-away  and 
neutron  wave  experiments  and,  in  general,  between  dynamic 
and  static  experiments  have  been  more  clearly  defined 
through  the  theoretical  work  presented  here  combined  with 
the  previous  work  of  Ohanian  (19). 

G.  In  connection  with  the  comparison  of  experi- 
mental work  previously  done  in  neutron  wave  propagation  in 
moderating  media,  the  following  extensions  have  been  made. 
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1.  A  neutron  generator  has  been  used  with  an 
electronically  varied  source  intensity  rather 
than  mechanical  sources.   This  increases  the 
frequency  range  where  data  can  be  obtained 
and  allov;s  for  simpler  expressions  for  the 
spatial  dependence  of  the  source. 

2.  The  data  analysis  is  based  upon  Fourier 
series  expansions  of  Fourier  transforms 
rather  than  upon  the  area  method.   The  Fourier 
method  used  here  allows  for  more  accurate 
source  normalizations  (independent  of  frequency 
harmonic  distortion)  and  more  exacting  extrac- 
tion of  the  amplitude  and  phase  values. 

3.  The  experiment,  because  of  the  above  extensions, 
has  been  applied  to  the  thermalizat ion  problem. 
Previous  experimentalists  concentrated  upon 
diffusion  effects. 

It  is  evident  from  the  results  presented  in  this 
dissertation  and  from  those  summarized  in  the  Review  of 
the  Literature  that  much  more  can  still  be  done  in  the 
area  of  neutron  wave  experimentation.   A  logical  extension 
of  the  theory  presented  here  is  to  make  the  more  exacting 
comparison  between  theory  and  experiment  by  generating 
a  "theoretical  experiment."   That  is,  to  synthesize  the 
frequency,  space,  and  energy  dependent  flux  predicted  by 
the  theory  for  several  w  and  z  values  and  compare  this 
result  to  the  experimental  measurements.   It  was  assumed 
here  that  the  behavior  of  the  complex  eigenvalue  associated 
with  the  measured  flux  was  well  represented  by  the  fundamen- 
tal eigenvalue.   This  assumption  was  shown  to  be  valid  for 
low  frequencies  for  the  experiment  presented  but  of 
questionable  validity  at  high  frequencies. 

The  experiment  could  be  improved  in  several  ways. 
First,  it  is  expected  that  less  energy  mode  contamination 
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would  be  present  if  a  larger  system  were  investigated.   In 
the  experiment  presented,  the  transverse  leakage  was  much 
larger  than  the  zero  intercepts  of  the  complex  inverse 
relaxation  lengths  associated  with  the  various  energy 
modes.   If  the  opposite  were  true,  differences  in  these 
intercepts  would  cause  a  proportionally  wider  soread  in 
the  various  p's  at  low  frequency,  thereby  causing  cross- 
over points  of  the  a  curves  versus  frequency  to  occur  at 
higher  frequencies.   Improvements  can  also  be  made  in 
the  error  analysis  which,  for  this  work,  could  have  been 
biased  by  the  lack  of  consistent  statistical  techniques. 
The  fact  that  little  variation  was  observed  in  the  p's 
as  determined  by  the  variance-weighted,  STEP-WISE  MULTIPLE 
REGRESSION   code  could  have  been  fortuitous.   More  con- 
fidence could  have  been  placed  in  this  result  if  the 
weights  had  been  more  rigorously  determined. 

Significant  improvements  can  be  made  in  the  attain- 
ment of  pulse  propagation  data  to  be  analyzed  by  the  method 
suggested  by  Moore.   The  refinements  which  appear  to  be 
most  profitable  are  decreasing  the  channel  width  of  the 
multi-channel  analyzer,  increasing  the  peak  to  background 
ratio   of  the  pulse,  and  decreasing  the  pulse  width. 

Much  can  also  be  done  to  improve  the  source  used 
for  the  experiment.   A  thermal  source  of  a  single  energy, 
as  can  be  obtained  by  use  of  a  Van  de  Graaff  accelerator, 
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is  surely  more  desirable  than  the  neutron  generator  and 
thermalizing  tank  used  for  this  work.   Also,  the  elimi- 
nation of  higher  spatial  modes  by  use  of  a  source  which 
contains  only  the  fundamental  spatial  mode  is  of  un- 
questionable value. 

The  extension  of  the  neutron  wave  technique  to 
other  areas  of  application  has  already  begun.   Some  of 
these  areas  are  listed  for  completeness: 

(1)  The  data  analysis  and  experimental  technique 
discussed  here  could  be  applied  to  fast  waves  and  fast 
systems.   In  fact,  a  good  starting  point  for  this  inves- 
tigation could  be  the  analysis  of  the  epicadmium  pulse 
propagation  data  as  detected  in  the  experiment  done  by 
Hartley. 

(2)  Some  work  has  been  done  with  respect  to  the 
propagation  of  neutron  waves  in  heterogeneous  assemblies. 
It  is  hoped  that  neutron  waves  will  be  sensitive  to  the 
lattice  parameters  associated  with  these  assemblies. 

(3)  The  most  obvious  extension  of  this  work  has 
been  applied  by  several  experimentalists.   It  is  the 
investigation  of  thermal  multiplying  systems  by  the 
neutron  wave  technique.   Still,  much  can  and  surely  will 
be  done  in  this  area  in  the  future. 
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APPENDIX 


RESULTS  OF  THE  DATA  ANALYSIS 

This  Appendix  is  divided  into  three  parts.   In 
Part  1  the  amplitude  and  phase  of  the  neutron  wave 
versus  distance  from  the  source  is  listed  for  all  of 
the  frequencies  analyzed.   The  error  associated  with 
these  data  has  been  estimated  to  be  about  l  per  cent. 
Generally,  results  of  the  wave  experiments  are  listed 
first  and  followed  by  points  obtained  from  the  pulse 
propagation  experiment.   These  data  were  used  on  the 
final  channel  chop  sequence  from  which  the  real  and 
imaginary  components  of  the  complex  inverse  relaxation 
length  were  determined.   The  data  have  been  normalized 
to  the  first  z  position  at  which  measurements  were  made 
which  was  always  within  five  and  usually  about  two 
centimeters  from  the  source  boundary. 

In  Part  2  the  real  and  imaginary  components  of 
the  complex  inverse  relaxation  length  (slopes  of  the 
curves  of  Part  1)  are  shown  along  with  their  error 
estimates  for  each  of  the  frequencies  analyzed.   The 
squared  complex  inverse  relaxation  length  is  shown  in 
Part  3  for  each  frequency  analyzed.   Error  estimates  are 
also  shown. 
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PART    1       AMHLITUDb    AND    PHASt    CURVCS 


RESULTS    OH  10.    CPS.    DATA    ANALYSIS 

POSITION                             AMPLUUDt  PHASt 

0.  O.IOOOt    01  0. 

O.bOlOt    01  0.6208E    00  0.7600E-02 

0.989Ut    01  0.^114E-00  0.1550t-0l 

O.UayE    02  0.2b39e-00  0.2^40E-0l 

0.19/8E    02  0.1770E-00  0.3240E-01 

0.2478t    02  0.1092E-00  0.'t240E-0l 

0.296dE    02  0.7470E-01  O.i>17OE-0l 

0.3467E    02  0.4901E-01  0.5970C-0i 

0.3a57E    02  0.3845E-01  0.65a0E-0l 

0.4443E    02  0.2366G-01  0.7670F-01 

0.4946E    02  0.1588E-01  0.8510E-01 

0.543tE    02  0.1038E-01  0.9460t-0l 

0.5934C    02  0.7392b-02  0.1022E-00 


RESULTS    OF  15.    CPS.    DATA    ANALYSIS 

POSITION  AMPLITUDE  PHASE 

0.  O.IOOOE    01  0. 

0.5010E    01  0.6208E    00  0.1141C-01 

0.9890E    01  0.4113E-00  0.2336E-QI 

0.L489E    02  0.2539E-00  0.3672E-01 

0.I978E    02  O.lirrOE-00  0.4a70E-01 

0.247aE    02  0.1092E-00  0.6371E-01 

0.2968E    02  0.7464E-01  0.7758E-01 

0.346i'E    02  0.4898E-01  0.8936E-01 

0.385/t    02  0.3a42E-0l  0.9883E-01 

0.4445C    02  0.2364E-01  0.1152E-00 

0.4946t    02  0.1586t-0l  0.1277E-00 

0.5434E    02  0.10B7E-01  0.l420t-00 

0.b934e    02  0.7383E-02  0.1533E-00 


RESULTS    OF  25.    CPS.    DATA    ANALYSIS 

POSITION  AMPLITUDE  PHASE 

0.  O.IOOOE    01  0. 

0.5010E    01  0.6156E    00  0.1909E-01 

0.9890E    01  0.407aE-00  0.3909E-01 

0.1489E    02  0.2bl6E-00  U.6139E-01 

0.1978C    02  0-1733E-00  0.8139E-01 

0.2473E    02  0.1081E-00  0.1064E-00 

0.296aE    02  0.7386E-01  0.1295E-00 

0.346/E    02  0.4843E-01  0.1495E-00 
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0.3857E  02 
0.444DE  02 
O.-^ry^bf  02 
0.bA34E  02 
0.5934E  02 


0.3800E-01 
0.2337t-0l 
0.l568fc-01 
0.1074E-01 
0.729Ifc-02 


0.1651E-00 
0.1924E-00 
0.2134^-00 
0.2370E-00 
0.2560t-00 


RESULTS  OF 

50.  CPS.  DATA 

ANALYSIS 

POSITION 

AMPLITUDE 

PHASE 

0. 

O.IOOOE  01 

0. 

O.^iOlOE 

01 

0.6199E  00 

O.3830t-0l 

0.9890E 

01 

0.4099t-00 

O.7B40t-0l 

0.i489L 

02 

0.2'S24t-00 

0.1232E-00 

0. 19780 

02 

0-1737E-00 

0.1637F-00 

0.?'»7afc 

02 

o.ioaoE-00 

0.21346-00 

0.296(it 

02 

O.73'i6E-0l 

0.2'>94e-00 

0.3467E 

02 

0.4820t-0l 

0.3001E-00 

0.385^E 

02 

n.3780E-0I 

0.3318E-00 

0.444'iE 

02 

0.2316E-01 

U.3862e-00 

0.49A6E 

02 

0.1556E-01 

0.4291E-00 

0.5434E 

02 

0.1061E-01 

u.476ie-oo 

0.3934b 

02 

0.7186t-02 

0.5i47E  00 

RESULTS  OF 

63.  CPS.  DATA 

ANALYSIS 

PUSITION 

AMPLITUDfc 

PHASF 

0. 

O.lOOOfc  01 

0. 

0.5010G 

01 

0.6l90t  00 

0.5047E-01 

0.9890^ 

01 

0.4086t-00 

0.1028E-00 

0. 1489b 

02 

0.2510E-00 

0.1611E-00 

o.i97ac 

02 

0.l744b-00 

0.2140E-00 

0.247ab 

02 

0.1069E-00 

0.2783E-00 

0.2968t 

02 

0.7264E-01 

0.3376E-00 

0.34671: 

02 

0.47b2b-0l 

0.3912E-00 

0.38!>7E 

02 

U.3724b-0l 

0.433lt-00 

0.4445b 

02 

U.2273E-01 

0.5035E  00 

0.4946E 

02 

0.i528E-0l 

0.5603E  00 

0.5434E 

02 

0.l039b-01 

0.6209E  00 

0.5934e 

02 

0.7017E-02 

0.6716E  00 

RESULTS  OF 

75.  CPS.  OAIA 

ANALYSIS 

PosiriuN 

AMPLITUDE 

PHASfc 

0. 

O.IOOOC  01 

0. 

O.bOlUb  01 

O.Aiaib  00 

0.5a40E-0l 

0.9890b  01 

0.4073b-00 

0.1187E-00 

0.1489b  02 

0.2498E-00 

0.1859E-00 

0.l97ab  02 

0.1733E-00 

U.2480t-00 

0.24/8C  02 

0.1060b-U0 

0.3209E-00 
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0.2968C 

02 

0.7185E-01 

0.3890E- 

-00 

0.3A6  7E 

02 

0.4694E-01 

0.4312E- 

-00 

0.3857t 

02 

0.3676t-0i 

0.5002E 

00 

0.^^4t>b 

02 

0.2240E-01 

o.saiot 

00 

0.A946E 

02 

0.1504E-01 

0.6473E 

00 

0.t)43<»t 

02 

0.l019t-0l 

0.7179E 

00 

0.b934E 

02 

0.6873E-02 

0.7738E 

00 

RESULTS  OF 

100.  CPS.  DATA 

ANALYSIS 

POSITION 

AMPLITUDE 

PHASE 

0. 

O.IOOOE  01 

0. 

O.bOlOc 

01 

0.6150E  00 

0.7745fc-0l 

0.9890E 

01 

0.403lt-00 

0.1572E-00 

0.1489t 

02 

0.2457t-00 

0.2460E-00 

0.1973E 

02 

0.1696E-00 

0.32  77E-00 

0.2478E 

02 

0.1030E-00 

0.4236E-00 

0.2968t 

02 

0.6933E-01 

0.5124E  00 

0.346  7t 

02 

0.4b04E-01 

0.5964t  00 

0.3a5/t 

02 

0.3bl6E-01 

0.6633t  00 

0.44A3E 

02 

0.2124E-01 

0.7690t  00 

0.4946E 

02 

0.1425E-01 

0.8602E  00 

0.5434b 

02 

0.956lfc-02 

0.9506E  00 

0.b934E 

02 

G-64UIE-02 

0.1030E  01 

RESULTS  OP 

110.  CPS.  UAIA 

ANALYSIS 

POSinON 

AMPLITUDE 

PHASt 

0. 

O.IOOOE  01 

0. 

0.4000t 

01 

0.7123K  00 

0.8506E-01 

o.889at: 

01 

0.4779E-00 

0.1788E-00 

0.i289L 

02 

0.3337b-00 

0.2335E-00 

0.l778t 

02 

0.2194E-00 

0.3622b-00 

0.2667fc 

02 

O.lOOlF-00 

0.5001E  00 

0.2173E 

02 

0.l446b-00 

0.4134E-00 

0.306^c 

02 

0.6564E-01 

0.!>902t  00 

0.1089E 

02 

-0.2848E  01 

-0.1794E-00 

0.5880E 

01 

-C.1747t  01 

-0.9476E-01 

O.iOOOC 

01 

-0.1142E  01 

~0.7700E-02 

0.4000E 

01 

0.6942E  00 

0.8920E-01 

0.8890t 

01 

0.4779E-00 

0.l7a8E-00 

0.1389c 

02 

0.2893E-00 

0.2830E-00 

0.1879E 

02 

0.1941b-00 

0.3795F-00 

0.23/8t 

02 

0.l2l)7fc-00 

0.4719E-00 

0.2768r 

02 

0.9796E-01 

0.5460G  00 

0.3356E 

02 

0.5897E-01 

0.6610E  00 

0.3857E 

02 

0.3948E-01 

0.7623fc  00 

0,434">b 

02 

0.2638E-01 

0.8604E  00 

0.484t>E 

02 

o.r/^gE-oi 

0.9486E  00 
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KfcSULTS  OF 
POSITIUN 
0. 

0.5010E  01 
0.9890t  01 
0.l<»89b  02 
0.1978E  02 
0.2478E  02 
0.2963E  02 
0.3467F  02 
0.3857t  02 
0.444bh  02 
0.4946E  02 
0.5434C  02 
0.593^E  02 


12*>.  CPS.  DAIA 

AMPLITUDE 
O.lOOOt  01 
0.61L2t  00 
0.397yE-00 
0.2407E-00 
0.1650E-00 
0.9940E-0I 
0.6637E-01 
0.'t276E-0l 
0.3319E-01 
0.1986E-01 
0.l32bE-01 
0.8792E-02 
0.5825E-02 


ANALYSIS 

PHASE 
0. 

0.9480E-01 
0.1923t-00 
0.3007t-00 
0.401AE-00 
0.5174fc  00 
0.6247t  00 
0.72y2E  00 
0.8l43t  00 
0.9419E  00 
0.1058E  01 
0.ll67fc  01 
0.1267E  01 


RESULTS  OF 
POSITION 
0. 

0.5010E  01 
0.9890fc  01 
0.l<tb9fc  02 
0.1978E  02 
0.2473t  02 
0.2968E  02 
0.3<r67E  02 
0.3857E  02 
0.4443C  02 
0.4946E  02 
0.5A34t  02 
0.'i93<»E    02 


150.  CHS.  DATA 

AMPLITUDE 
O.IOOOE  01 
0-60746  00 
0.3926E-00 
0.2336t-00 
0.1602E-00 
0.9583fc-01 
0.63bOE-01 
0.4049b-01 
0.3118E-01 
0.1847E-01 
0.1223E-OI 
0.8017E-02 
0.524AE-02 


ANALYSIS 

PHASE 
0. 

0.1104E-00 
0.2240E-00 
0.3502E-00 
0.4678E-00 
0.6013fc  00 
0.7260f  00 
U.8487t  00 
0.94y3t  00 
0.1097E  01 
0.1236E  01 
0.1360E  01 
0.1480E  01 


RESULTS  OF 
POSITION 
0. 

0.5010E  01 
0.9890b  01 
0.lA89b  02 
0.L978E  02 
0.247dE  02 
0.296at  02 
0.3467E  02 
0.38i)7E  02 
O.^A'tbb  02 
0.4946F  02 
0.5434b  02 


175.  CPS.  DATA 

AMPLITUDE 
O.IOOOE  01 
0.6039E  00 
0.3878E-00 
0.23lOb-00 
0.1559E-00 
0.y264b-0l 
0.6096E-01 
0.3849b-0l 
0.2937E-01 
0.1727E-01 
0.1133E-01 
0.7340E-02 


ANALYSIS 

PHASb 
0. 

0.1252E-00 
0.2e)4lt-00 
0. 3971E-00 
0.5303E  00 
0.6822E  00 
0.8222E  00 
0.9609E  00 
0.1075E  01 
0.l240t  01 
0.1400E  01 
0.1537fc  01 
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0.b934E  02        0.4738t-02        0.1673E  01 


KfcSULTS  OF 

200.  CPS.  DATA 

ANALYSIS 

pusinuN 

AMPLITUDE 

PHASE 

0. 

O.lOOOfc  01 

0. 

0.3890E 

01 

0.7105E  00 

0.6277E- 

-01 

0.7890E 

01 

0.4723E-00 

0.9343E- 

-01 

0.117y& 

02 

0.33b9E-00 

0.2030E- 

-00 

0.i6f8C 

02 

0.2144E-00 

0.2297E- 

-00 

0.2067E 

02 

0.1493E-00 

0.3471t- 

-00 

0.2b67E 

02 

0.94i6E-0l 

0.3471E- 

-00 

0.29i>6E 

02 

0.6749E-01 

0.5120E 

00 

0.3456E 

02 

0.42b6E-0l 

0.5071E 

00 

0.1890E 

01 

0.8773E  00 

0.6160E 

00 

0.6900E 

01 

0.3269E  00 

0.6'i54E 

00 

0.11780 

02 

0.33596-00 

0.7a45t 

00 

0.167at 

02 

0.1987E-00 

0.8256E 

00 

0.2i67t 

02 

0.1331E-00 

0.9024E 

00 

0.2667E 

02 

0.7861E-01 

0.9821E 

00 

0.3l5rE 

02 

0.5138E-01 

0.1070E 

01 

0.3656t 

02 

0.3214E-01 

0.1136E 

01 

C'+O^bF 

02 

0.2435E-01 

0.1262E 

01 

Q.A634F 

02 

0.1422E-01 

0.1446E 

01 

0.bi35b 

02 

0.92b5E-02 

0.1623F 

01 

0.b623t 

02 

0.5932t-02 

0.1772E 

01 

0.6123E 

02 

0.3786b-02 

0.1924F 

01 

RESULTS  OF 

225.  CPS.  DAI 

rA  ANALYSIS 

POSITION 

AMPLITUDE 

PHASE 

0. 

O.IOOOE  01 

0. 

0.5010E 

01 

0.5966E  00 

0.1555E- 

-00 

0.9890e 

01 

0.3774E-00 

0.3145F- 

-00 

0.1489F 

02 

0.2214E-00 

0.4915E- 

-00 

0.l97at 

02 

0.1471E-00 

0.65'y-3t 

00 

0.247db 

02 

G.8623E-01 

0.8445E 

00 

0.296dfc 

02 

0.5591t-01 

0.10i8E 

01 

0.346/e 

02 

0.3465E-01 

0.1188E 

01 

0.3857E 

02 

0.2609t-01 

0.1326E 

01 

0.444bE 

02 

O.1513E-01 

0.1531E 

01 

0.4946E 

02 

0.9795t-02 

0.1724E 

01 

0.5434E 

02 

0.6210t-02 

0.1886E 

01 

0.5934E 

02 

0.3929t-02 

0.2054E 

01 

RESULTS  OF    2bO.  CPS.  DATA  ANALYSIS 
POSITION  AMPLITUDE  PHASt 

0.  O.IOOOE  01        0. 
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0.5010C 

01 

0.5918E    00 

0,1698E- 

-00 

0.98906 

01 

0.57106-00 

0.3422E- 

-00 

0.l489h 

02 

0.21546-00 

U.5355E 

00 

0.19/3E 

02 

0.1A18E-00 

0.7139E 

00 

0.2^7at 

02 

0.82316-01 

0.9211E 

00 

0.296ab 

02 

0.5281E-01 

O.IUIE 

01 

0.  W67fc 

02 

0.32376-01 

0.12986 

01 

0.3857E 

02 

0. 24206-01 

0.1450E 

01 

O.^^A'iC 

02 

0.13926-01 

0.1676E 

01 

o.^y^dt 

02 

0.39626-02 

0.18876 

01 

0.!><»3<>t 

02 

0.36116-02 

0.2063E 

01 

0.5934t 

02 

0.3b256-02 

0.22!>0E 

01 

RESULTS    06 

275.    CPS.    DATA 

ANALYSIS 

POSITION 

AMPLITUDE 

PHASE 

0. 

0.10006    01 

0. 

O.bOlOE 

01 

0.58686    00 

0.1820E- 

-00 

0.9890C 

01 

0.36466-00 

0.3659r- 

-00 

0.14896 

02 

0.20946-00 

0.5731E 

00 

0.19786 

02 

0.13646-00 

0.7639E 

00 

0.2478E 

02 

0.78306-01 

0.90626 

00 

0.29686 

02 

0.4964G-01 

0.1190E 

01 

0.34676 

02 

0.3000E-01 

0.1392E 

01 

0.3857E 

02 

0.22266-01 

0.1559E 

01 

0.44456 

02 

0.1263E-01 

0.1806E 

01 

0.49466 

02 

0.80a8E-02 

0.2035E 

01 

0.54346 

02 

0.49756-02 

0.2224E 

01 

0.59346 

02 

0.31036-02 

0.2433E 

01 

RESULTS  OF    290.  CPS.  OAIA  ANALYSIS 

POSITION           AMPLITUDE  PHASE 

0.  0.10006  01  0. 

0.1890E  01  0.82156  00  0.7499E-01 

0.3890E  01  0.6711F  00  0.1279E-OO 

0.7840E  01  0.4273E-00  0.2957E-00 

0.98906  01  0.35426-00  0.3704F-0O 

0.1178E  02  0.31776-00  0.46326-00 

0.1678E  02  0.1723E-00  0.6743F  00 

0.1878E  02  0.14386-00  0.7390E  00 

0.2067E  02  0.1175E-00  0.8041E  00 

0.25676  02  0.71896-01  0.1032E  01 

0.18906  01  0.87966  00  0.8481E-01 

0.O9006  01  0.51386  00  0.2732E-00 

0.11786  02  0.31776-00  0.4632F-00 

0.16736  02  0.18136-00  0.6776E  00 

0.21676  02  0.1174E-00  0.8747E  00 

0.26676  02  0.67006-01  0.11046  01 

0.3157E  02  0.42186-01  0.13156  01 
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0.36b6e 

02 

0.2527fc-01 

0.li>25E 

01 

0.4046t 

02 

0.1862E-01 

0.1699E 

01 

0.4634h 

02 

0.l047t-0l 

0.19'j6I: 

01 

0-513'jb 

02 

0.6673E-02 

0.2198E 

01 

0.5623t 

02 

O.^O-iQE-O? 

0.2392E 

01 

0.6123E 

02 

0.2519E-02 

0.261AE 

01 

KIZSULFS    OK 
PUSITION 
0. 

0-5010E  01 
0.9890t  01 
0.1489F  02 
0.197aE  02 
0.2473E  02 
0.2968E  02 
0.3467E  02 
0.3857b  02 
O.^^'iiijt  02 
0.4946E  02 
0.b434b  02 
0.5934b  02 


300.  CPS.  DATA 

AMPLITUDE 
O.LOOOE  01 
0.5827E  00 
0.3592E-00 
0.2043t-00 
0.1318E-00 
0.7494b-01 
0.4698E-01 
0.2799E-01 
0.2052b-01 
0.1148b-01 
0.7233E-02 
0.A39ab-02 
0.2721E-02 


ANALYSIS 

PHASE 
0. 

0.1927E-00 
0.3866E-00 
CbO'iSE  00 
0.8064e  00 
0.1041E  01 
0.1235E  01 
0.1469E  01 
0.lo48E  01 
0.1911E  01 
0.2160b  01 
0.2357t  01 
0.2587E    01 


RESULTS  OF 
POSITION 
0. 

O.bOlOE  01 
0.9890b  01 
0.1489b  02 
0.1978C  02 
0.247aE  02 
0.2968L  02 
0.346?b  02 
0.3857b  02 
0.4445b  02 
0.4946b  02 
0.5434b  02 
0.5934E  02 


325.    CPS.    DATA 

AMPHTUDb 
O.IOOOE  01 
0.5793G  00 
0.3549b-00 
0.2002E-00 
0.12a3b-00 
0.7241b-01 
0.4503E-01 
0.26b2E-01 
0.1922E-01 
0.1063F-01 
0.6664E-02 
0.3965E-02 
0.2439E-02 


ANALYSIS 

PHASE 
0. 

0.2037E-00 
0.4080E-00 
0.6385E  00 
0.8488E  00 
0.1095E  01 
0.1319E  01 
0.1541E  01 
O.1730E  01 
0.2006E  01 
0.2268t  01 
0.2470b  01 
0.2718E  01 


RESULTS  OF 

POSITION 
0. 

0.2000E  01 
0.3890b  01 
0.5900b  01 
0.7900C  01 


350.  CPS.  DATA 

AMPLITUDE 
O.IOOOE  01 
0.7901E  00 
0.6492E  00 
0.5592E  00 
0.4501E-00 


ANALYSIS 

PHASE 
0. 

0.1532t-01 
0.2546E-01 
0.6960E-01 
0.8690E-01 
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0.8890t  01  0.^210E-00  0.1095E-00 

O.iOaJt  02  0.3433E-00  0.1157E-00 

0.1278E  02  0.2766E-00  0. 1603E-00 

0. 1479b  02  0.2229E-00  0.160AE-00 

0.1679E  02  0.l804t-Q0  0.1B93E-00 

0.  0.1016E  01  0.2398E-00 

0.2000E  01  0.85ait  00  0.2449E-00 

0.3890E  01  0.7085E  00  0.3363E-00 

0.5900E  01  0.5)909E  00  0.3277E-00 

O.B890t  01  0.4292fc-00  0.3316E-00 

0.1089E  02  0.3532t-00  0.3779E-00 

0.1278E  02  G.2766b-00  0.3799E-00 

0.1478t  02  0.2223t-00  0.4543E-00 

0.1778E  02  0.159afc-00  0.4358E-00 

O.lOOOfc  01  O.lOAbb  01  0.547<^E  00 

0.2890b  01  0. 7551b  00  0.5'^74E  00 

0. 4890b  01  0.6199b  00  0.5474E  00 

0./900t  01  O.^'tlOb-OO  O.S-VS^fE  00 

0.289Ub  01  0.7892E  00  Q.b'iJ^E    00 

0.79006  01  0.4544E-00  0.7302E  00 

0.1278E  02  0.2767E-00  0.7843E  00 

O.niQE    02  0.1548E-00  0.7908E  00 

0.2267E  02  0.9863E-01  0.1012b  01 

0.2767b  02  0.3536b-01  0.l272t  01 

0.32b7b  02  0.3<»l9b-01  O.l-iOSE  01 

0.37;>6E  02  0.2000E-01  0.1740E  01 

0.4146E  02  0.1437E-01  0.1937E  01 

O.A734b  02  0.7886b-U2  0.2226f  01 

0.523bE  02  0.4900E-02  0.2498b  01 

0.b723E  02  0.2887E-02  0.2706E  01 

0.6223E  02  0.1781t-02  0.29?2E  01 


RESULTS  OF 

375.  CPS.  DATA 

ANALYSIS 

PUSIT[ON 

AMPLITUDE 

PHASE 

0. 

O.lOUOb  01 

0. 

0.'3010t 

01 

0.5714E  00 

0.2280E- 

-00 

0.9890E 

01 

0.3452b-00 

U.4555E- 

-00 

0.1489E 

02 

0.1912E-00 

0.7122b 

00 

0.1978E 

02 

0.1210E-00 

0.9456E 

00 

0.2478E 

02 

0.6739t-01 

0.1222E 

01 

0.296d£ 

02 

0.4l26b-0l 

0.1472E 

01 

0.346/E 

02 

0.2398b-01 

0.1719E 

01 

0.3857b 

02 

0.1713b-01 

0.1924E 

01 

0.44A5b 

02 

0.9336E-02 

0. 2231b 

01 

0.4946t 

02 

0.5784E-02 

0.2511E 

01 

0.5434t 

02 

0.3372b-02 

0.2737E 

01 

0.5934E 

02 

0.2l08h-02 

0.3025E 

01 
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RCSULIS  OF 

AOO.  CPS.  DAIA 

ANALYSIS 

POSITION 

AMPLITUDE 

PHASE 

0. 

O.IOUOE  01 

0. 

0.2000E 

01 

0.80^2b  00 

-0.52l7t- 

■01 

0.3890f- 

Oi 

0.6286E  00 

0.4326E- 

-01 

o.'fBgot 

01 

0.5639t  00 

-0.1t74t- 

-01 

0.8a90t 

01 

0.3472E-00 

0.6108E- 

-01 

0.l089t 

02 

0.2772E-00 

0.2606E- 

-01 

0.1278E 

02 

0.2237E-00 

0.3680E- 

-00 

0.i373t 

02 

0.2060E-00 

0.6496E- 

-01 

0.1778G 

02 

0.1303E-00 

-0.4567E- 

-01 

0. 

0.8268t  00 

0.1t.25E- 

-00 

O.IOOOE 

01 

0.7333t  00 

0.1150E- 

-00 

0.2000fc 

01 

0.6580E  00 

0.1802E- 

-00 

0.2890E 

01 

0.5797E  00 

0.1317E- 

-00 

0.3890E 

01 

0.b207E  00 

0.1783E- 

-00 

0.38906 

01 

0.*iA03E  00 

0.2367E- 

-00 

0.4890E 

01 

0.4940E-00 

0.8283E- 

-01 

o.^9ooi: 

01 

O.^^SSE-OO 

0.3862t- 

-00 

0.7900L 

01 

0.3651E-00 

0.3298E- 

-00 

0.8890C 

01 

0.3297E-00 

0.8058E 

00 

0.9890C 

01 

0.2980E-00 

0.3741E- 

-00 

0.1089E 

02 

0.2a31E-O0 

0.4?85E- 

-00 

0.1l78t 

02 

0.2^61E-00 

0.4171E- 

-00 

0.i278C 

02 

0.2237E-00 

0.A443E- 

-00 

0.1378E 

02 

0.2042b-00 

0.5b87E 

00 

0.1478E 

02 

0.1769E-00 

0.5138E 

00 

0.167BC 

02 

0.l4b^E-00 

0.3206t- 

-00 

0.i778E 

02 

0.12b9t-00 

0.6i80E 

00 

0.8890e 

01 

0.3160E-00 

0.5587E 

00 

0.200GE 

01 

0.6408E  00 

0.b930C 

00 

0. 

0.7659fc  00 

0.6618E 

00 

0.2890E 

01 

0.6600E  00 

0.8280E 

00 

0.7900E 

01 

0.373rE-00 

0.7686E 

00 

0.i273E 

02 

0.2237fc-00 

0.31.33E 

00 

0.1778E 

02 

0.1223E-00 

0.8259E 

00 

G.2267E 

02 

0.7682E-01 

0.107lt 

01 

0.2767E 

02 

0.^227E-01 

0.1363E 

01 

0.3257E 

02 

0.2556E-01 

0.1626E 

01 

0.37!J6e 

02 

0.1470E-01 

0. 1892E 

01 

0.4146E 

02 

0.1041t-01 

0.2109F 

01 

0.''t734t 

02 

0.5622E-02 

0. 

RESULTS  OF 
POSITION 
0. 

O.'iOlOt  01 
0.9890E  01 
0.1489E  02 
0.1978E  02 


425.  CPS.  DATA 

AMPLITUDE 
O.IOOOE  01 
0.'J620E  00 
0.3337F-00 
0.1806E-00 
0.1122E-00 


ANALYSIS 

PHASE 
0. 

0.2449E-00 
0.4904E-00 
0.7653E  00 
0.1020E  01 
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0. 

.2478E 

02 

0. 

,296BC 

02 

0, 

.3A67E 

02 

0. 

.3bb^E 

02 

0. 

.4445E 

02 

0. 

,^946E 

02 

0. 

.b43'»E 

02 

U.609<rE-0l 
0.-3635E-01 
0.2053b-01 
U.l^'tlfc-Ol 
0-76^8t-02 
0.47l6b-02 
0.2'>8dF-02 


0.1322E 
0.1'393E 
0.1875E 
0.2106E 
0.2449E 
0.2762fc 
0.3033E 


01 
01 
01 
01 
01 
01 
01 


RESULTS  OF 
PUSITION 
0. 

O.lOOOt  01 
0.2000E  01 
0.2890E  01 
0.3890E  01 
O.'^SQOE  01 
0.^900E  01 
0. f900b  01 
0.8890E  01 
0.9a90E  01 
O.IOOOC  01 
0.2890E  01 
0.7900E  01 
0.127at  02 
0.l77aC  02 
0.226?E  02 
0.276/b  02 
0.3257E  02 
0.37b6E  02 
O.^lAoG  02 
0.^73^t  02 
0.5235t  02 
0.5723b  02 


450.  CPS.  DATA 

AMPLITUDE 
O.lOOOb  01 
0.907'Jb  00 
0.8010E  00 
0.7275E  00 
0.6750E  00 
0.5934E  00 
0.bl70E  00 
0.4075E-00 
0.3b56E-00 
0.3230E-00 
0.9110G  UO 
0./282E  00 
O.^OrSE-OO 
0.240bE-00 
0.1292E-00 
0.797i>E-01 
0.4286E-01 
0.2533E-01 
0.1401E-01 
0.9718b-02 
0.b086E-02 
0.30:>OE-02 
0.1616b-02 


ANALYSIS 

PHASE 
0. 

0.3930fc-0l 
0.7820E-01 
0.25fbE-0l 
0.  l36lb-00 
0.1458E-00 
0.  1640E-00 
0.2206E-00 
0.281  rE-00 
0.3977E-00 
0.3V77E-00 
0.4243E-00 
0.<»840E-00 
0.6'iObb  00 
0.9314E  00 
0.1194E  01 
0.1502E  01 
0.1777E  01 
0.2068E  01 
0.2306E  01 
0.2659E  01 
0.2989E  01 
0.3258E  01 


KbSULTS  Ob 
POSITION 
0. 

O.IOOOC  01 
0.3890C  01 
0.^90UE  01 
0.  ?9UJb  01 
0.a89Ub  01 
0.1089E  02 
0.2890C  01 
0.7900E  01 
0.1278E  02 
0.177BC  02 
0.226K    02 


451.  CPS.  DATA 

AMPLITUDE 
O.IOOOE  01 
0.?817b  00 
0.6108E  00 
0.4768b-00 
0.3724C-00 
0.3262E-00 
0.2617E-00 
G.66b7b  00 
0.3724E-00 
0.2197E-00 
0.1180E-00 
0.7284b-0l 


ANALYSIS 

PHASE 
0. 

0.8'Jb7E-0l 
0.1056b-00 
0.1762E-00 
0.3260f-00 
0.3b78b-00 
0.3')78b-00 
0.3i76b-00 
0.5307E  00 
0. 6109b  00 
0.8921K  00 
O.llbSt  01 
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0. 27670 

02 

0.3913E-01 

0.1463E 

01 

0.3257h 

02 

0.2313E-01 

0.1738E 

01 

0.3756t 

02 

G.1278E-01 

0.2030t 

01 

0.4l^6E 

02 

0.8855t--02 

0.2269E 

01 

0.4734t 

02 

0.A636E-02 

0.2621t 

01 

O.'iZibE 

02 

0.2814E-02 

0.2952E 

01 

O.bJliL 

02 

0.l<»69E-02 

0.3220E 

01 

RESULTS  OF 

475.  CPS.  DATA 

ANALYSIS 

POSiriUN 

AMPLITUDE 

PHASE 

0. 

O.IOOOE  01 

0. 

O.'iOlOC 

01 

0.55aOE  00 

0.2613E- 

-00 

0.9890F 

01 

0.32/6E-00 

0.b230E 

00 

0.l489t 

02 

0.1752E-00 

0.8121t 

00 

0.l9/6t 

02 

0.l076t-00 

0.1085t 

01 

0.2478E 

02 

0.57'.4t-01 

0.l3y9E 

01 

0.2968E 

02 

0.3388b-01 

0.1681t 

01 

0.3467C 

02 

0.1848E-01 

O.iy/bE 

01 

0.3b57t 

02 

0.1269E-01 

0.22i8E 

01 

0.4445t 

02 

0.6635E-02 

0.2b77E 

01 

0.4946C 

02 

0.3951E-02 

0.2yL3E 

01 

0.b434E 

02 

0.2008E-02 

0.31!>5E 

01 

RESULTS  OF 

500.  CPS.  DAIA 

ANALYSIS 

POSITION 

AMPLITUDE 

PHASi 

_ 

0. 

O.IOOOE  01 

0. 

O.IOOOE 

01 

0.9222E  00 

0.4976E- 

-01 

o.20oor 

01 

0.8432fc  00 

0.l454t- 

-00 

0.2890fc 

01 

0.7665fc  00 

0.8940E- 

-01 

0.3890t 

01 

U.6862E  00 

0.1559E- 

-00 

0.4890fc 

01 

0.6081t  00 

0.1709E- 

-00 

0.5900E 

01 

0.5374E  00 

0.2Q50t- 

-00 

0.b900E 

01 

0.5323E  00 

0.27iOE- 

-00 

0.7900E 

01 

0.4229E-00 

0.3154E- 

-00 

0.88900 

01 

0.3882E-00 

0.3240E- 

-00 

0.9890C 

01 

0.3459E-00 

0.3390E- 

-00 

0.1089E 

02 

0.3l40b-00 

0.4395E- 

-00 

0.il78E 

02 

0.2805t-00 

0.4443E- 

-00 

0.12/3fc 

02 

0.2530E-00 

0.4767E- 

-00 

0.590ab 

01 

0.5653E  00 

0.5461E 

00 

0.2000E 

01 

0.9109t  00 

0.5911E 

00 

0.2a90E 

01 

0.7804E  00 

0.6472E 

00 

0./900C 

01 

0.4336E-00 

0.7I57E 

00 

0.1278E 

02 

O.2530t-00 

0.7lb7t 

00 

0.177aE 

02 

0.1345E-00 

0.10l6t 

01 

0.2267E 

02 

0.a2l7E-0l 

0. 1303E 

01 

0.2767E 

02 

0.4375E-01 

0.1627E 

01 

0.3257h 

02 

0.2566t-0l 

0.1923F 

01 
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0.3756t  02  0.1390E-01  0.2221E  01 

0.4l46t  02  0.95^6E-02  0.2474E  01 

0.A734E  02  0.b0l4E-02  0.2B48E  01 

0.'>723E  02  O.lt)0rE-02  0.3403E  01 


RESULTS  OF 

510.  CPS.  UAI/1 

i    ANALYSIS 

POSITION 

Ah'PLlTUDfc 

PHASE 

0. 

O.IOOOE  01 

0. 

0.2000t 

01 

0.7790E  00 

O.IO^IE- 

•00 

0.3890E 

01 

0.6276E  00 

O.lldOE- 

■00 

0.5900t 

01 

0.4718E-00 

0.1f90E- 

■00 

0.790Ut 

01 

0.3669b-00 

0.2981b- 

■00 

0.8890L 

OL 

0.3255b-00 

0. 4131b- 

■00 

0.I089C 

02 

0.2b99E-00 

0.3963b- 

-00 

0.127bE 

02 

0.2067E-00 

0.4b34E- 

-00 

0.1479E 

02 

0.1653E-00 

0.5661b 

00 

0.1679t 

02 

0.1268E-00 

0.6715E 

00 

0.2890b 

01 

0.6411E  00 

0.6715E 

00 

0. 7900b 

01 

0.35'5JE-O0 

0.8031E 

00 

0.1278t 

02 

0.2067E-00 

0.8895E 

00 

0.17/dt 

02 

0.109bE-00 

0.9768E 

00 

0.2267F 

02 

0.6666E-01 

0.1270E 

01 

0.2767fc 

02 

0.3541E-01 

0.l'598b 

01 

0.3257t 

02 

0.2070fc-01 

0.1901E 

01 

0.3756E 

02 

0.1119E-01 

0.2204E 

01 

0.^I46fc 

02 

0.7685E-02 

0.24'39F 

01 

0.4734E 

02 

0.4046t-02 

0.28A7F 

01 

0.'>2i5t 

02 

0.2372E-02 

0.3176b 

01 

0.t)723C 

02 

O.l223fc-02 

0.3401E 

01 

RESULTS  OF 

5?5.  CPS.  DAI/ 

k  ANALYSIS 

POSITION 

AMPLITUDE 

PHASt 

0. 

0. 1000b  01 

0. 

O.SOIOE  01 

0.'>516E  00 

0.2851E- 

■00 

0.9d90E  01 

0.3194E-00 

0.5652E 

00 

0.1489E  02 

0.1682E-00 

0.a772£ 

00 

0.l97aE  02 

0.1016E-00 

0.  1180E 

01 

0.2478b  02 

0.5375b-01 

O.lbiSb 

01 

0.2968b  02 

0.3ll8b-01 

0. 1827b 

01 

0.3467b  02 

0.1680E-01 

0.2138E 

01 

0.3857t  02 

0.1154E-01 

0.2404E 

01 

0.44^'3b  02 

0.6067E-02 

0.2813E 

01 

0.4y46L  02 

0.35t>lb-02 

0.3141E 

01 

0.'3'(»3^E  02 

0.1834b-02 

0.3377E 

01 

RESULTS  OF    550.  CPS.  DATA  ANALYSIS 
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POSITION 

AMPLITUDE 

PHASE 

0. 

O.IOOOE  01 

0. 

0.5010t 

01 

0.5^69E  00 

0.2938E- 

-00 

0.989Ub 

Oi 

0.3146E-00 

0.5801E 

00 

Cl^b-Jt 

02 

0.1639E-00 

0.9001E 

00 

0.1978E 

02 

0.97'>2E-01 

0.1213F 

01 

0.24f8b 

02 

0.Sl09E-0i 

0.l'>b7t 

01 

0.2968t 

02 

0.2909E-01 

0.1881E 

01 

0.3^6/E 

02 

0.l549t-0l 

0.2207E 

01 

0.3b57E 

02 

0.1054E-01 

0.2500E 

01 

0.4445t 

02 

0.546 3E-02 

0.2947E 

01 

0.-^946^ 

02 

0.3La6E-02 

0.32  76fc 

01 

0.t)434E 

02 

O.l679t-02 

0.3553E 

01 

RESULTS  OF 

575.  CPS.  DAIA 

ANALYSIS 

POSITION 

AMPLITUDE 

PHASE 

0. 

O.IOOOE  01 

0. 

O.SOlOt 

01 

0.5434E  00 

0.3002E- 

-00 

0.9890E 

01 

0.3116E-00 

0.5929fc 

00 

0.1489E 

02 

0.1610E-00 

0.9177b 

00 

0.1973t 

02 

0.9444t-0l 

0.1235E 

01 

0.2478t 

02 

0.49156-01 

0.1583E 

01 

0.2963E 

02 

0.2750E-01 

0.1908E 

01 

0.3467t 

02 

0.1440e-0l 

0.2242E 

01 

0.3a57E 

02 

0.9630E-02 

0.2561E 

01 

0.4445E 

02 

0.4799C-02 

0.3033t 

01 

0.4946E 

02 

0.2757E-02 

0.3362b 

01 

0.b434E 

02 

0.1436E-02 

0. 36926 

01 

RESULTS  OF 

600.  CPS.  DATA 

ANALYSIS 

pusirioN 

AMPLITUDE 

PHASE 

0. 

O.lOOOb  01 

0. 

0.5010C 

01 

0.5413E  00 

0.30/5E- 

-00 

0.9890E 

01 

0.3098b-00 

0.6100E 

00 

0.l4b9E 

02 

0.1594E-00 

0.9403b 

00 

0.1978E 

02 

0.9278E-01 

0.1260b 

01 

0.2478E 

02 

0.4824t-01 

0.1613b 

01 

0.296aE 

02 

0.26a2E-0I 

0.1936E 

01 

0.346?E 

02 

0-l389b-0l 

0.227ie 

01 

0.3b57E 

02 

0.9030b-02 

0.2601E 

01 

0.4445r 

02 

0.4318b-02 

0.3071E 

01 

0.4946b 

02 

0.2445E-02 

0.3378E 

01 

0.54  34b 

02 

0.1243b-02 

0.3f47E 

01 

RESULTS  Ob    625.  CPS.  DATA  ANALYSIS 
POSITION  AMPLITUDE  PHASE 


163 


0. 

O.IOOOE  01 

O. 

0.'>OlOb 

01 

0.5393E  00 

0.3171E- 

-00 

0.9890t 

01 

0.3074t-00 

0.6331E 

00 

O.l^eJL 

02 

0.1575fc-00 

0.9724t 

00 

0.19/8E 

02 

0.9137C-01 

0.1298E 

01 

0.2478t 

02 

0.47b3E-0l 

O. 1664E 

01 

Q.29bHr. 

02 

0.2645C-01 

0.l993t 

01 

0.3^67t 

02 

0.1371E-01 

0.2330E 

01 

0.38571: 

02 

0.8?llf-02 

0.2663E 

01 

0.^44'jE 

02 

0.4206E-02 

0.3122b 

01 

0.<^946C 

02 

0.23*j8E-02 

0.33/5E 

01 

0.5434t 

02 

O.H9bE-02 

0.3r77E 

01 

RESULTS  01- 

6b0.  CPS.  DAIA 

ANALYSIS 

PusiriuN 

AMPLITUDE 

PHASE 

0. 

0. 1000b  01 

0. 

0.3010t 

01 

0.b360E  00 

0.3273E- 

-00 

0.9890E 

01 

0.3028E-00 

0.6565E 

00 

O.l^BiE 

02 

0.1542E-00 

0. 1007b 

01 

0.197at 

02 

0.8890t-01 

0.1344E 

01 

0.247ab 

02 

0.4602E-01 

0.1726E 

01 

0.2968E 

02 

0.25b2E-01 

0.2072E 

01 

0.3467E 

02 

O.l328b-01 

0.2424E 

01 

0.38b/L 

02 

0.83'>6E-02 

0.2767E 

01 

0.444*56 

02 

0.4104E-02 

0.3235b 

01 

0.4946C 

02 

0.2400E-02 

0.34  36b 

01 

0.b434E 

02 

0.1240E-02 

0.3901E 

01 

RESULTS  OF 

675.  CPS.  UAIA 

ANALYSIS 

pusi riuN 

AMPLITUDE 

PHASE 

0. 

0. 1000b  01 

0. 

0.5010E 

01 

0.5320E  00 

0.3352E- 

-00 

0.98yOE 

01 

0.2971E-00 

0.6731b 

00 

0.  1469b 

02 

0.1500E-00 

0.1032b 

01 

O.I97aE 

02 

0.8545b-01 

0.1381E 

01 

0.2478L 

02 

0.4369b-0l 

0.1776E 

01 

0.296dE 

02 

0.2384b-0l 

0. 2141b 

01 

0.3467b 

02 

0.1234E-01 

0.2bl6E 

01 

0.38b7L 

02 

0.7683E-02 

0.2865E 

01 

0.4445C 

02 

0.3788b-02 

0.3386E 

01 

0.4946E 

02 

0.2387E-02 

0.3562b 

01 

0.5434E 

02 

0.1239E-02 

0.4132F 

01 

RESULTS  Ob    700.  CPS.  DATA  ANALYSIS 
PGSIIION  AMPLIIUUE  PHASE 

0.  0. 1000b  01        0. 
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0.2000E 

01 

0.7747t  00 

0.102t>fc- 

-00 

0.3890F: 

01 

0.5g49E  00 

0.ll95fc- 

-00 

0.5900fc: 

01 

0.4592t-00 

0.1853E- 

-00 

0./900E 

01 

0.3514E-00 

0.  I'iOSE- 

-00 

0.8890E 

01 

0.3098E-00 

0.2703E- 

-00 

0.10891: 

02 

0.2391E-00 

0.3a64E- 

-00 

0.l27dC 

02 

0.ia29E-00 

0.5185E 

00 

0.2890E 

01 

0.6349E  00 

0.4912E- 

-00 

O.IOOOE 

01 

0.8436E  00 

0.5893E 

00 

0.2890E 

01 

0.6244E  00 

0.70blfc 

00 

0.  /900t 

01 

0.3304E-00 

0.8338E 

00 

0.t273E 

02 

0.1830E-00 

0.8338E 

00 

0.l778t 

02 

0.9175E-01 

0.ll98t 

01 

0.2267E 

02 

0.bl49E-0l 

0.1553E 

01 

0.2767E 

02 

0.2600E-01 

0.19'>2E 

01 

0. 3257b 

02 

0.1384E-01 

0.2324b 

01 

0.3756t 

02 

0.7028E-02 

0.2716E 

01 

0.4146E 

02 

0.4261E-02 

0.30a2fc 

01 

0./t734E 

02 

0.2031E-02 

0.36'j5£ 

01 

0.t>23'iE 

02 

0.1404E-02 

0.3849E 

01 

RESULTS  Oh 

72'5.  CPS.  DATA 

ANALYSIS 

POSITIGN 

AMPLITUDE 

PHASE 

0. 

O.IOOOE  01 

0. 

0.5010E 

01 

0.5277E  00 

0.3474E- 

-00 

0.9a90E 

01 

0.2911fc-00 

0.69b8E 

00 

0.1489E 

02 

0.14b6fc-00 

0. Iu67fc 

01 

0.197at 

02 

0.8083E-01 

0.1422t 

01 

0.247ilt 

02 

0.4068E-01 

0.1821F 

01 

0.2963E 

02 

0.2l33b-Ol 

0.2188fc 

01 

0.3467E 

02 

0.1066E-01 

0.2;>78E 

01 

0.385/C 

02 

0.6286E-02 

0.2950E 

01 

0.444'5t 

02 

L.2764E-02 

0.3520E 

01 

0. 49466 

02 

0.2061E-02 

0.3799fc 

01 

RESULTS  01- 

750.  CPS.  DATA 

ANALYSIS 

POSITION 

AMPLITUDE 

PHASt 

0. 

O.lOOOfc  01 

0. 

O.SOiOE 

01 

0.3266E  00 

0.3b69E- 

-00 

0.9890E 

01 

0.2895b-00 

0.7158E 

00 

0.1489E 

02 

0.1446t-00 

0.1099E 

01 

0.197bE 

02 

0.7981E-01 

0.l4b6b 

01 

0.2478C 

02 

0.403t)E-01 

0.1863e 

01 

0.2968E 

02 

0.2113E-01 

0.2228E 

01 

0.3467E 

02 

0.1062E-01 

0.2610E 

01 

0.385fE 

02 

0.6186E-02 

0.2968t 

01 

0.'f445r 

02 

0.2629E-02 

0.3465t 

01 

0.4946b 

02 

0.1936E-02 

0.3865t 

01 
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RESULTS  OF 

775.  CPS.  DATA 

ANALYSIS 

PUSITION 

AMPLITUDE 

PHASE 

0. 

O.IOOOE  01 

0. 

0.501UE  01 

0.b242E  00 

0.3680E- 

-00 

0.9890r.  01 

0.2860E-00 

0.7399E 

00 

0.l489t  02 

0.1418E-00 

0.ll38r 

01 

0.l9/8r  02 

0.7797E-01 

0.l501b 

01 

0.247at  02 

0.3940fc-01 

0.1928E 

01 

0.2968fc  02 

0.2066fc-01 

0.2'i07F 

01 

O.i^bft    02 

G.l064b-0l 

0.2698E 

01 

0.3857E  02 

0.6H06-02 

0.3049E 

01 

O.A44^>t  02 

0.2a07E-02 

0.3493E 

01 

0.4946b  02 

0.1901E-02 

0.3949t 

01 

RESULTS 

OH    800.  CPS. 

DATA  ANALYSIS 

POSITION 

AMPLirUDb 

POSITION 

PHASE 

0. 

O.IOOOE  01 

0. 

0. 

O.IOOOE 

01 

0.8794E  00 

O.IOOOE 

01 

0.5856C-01 

0.1890b 

01 

0.7622E  00 

0.2000E 

01 

0.1408b-00 

0.2890E 

01 

0.6626E  00 

0.8890E 

01 

0.6665E  00 

0. 38906 

01 

0.5958t  00 

U.2890E 

01 

O.iyb9E-0O 

0.4900C 

01 

0.5060b  00 

0.3890E 

01 

0.2709E-00 

0.b900E 

01 

0.384  3r-00 

0.4890E 

01 

0.3579E-00 

0. /b90E 

01 

0.3384E-00 

0.5900E 

01 

0.4445E-00 

0.ia90t 

01 

0.5266G  00 

0.2890t 

01 

0.2253t-00 

0.6900b 

01 

0.3969E-00 

0. J9D0E 

01 

0.6029e  00 

0.1178E 

02 

0.1576E-00 

0.7900E 

01 

0.6029E  00 

0.1678L 

02 

0.7231b-0l 

0.1278f 

02 

0.9849E  00 

0.2167E 

02 

0.3945E-01 

0.177ab 

02 

0.1396fc  01 

0.2667E 

02 

0.1964E-01 

U.2267E 

02 

0.1767E  01 

0.3157b 

02 

0.1017E-01 

0.2767E 

02 

0.2215E  01 

0.i656E 

02 

0.5360G-02 

0.3257E 

02 

0.2618E  01 

0.4046E 

02 

0.3090b-02 

0.3756b 

02 

0.3044b  01 

0.4634t 

02 

0.1550e-02 

0.4146E 

02 

0.3416t  01 

0.5135C 

02 

0.9854t-03 

0.4734E 

02 

0.3872E  01 

0.5235E 

02 

0.4314E  01 

RESULTS 

OF    825.  CPS. 

DATA  ANALYSIS 

PUSH  lUN 

AMPLITUDE 

POSITION 

PHASE 

0. 

O.IOOOE  01 

0. 

0. 

0.5010C  01 

0.5174t  00 

0.5010b  01 

0.3845b-00 

0.9a90E  01 

0.2757b-00 

0.9890E  01 

0.7712E  00 

0. 1489b  02 

0.1331E-00 

0. 1489b  02 

0.1187b  01 

0.1978E  02 

0.7185E-01 

0.1978E  02 

0.1563b  01 

0.2478E  02 

0.3488t-01 

0.2478E  02 

0.2017C  01 
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0.2968t  02 
0.3467E  02 
0.38b7E  02 
0.44^5E  02 
0.4946C  02 
0.e)434E  02 
0.5934t  02 


0.l748t-0l 
0.9211E-02 
0.5257E-02 
0.2810e-02 
0.17556-02 
0.9205F-03 
0.4665t-03 


0.2968E  02 

0.3467E  02 

0.3857E  02 

0.444!>E  02 

0,4946E  02 


0.2440E  01 
0.2919E  01 
0.3300E  01 
0.3848E  01 
0.4264E  01 


RtiSULTS 

Oh    850.  CPS. 

DATA  ANALYSIS 

POSITIUN 

AMPLirUDE 

POSITIUN 

PHASE 

0. 

O.IOOOE  01 

0. 

0. 

0.5010E  01 

0.5156E  00 

0.5010E 

01 

0.3914E-00 

0.9890t  01 

0.2728E-00 

U.9890E 

01 

0.7816t  00 

0.1489E  02 

0.1108E-00 

0.1489E 

02 

0.1199b  01 

0.1978E  02 

0.7012E-01 

0.1978E 

02 

0.l576t  01 

0.247Bt  02 

0.3348E-01 

U.2478fc 

02 

0.2018E  01 

0.2968E  02 

0.1618G-01 

0.2968E 

02 

0,2436E  01 

0.3467E  02 

0.b259F-02 

G.3467E 

02 

0.2959E  01 

0.3857E  02 

0.4530E-02 

0.3857E 

02 

0.3344E  01 

0.4445E  02 

0.2451t-02 

0.4445t 

02 

0.4012t  01 

0.4946h 

02 

0.4432E  01 

RESULTS 

OF    875.  CPS. 

DATA  ANALYSIS 

POSITION 

AMPLITUDE 

POSITION 

PHASE 

0. 

O.IOOOE  01 

U. 

0. 

0.5010E 

01 

0.5144E  00 

0.5010E 

01 

0.4013E-00 

0.9890E 

01 

0.2709E-00 

0.9890fc 

01 

0.79856  00 

0.1489E 

02 

0.1297F-00 

0.1489t 

02 

0.1222E  01 

0.1978t 

02 

0.693aF-01 

0.19788 

02 

0.1602E  01 

0.2478t 

02 

0.3331E-01 

0.2473E 

02 

0.2032k  01 

0.2968E 

02 

0.1581fc-01 

0.2968t 

02 

0.2434E  01 

U.3467C 

02 

0.7791F-02 

0.3467E 

02 

0.2961E  01 

0.3857E 

02 

0.4165G-02 

0.3857E 

02 

0.3338E  01 

0.4445E 

02 

0.209aE-02 

0.4445E 

02 

0.4074E  01 

0.4946E 

02 

0.1330t-02 

RESULTS  OF 

900.  CPS.  DATA 

ANALYSIS 

POSITION 

AMPLITUDE 

PHASE 

0. 

O.IOOOE  01 

0. 

O.IOOOE  01 

0.84106  00 

0.7390E-01 

0.2000E  01 

0.7154E  00 

0.1439E-00 

0.2890E  01 

0.6225E  00 

0.2292E-00 

0.3890fc  01 

0,5508E  00 

0.2147E-00 

0. 48906  01 

C.4815E-00 

0.2907F-00 

0.590aL  01 

0.40706-00 

0.3625E-00 

0.7900F  01 

0.3031E-00 

0.4797E-00 
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0.2890t 

01 

0.5923fc  00 

o.62aot 

00 

o.79oae 

01 

0.3O31E-0O 

0.6281E 

00 

0.12/ae 

02 

0.l'>85fc-00 

0.l036fc 

01 

0.l77Bt 

02 

0.7559E-01 

0.1^73t 

01 

0.2i'67C 

02 

0.4036t-01 

0.l862t 

01 

0.2767t 

02 

0.l967fc-01 

0.2296E 

01 

0.3257E 

02 

0.9299E-02 

0.2695t 

01 

0.37561: 

02 

0.<»b66E-02 

0.321^E 

01 

O.^i^ht 

02 

0.2'»93E-02 

0.3568E 

01 

0.4734E 

02 

0.ll70t-02 

0.4287E 

01 

RESULTS 

OF    92b.  CPS. 

DATA  ANALYSIS 

POSITION 

AMPLITUDE 

POSITION 

PHASt 

0. 

O.IOOOE  01 

0. 

0. 

O.bOlOE  01 

0.5073E  00 

0.5010t 

01 

0.4239&-00 

0.9890h  01 

0.2619E-00 

0.9890t 

01 

0.8396E  00 

0.1489C  02 

0.1235E-00 

a.l4b9E 

02 

0.l29iE  01 

0.1978E  02 

0.6346E-01 

0.1978E 

02 

0.1692E  01 

0.2478b  02 

0.32l2r-01 

0.2A78E 

02 

0.2149E  01 

0.2y68tr  02 

0.1499E-01 

0.2968E 

02 

0.2557E  01 

0.3467t  02 

0.7'503t-02 

0.3467E 

02 

0.3091E  01 

0.3857fc  02 

0.4104E-02 

0.3857E 

02 

0.3445t  01 

0.4445C  02 

0.2038E-02 

RESULTS 

OF    950.  CPS. 

DATA  ANALYSIS 

POSITION 

AMPLITUDE 

POSITIUN 

PHASE 

0. 

O.IOOOE  01 

0. 

0. 

0.5010C 

01 

0.2589E-00 

0.5010E 

01 

0.4301E-00 

0.9890E 

01 

0.13l9n-00 

0.9890E 

01 

0.8480t  00 

0.l4a9t 

02 

0.6120E-01 

0.1489E 

02 

0.I306E  01 

0.1978C 

02 

0.3204E-01 

U.1978E 

02 

0.1/15E  01 

0.2478t 

02 

0.1553E-01 

U.2478C 

02 

0.2196E  01 

0.296flb 

02 

0. 7018E-02 

0.2968E 

02 

0.2608E  01 

0.i467L 

02 

0.3547E-02 

0.3467E 

02 

0.3186E  01 

0.i857E 

02 

0.l930e-02 

0.38b7t 

02 

0.3589E  01 

0,4445t 

02 

0.1043E-02 

0.4946t 

02 

0.3889P;-03 

RESULTS 

OF    975.  CPS. 

DATA  ANALYSIS 

posiriuN 

AMPLITUDE 

POSITION 

PHASE 

0. 

O.IOOOE  01 

0. 

0. 

O.bOlUb  01 

0.5004E  00 

O.SOlOfc  01 

0.4349fc-00 

0.9890b  01 

0.2533L-00 

0.9890E  01 

0.8523E  00 

0.1489L  02 

0.116  3E-00 

0.1489E  02 

0.1310E  01 

0.1978C  02 

0.6019E-01 

0.1978E  02 

0.l720t  01 

0.2478E  02 

0.2853E-01 

U.2473E  02 

0.2208E  01 
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0.2968b  02 
0.3467E  02 
0.3857b  02 
0.4A45C  02 


0.1246E-01 
0.6l23b-02 
0.32l!)E-02 
0.179'>C-02 


0.2968E  02 
0.3467E  02 
0.3857E  02 


0.2593E  01 
0.322lt  01 
0.3643E  01 


RESULTS  Of-   1000.  CPS.  DATA  ANALYSIS 


POSIIION 

AMPLITUDE 

PHASE 

0. 

O.IOOOE  01 

0. 

Q.5010E  01 

0.4991E-00 

0.4425t- 

-00 

0.989Ub  01 

0.2b20E-00 

0.8630E 

00 

0. 1489b  02 

0.1156b-00 

0.1322E 

01 

0.1978E  02 

0.5960E-U1 

0.1734E 

01 

0.2473b  02 

0.2790E-01 

0.2215C 

01 

0.29630  02 

0.1224E-01 

0.2563E 

01 

0.346/E  02 

0.'>726E-02 

0. J194E 

01 

0.i857C  02 

0.2786E-02 

0.3563E 

01 

0.444''ib  02 

0.1396E-02 

0.4698t 

01 

RESULTS  OF   1025.  CPS.  DATA  ANALYSIS 


POSITION 

AMPLITUDE 

POSITION 

PHASC 

0. 

0.10006  01 

0. 

0. 

0.5010E 

01 

0.4971G-00 

0. 5010b 

01 

0.4530E- 

-00 

0.9890C 

01 

0.2498r-00 

0.9890E 

01 

0.8320E 

00 

0.14a9E 

02 

0.1146E-00 

0.1489E 

02 

0.1352e 

01 

0.1978E 

02 

0.5901t-01 

0.19  78E 

02 

0.17736 

01 

0.247BE 

02 

0.27600-01 

0.2478b 

02 

0.2253E 

01 

0.2968E 

02 

0.1253F-01 

0.2968E 

02 

0.25896 

01 

0.3467t 

02 

0.5817E-02 

0.3467E 

02 

0.3194C 

01 

0.3857E 

02 

0.2785E-02 

0.3857E 

02 

0.3^91E 

01 

0.4445E 

02 

0.1048C-02 

0.4946E 

02 

0.5930F-03 

RESULTS  OF   1050.  CPS.  DATA  ANALYSIS 


POSITION 

AMPLITUDE 

PHASE 

0. 

O.IOOOE  01 

0. 

0.5010C 

01 

0.4933E-00 

0.4632E- 

-00 

0.9890t 

01 

0.245^E-00 

0.9003C 

00 

0.1489b 

02 

0.1117E-00 

0.1388E 

01 

0.  1978H 

02 

0.56986-01 

0.1822E 

01 

0.247aF 

02 

0.2675E-01 

0.2316E 

01 

0.2968F 

02 

0.1241E-01 

0.2677E 

01 

0.3467b 

02 

0.59256-02 

0.32816 

01 

0.38576 

02 

0.3O31E-02 

0.3510E 

01 
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RESULTS  OF   1075.  CPS.  DATA  ANALYSIS 


POSITION 

AMPLITUDE 

PHASE 

0. 

O.IOUOE  01 

0. 

O.tjOlOE 

01 

0.489lb-00 

0.4693E- 

-00 

0.9890E 

01 

0.2404t-00 

0.9083E 

00 

O.U89(: 

02 

0.10/3E-00 

0.1A09E 

01 

0.1978E 

02 

0.b378E-0l 

0.18'J2E 

01 

o.p-vTat 

02 

0.248i>E-01 

0.2i63t 

01 

O.*^960b 

02 

0.I146t-Ul 

0.2764E 

01 

0.346?E 

02 

0.5590fc-02 

0.3^09E 

01 

0.3857E 

02 

0.3116E-02 

0.362^E 

01 

RESULTS  or   1100.  CPS.  DAIA  ANALYSIS 


POSITION 

AMPLITUDIr 

PHASE 

0. 

O.lOOOt  01 

0. 

0.30  lot 

01 

0.4868E-00 

0.^732E- 

-00 

0. 98906 

01 

0.2378E-U0 

0.90976 

00 

0.  1489C 

02 

0.10^2t-00 

0.1410E 

01 

0.l973t 

02 

0.t)125E-01 

0.1849r 

01 

0.2478t 

02 

0.2299E-01 

0.2360t 

01 

o.?96ae 

02 

0.1023E-01 

0.2?84E 

01 

0.3^6 7 fc 

02 

0.48't6r-02 

0.3489E 

01 

0.3857E 

02 

0.2943E-02 

0.3729E 

01 

RESUI 

lTS  Oh   1125.  CPS. 

DATA  ANALYSIS 

POSITION 

AMPLITUDE 

POSITION 

PHASE 

0. 

O.IOOOG  01 

0. 

0. 

0.5010E  01 

O.A863E-00 

0.5010E 

01 

0.4798E-00 

0.9890C  01 

0.2378t-00 

0.9890b 

01 

0.9174E  UO 

0.1489h  02 

0.103311-00 

0.1489E 

02 

0.1413E  01 

0.1978t  02 

0.5054E-01 

U.1978E 

02 

0.1841E  01 

0.2478t  02 

0.222bfc-0l 

0.2A78E 

02 

0.2329E  01 

0.2968e  02 

0.9652b-02 

0.2968fc 

02 

0.2743t  01 

0.3A67E  02 

0.4192C-02 

0.3467E 

02 

0.3448E  01 

0. 38571-  02 

0.2744F-02 

0.38b7C 

02 

0.3753E  01 

0.4445E  02 

0.1242E-02 

RESULTS 

OF   1150.  CPS. 

DAIA  ANALYSIS 

POSITION 

AMPLITUDE 

POSITION 

0. 

O.lOOOt  01 

0. 

0.5010E  01 

0.4856C--00 

0.5010b  01 

0.9890E  01 

0.2375E-00 

0.9890E  01 

0.1489b  02 

0.1030t-00 

0.1489b  02 

0. 1978b  02 

0.5070E-01 

G.1978E  02 

0.2478t  02 

0.2254F-01 

0.2478E  02 

0.2968b  02 

0.9861C-02 

0.2968E  02 

PHASt 

0. 

0. 

.4915E- 

-00 

0, 

,9378E 

00 

0. 

.1436b 

01 

0. 

,1862E 

01 

0. 

2329E 

01 

0. 

2739E 

01 
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0.3467C  02  0.4ll8f:-02        0.3467E  02        0.33A9E  01 

0.33bfc    02  0.2801t-02        0.38b7h  02        0.3745e  01 

0.4445e  02  0.1119E-02 

0.49466  02  0.5139C-0J 


RESULTS  OF   1175.  CPS.  DATA  ANALYSIS 


POSITION 

AMPLITUDE 

POSITION 

PHASE 

0. 

O.IOOUE  01 

0. 

0. 

0.5010E  01 

0.4824E-00 

0.5010E 

01 

0.5050E  00 

0.9890C  01 

0.2345r-00 

0.9890E 

01 

0.9638E  00 

0.l489t  02 

0.1012E-00 

0.1489E 

02 

0.1471E  01 

0.1978E  02 

0.4997E-01 

0.1978E 

02 

0.1909E  01 

0.2478C  02 

0.2253t-0l 

0.2478E 

02 

0.2376E  01 

0.2968E  02 

0.1002E-01 

0.2968t 

02 

0.2822E  01 

0.3467c  02 

0.4433E-02 

0.3467t 

02 

0.3347E  01 

0.3b57t  02 

0.3054E-02 

0.3837E 

02 

0.382aE  01 

0.44450  02 

0.1135E-02 

RESULTS  OF   1200.  CPS.  DATA  ANALYSIS 


POSITION 

AMPL ITUDE 

POSITION 

PHASE 

0. 

O.IOOOE  01 

0  . 

0. 

0.50l0t  01 

0.477bE-00 

0.5010E 

01 

0.5149E  00 

0.9890E  01 

0.2293E-00 

0.9890E 

01 

0,9817E  00 

0.i489C  02 

0.9752E-01 

0.1489t 

02 

0.1496E  01 

0.l978h  02 

0.4776E-01 

0.1978E 

02 

0.1949E  01 

0.2478E  02 

0.2151E-01 

0.2478E 

02 

0.2426E  01 

0.2968E  02 

0.9386E-02 

0.2968E 

02 

0.2944t  01 

0.3467E  02 

0.4550E-02 

0.34t.?E 

02 

0.3440E  01 

0.3857E  02 

0.3160E-02 

0.3857E 

02 

0.4017E  01 

0.4445t  02 

0.1291E-02 

RESULTS  OF 

1229.  CPS.  DATA  ANALYSIS 

POSITION 

AMPLITUDE 

PHASE 

0. 

O.IOOOE  01 

0. 

O.IOOOE 

01 

0.8465t  00 

0.8747E- 

-01 

0.2000E 

01 

0.7240E  00 

0.I639E- 

-00 

0.2890E 

01 

0.6154E  00 

0.2841fc- 

-CO 

0.3890E 

01 

0.5116E  00 

0.2aAOE- 

-00 

0.2890E 

01 

0.61b4E  00 

0.39  30E- 

-00 

0.7900t 

01 

0.2911E-Q0 

0.8048E 

00 

0.1278E 

02 

0.l3a2E-00 

0.  I273t 

01 

0.1773t 

02 

0.b789fc-01 

0.  W82E 

01 

0.2267E 

02 

0.2781t-01 

0.2238E 

01 

0.276^t 

02 

0.1227r-0l 

0.2711E 

01 

0.3257c 

02 

0.4850E-02 

0.3299t 

01 

0.3756E 

02 

0.25/6E-02 

0.3814E 

01 
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0.4l46t  02        0.1744E-02        0.<f348E  01 


RbSULTS  OF 

12')0.  CPS.  DATA 

ANALYSIS 

PUSITIUN 

AMPLITUUE 

PHASE 

0. 

O.lOOOfc  01 

0. 

CSOlUf  01 

0.47l7t-00 

0.5253E 

00 

0.9890b  01 

0,2234t-00 

0.99i>5E 

00 

0.148911  02 

0.9349E-01 

0.1497E 

01 

0.197at  02 

0.4464fc-0l 

0. 1947b 

01 

0.247dt  02 

0.1962E-01 

0.239bE 

01 

0.296aE  02 

0.7029E-02 

0.2974E 

01 

0.i467t  02 

0.i893b-02 

0.3b08E 

01 

RESULTS 

OF   1275.  CPS. 

DATA  ANALYSIS 

POSITION 

AMPLITUDE 

POSITION 

PHASb 

0. 

O.lOOOb  01 

0. 

0. 

0.5010E  01 

0.4707E-00 

U. 5010b  01 

0.5349E 

00 

0.9890b  01 

0.2224E-00 

0.9B90E  01 

0.1015E 

01 

0.1489b  02 

0.9383b-01 

0.1489E  02 

0.1516C 

01 

0.19/8L  02 

0.4508r-0l 

0.197aE  02 

0. 1964b 

01 

0.24f8E  02 

0.2031b-01 

0.2470E  02 

0.2384E 

01 

0.2968b  02 

0.6850C-02 

U.296ab  02 

0.2904E 

01 

0.3467b  02 

0.3988b-02 

0.3467L  02 

0.3434E 

01 

0.3857b  02 

0.1922E-02 

0.4445b  02 

0.1133t-02 

RESULTS 

OF   1300.  CPS. 

DATA  ANALYSIS 

PUSITIUN 

AMPLITUDE 

POSITION 

PHASt 

0. 

O.IOOOE  01 

U. 

0. 

O.-iOlOb  01 

0.4683E-0J 

0.5010E  01 

0.5479b 

00 

0.9890E  01 

0.2188E-00 

0.9890E  01 

0.1C42E 

01 

0.1489t  02 

0.9282E-0i 

U.l489t  02 

0.1557b 

01 

0.1978C  02 

0.4496b-01 

0.l97bt  02 

0.2018E 

01 

0.24 rat  02 

0.21O9E-0I 

0.24/8E  02 

0.2439E 

01 

0.2968b  02 

0.7044C-02 

U.2968F  02 

0.2937E 

01 

0.3467E  02 

0.4503F-02 

0.3467E  02 

0.3458E 

01 

0.38b7L  02 

0.2048b-02 

0.4445b  02 

0.9658b-03 

0.4946E  02 

0.3748E-03 

RESULTS  OF   1325.  CPS.  DAIA  ANALYSIS 

POSITIUN  AMPLITUDE  POSITION  PHASE 

0.  O.IOOOE  01         0.  0. 

0. 5010b  01  0.4642E-00         0.5010t  01         0.5597t  00 
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0.9890E 

01 

0.2122E-00 

0.9890E 

01 

0.1061E 

01 

0.iA896 

02 

0.8931E-01 

0.1489E 

02 

0.1595E 

01 

0.1978E 

02 

0.4301E-01 

0.197bE 

02 

0.2080E 

01 

0.^478E 

02 

0.2080C-01 

U.2478E 

02 

0.2526b 

01 

0.2968E 

02 

0.6672t-02 

0.2968fc 

02 

0.3055E 

01 

0.3467C 

02 

0.48656-02 

0.3467t 

02 

0.360Afc 

01 

U.i8t>7t 

02 

0.2300G-02 

0.44^bb 

02 

O.lOlOfc-02 

0.4946t 

02 

0.4343E-03 

RESULTS 

or   i3b0.  CPS. 

DATA  ANALYSIS 

POSITION 

AMPLITUDE 

POSITION 

PHASE 

0. 

O.IOOOE  01 

0. 

0. 

0.5010fc  01 

0.4601E-00 

0.50l0fc  01 

0.5687fc  00 

0.9890E  01 

0.2058b-00 

U.9890E  01 

0.1062fc  01 

0.1489E  02 

0.8509e-0l 

0.1489E  02 

0.1605E  01 

0. 1978b  02 

0.4005E-01 

0.19/8E  02 

0.2107E  01 

0.2478b  02 

0.194"3G-01 

0.2478b  02 

0.258aE  01 

0.2968E  02 

0.3135E  01 

0.3467t  02 

0.3783E  01 

RESULTS  OF   1375.  CPS.  DATA  ANALYSIS 


POSITION 

AMPLITUDE 

POSITION 

PHASE 

0. 

O.IOOOC  01 

0. 

0. 

O.bOlOb  01 

0.4573b-00 

U.bOlOE  01 

0.5783E  00 

0.9890E  01 

0.2030E-00 

0.9B90t  01 

0.1056b  01 

0.1489E  02 

0.828/E-Ol 

0.1489E  02 

0.l595t  01 

0.1978E  02 

0.380'jb-Ol 

0.1978E  02 

0.2092E  01 

0.2478E  02 

0.1819E-01 

0.2478b  02 

0.2594t  01 

0.2968E  02 

0. 3016b  01 

RESULTS  OF   1400.  CPS.  DAIA  ANALYSIS 


POSITION 

AMPLlTUOb 

PHASb 

0. 

O.IOOOE  01 

0. 

0.5010E  01 

0.4541b-00 

0.5923E  00 

0.9690E  01 

0.2026t-00 

0.1061E  01 

0. 1489b  02 

0.8265E-0i 

0.1599b  01 

0.i973b  02 

0.377bE-01 

0.2082t  01 

0.24J'3C  02 

O.iir91b-0l 

0.2383E  01 

RESULTS 
POSir ION 
0. 
0.5010b  01 


OF   1425.  CPS.  DATA 
AMPLITUDE 
O.lOOOb  01 
0.4480E-00 


ANALYSIS 

PHASb 

0. 

0.60a3E  00 
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0.y090£  01 
0.1489E  02 
0.1978t  02 
0.247at  02 


0.2009t-00 
0.8197E-01 
0.^780t-0l 
0.1813E-01 


0.10/9E  01 

0.1628E  01 

0.2lllb  01 

0.2616E  01 


KESULFS  OF 

I4b0.  CPS.  DAIA 

ANALYSIS 

PUSITIUN 

AMPLirUOE 

PHASE 

0. 

O.IOOOE  01 

0. 

0.50l0fc  01 

O.A390fc-00 

0.6i90E 

00 

0.9a90C  01 

0.1963E-00 

0.1097F 

01 

0.1489fc  02 

0.790<tfc-01 

0.16b9E 

01 

0.1973E  02 

0.3662E-01 

0.2162E 

01 

0.^478b  02 

0.1  773t-0l 

0.2691E 

01 
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PART  2        IHE  COMHLtX  INVERSt  RELAXAFiON  LENGTH 


fkequency   omega 

RHO-«£AL 

ERRIJR- 

-RE 

RHO-IMAG 

ERRUR- 

-IM 

.I000bu2 

.6283e 

02 

.7995b- 

-01 

.lOOOE- 

-02 

.lb05E- 

-02 

.lOOOE- 

-04 

.I500b02 

.9425t 

02 

. 79y9E- 

-01 

.lOOOt- 

-02 

.26aOE- 

-02 

.1000b- 

-04 

.2500fc02 

.1571b 

0  3 

.8010E- 

-01 

.lOOOE- 

-02 

.4510E- 

-02 

.500ub- 

-04 

.50OO11O2 

.3i42L 

03 

.aO60E- 

-01 

.lOOOE- 

-02 

.9075F- 

-02 

.lOOOE- 

-03 

.65G0e02 

.4084E 

03 

.8105E- 

-01 

.lOOOE- 

-02 

.li75E- 

-01 

.3000E- 

-03 

W500fc02 

.47i2t 

03 

.8180E- 

-01 

.lOOOE- 

-02 

.1364E- 

-01 

.3000E- 

-03 

.  L000EU3 

.6283E 

03 

.8255E- 

-01 

.lOOOE- 

-02 

.i730E- 

-01 

.3000E- 

-0  3 

.il00cu3 

.69llh 

03 

.a350£- 

-01 

.90ooe- 

-03 

.1932E- 

-01 

.400Uh- 

-03 

.1250E03 

. 7854E 

03 

.8420E- 

-01 

.9000E- 

-03 

.2l7aE- 

-01 

.5000b- 

-03 

.I500fc03 

.9425b 

03 

.85g0b- 

-01 

.8000E- 

-03 

.25756- 

-01 

.3000E- 

-03 

.17bOLQ3 

.liOOE 

04 

.8795E- 

-01 

.800UE- 

-03 

.2902E- 

-01 

.200QE- 

-03 

-2000t03 

.1257C 

04 

.9000E- 

-01 

.6000E- 

-03 

.3195E- 

-01 

.4000E- 

-03 

.22'jOE03 

.1414E 

04 

.9180E- 

-01 

.6000E- 

-03 

.3562E- 

-01 

.3000E- 

-03 

.2'>00t03 

.1571b 

04 

.9350b- 

-01 

.AOOOE- 

-03 

.3915E- 

-01 

.2000b- 

-03 

.27bOE03 

.172aE 

04 

.9590E- 

-01 

.5000E- 

-03 

.4201E- 

-01 

.3000t- 

-03 

.2'V00i:03 

.i822E 

04 

.9740E- 

-01 

.5000b- 

-03 

.4372E- 

-01 

.4000t- 

-0  3 

.J000E03 

.I885E 

04 

.9810E- 

-01 

.6000b- 

-03 

.4492E- 

-01 

.300uE- 

-03 

. i250tU3 

.2042E 

04 

.lOOiE- 

-00 

.iOOOE- 

-02 

.4704E- 

-01 

.3000b- 

-03 

.3'j00l03 

.2199E 

04 

.1025E- 

-00 

.1000b- 

-02 

.4871E- 

-01 

.4000t- 

-03 

.3750L03 

.2356E 

04 

. 1044b- 

-00 

.8000E- 

-03 

.5190E- 

-01 

.4000E- 

-03 

.4000t03 

.2513b 

04 

.1071E- 

-00 

.  lOOOb- 

-02 

.5235E- 

-Oi 

.4000E- 

-03 

.'i250t03 

.2670t 

04 

.1099E- 

-00 

.loooe- 

-02 

.55  78E- 

-01 

.8000E- 

-03 

.4500fc03 

.2827b 

04 

.llllE- 

-00 

.lOOOb- 

-02 

.5823E- 

-01 

.7000tr- 

-03 

.4510L03 

.2834E 

04 

.1118E- 

-00 

.iOOOE- 

-02 

.5836E- 

-Oi 

.9000b- 

-03 

.4750C03 

.2985E 

04 

.ir39b- 

-00 

.1500b- 

-02 

.60  43E- 

-01 

.1000b- 

-02 

.""jOOOEOS 

.3142E 

04 

.1147E- 

-00 

.lOOOE- 

-02 

.6141E- 

-01 

.lOOOE- 

-02 

.^l00bU3 

.3204b 

04 

.li52E- 

-00 

.2000E- 

-02 

.6191E- 

-01 

.i200lr- 

-02 

."3250103 

.3299b 

04 

.il58E- 

-00 

.3000t- 

-03 

.638aC- 

-01 

.6000E- 

-03 

.5500E03 

.3456E 

04 

.11 r9E- 

-00 

.8000E- 

-03 

.6O03E- 

-01 

.1500b- 

-02 

.:>7bOC03 

.3613C 

04 

.1209E- 

-00 

.7000b- 

-03 

.6/89E- 

-01 

.1500E- 

-02 

.6000E03 

.3770E 

04 

.1228E- 

-00 

.lOOOb- 

-02 

.6955E- 

-01 

.1500E- 

-02 

.6250E03 

.3927E 

04 

.1242E- 

-00 

.1000b- 

-02 

.7132E- 

-01 

.  lOOOE- 

-02 

.6'J00h03 

.4084b 

04 

.1240E- 

-00 

.lOOOb- 

-02 

.7255E- 

-01 

.l20Ut- 

-02 

.67;>0E03 

.4241E 

04 

.1247E- 

-00 

.7000E- 

-03 

.7382E- 

-01 

.1500E- 

-02 

.7000h03 

.439ab 

04 

.1277E- 

-00 

.lOOOb- 

-02 

.7546E- 

-01 

.15001:- 

-02 

.7250h03 

.4555b 

04 

.1301E- 

-00 

.1000b- 

-02 

.7721F- 

-01 

.2000E- 

-02 

.7500E03 

.4712E 

04 

.1314E- 

-00 

.2000E- 

-02 

.785iE- 

-01 

.2000b- 

-02 

.7750E03 

.4869E 

04 

.1314E- 

-00 

.700UE- 

-02 

.7y29E- 

-01 

.2000t- 

-02 

.8000h03 

.5027C 

04 

.1358E- 

-00 

.lOOOE- 

-02 

.8237E- 

-01 

.lOOuL- 

-02 

.8250t03 

.5184b 

04 

.1343E- 

-00 

.3000E- 

■02 

.8280E- 

-01 

.2000E- 

-02 

.B'i00E03 

.5341E 

04 

.13a7E- 

-00 

.2000b- 

-02 

.8444E- 

-01 

.30006- 

-02 

.B750E03 

.549bE 

04 

.1398E- 

-00 

.2000E- 

-02 

.8620E- 

■01 

.3000E- 

-02 

.9000h03 

.5655b 

04 

.1415E- 

-00 

.lOOOE- 

-02 

.a750E- 

-01 

.2000c- 

-02 
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.9250003 
,9i>U0h03 
,97bOL03 
,  1000CU4 
,  l02btU^ 
,  10'i0r04 
,l075t0'^ 
.ilOOCO^ 

,  il50E0^ 
,  117!>L04 
.  l200tU4 
, I229fc04 
,12t>0fc0^ 
,l275fcOA 
.  1300C04 
r32bt04 

,i3^oi;o^ 

l^75fcO<^ 
:l400b0^ 

I425h04 
,i4t>0t04 


.5812E  04 

.bybye  04 

.6126f:  04 
.6283E  04 
.6440b  04 
.6!)97C  04 
.6754t  04 
.6911E  04 
.7069t  04 
.7226r  04 
.7383E  04 
.7540C  04 
.7722b  04 
.7854C  04 
-aOllt  04 
.bl68F  04 
.8325E  04 
.8482C  04 
.8639b  04 
.a796E  04 
.8954b  04 
.9ilih  0^ 


1426E-00 
1453E-00 
1467E-00 
1483E-00 
l501b-00 
1500E-00 
1494E-00 
1523E-00 
1533E-00 
1536E-00 
1529E-00 
1525E-00 
1580E-00 
i609E-00 
1586E-00 
1588E-00 
i574E-00 
1611E-00 
1635E-00 
1640E-00 
1639E-00 
1655E-00 


.inOOE-02 
.lOOUb-02 
.2000b-02 
.2000E-02 
.2200E-02 
.2300E-02 
.2300E-02 
.2400t-02 
.2600E-02 
.2600E-02 
.2800E-02 
.3000E-02 
.3200b-02 
.3200E-02 
.3400E-02 
.3500E-02 
.3500E-02 
.3800fc-02 
.3800E-02 
.3900E-02 
.'tOOOE-02 
.4000E-02 


.8706E-01 
.9030E-01 
.910'5E-Ol 
.9200E-01 
.9028E-01 
.9220E-01 
.g'iOlb-Ol 
.9650E-01 
.9600E-01 
.9b42E-01 
.97lbE-Ol 
.  1006E-00 
.  1030E-00 
.9977E-0i 
.9fl95E-01 
.lOOlE-00 
.1038E-00 
.  1070E-00 
.1U41E-00 
.  lUb5E-00 
.1071E-00 
.1097b-00 


.lOOOE-02 
.1500E-02 
.l500b-02 
. l800b-02 
.l80Ut-02 
.200Ut-0^ 
.2300E-02 
.2300E-02 
.2700t-02 
.2800E-02 
.2800e-02 
.3000E-02 
.3200E-02 
.i300E-02 
.3300E-02 
.340ub-02 
.3400E-02 
.3700L-02 
.370Ub-02 
. 3800fc-02 
.4000E-02 
.4000b-02 
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PART  3   THC  SQUARED  CUMPLEX  INVERSE  RELAXATION  LENGTH 


UMCGA   OMEGA**^   RH0ft«2-IM  ERRQR-IM   RH0»»2-RE  ERROR-RE 


.6283t02  .39485  04  .2886E-03  .5209E-05  .6389E-02  .159yE-03 

.9425t02  .a883b  04  .4287E-03  .696GE-0b  .6391E-02  .  l600t-03 

.1571E03  .2467E  05  .722'JE-03  .1703t-04  .6396E-02  .i607E-03 

.3142E03  .9870E  05  .I463E-02  .3427E-04  .64I4E-02  .163UE-03 

.40a4r03  .i668E  06  .1903E-02  .7213E-04  .643IE-02  .169lt-03 

.47l2t03  .222Ih  06  .2232E-02  .7636E-04  .6505E-02  .1713E-03 

.6283C03  .3948h  OG  .2856E-02  .8412E-04  .6blbe-02  .l75oE-03 

.69llt03  .477/E  06  .3226E-02  .10160-03  .6b99E-02  .i65d£-03 

./8b4t03  .6i68£  06  .3668E-02  .i234E-03  .6ol5E-02  .1733C-03 

.9425h03  .88830  06  .4424E-02  .92  74£~04  .67I6E-02  .15290-03 

.1100E04  .1209E  07  .5105E-02  .8161E-04  .68930-02  .1523E-03 

.12'>7L04  .1579E  07  .57t>lE-02  .11030-03  .70790-02  .1336E-03 

.14141:04  .1999E  07  .6540E-02  .9/82E-04  .71580-02  .131bt-03 

.1571E04  .2467E  07  .7321E-02  .8438E-04  .7zlOE-02  .  1279E-03 

.1728004  .2986F  07  .80570-02  .99550-04  .74320-02  .12110-03 

.1822004  .3320E  07  .85170-02  .1216E-03  .75750-02  .13240-03 

.ia85L04  .3553E  07  .8813E-02  .1128E-03  .7606E-02  .14470-03 

.^042004  .4170fc  0/  .9417E-02  .15410-03  .76070-02  .22840-03 

.2199004  .4836b  07  .99860-02  .1794E-03  .8134E-02  .24400-03 

.2356E04  .5552E  0/  .1084E-01  .1666E-03  .82060-02  .20860-03 

.2513004  .63170  07  .1121E-01  .19040-03  .87300-02  .25610-03 

.2670004  .71310  07  .12260-01  .28740-03  .89730-02  .30910-03 

.2827fc04  .79940  07  .12940-01  .27200-03  .895rE-02  .303bt-03 

.2834004  .80300  0/  .13050-01  .3180E-03  .9102E-02  .32870-03 

.2985004  .89070  07  .13760-01  .40900-03  .93170-02  .46250-03 

.3142004  .98700  0^  .14090-01  .35220-03  .9382E-02  .35220-03 

.3204004  .10270  08  .14260-01  .52410-03  .9433E-02  .6093E-03 

.3299004  .10880  08  .1479E-01  .2412E-03  .93290-02  .26190-03 

.3456004  .11940  08  .15580-01  .45950-03  .9549E-02  .38680-03 

,3613004  .13050  08  .16410-01  .457/0-03  .luOlF-Ol  .37290-03 

.3770004  .14210  03  .17080-01  .507b£-03  .1024E-01  .45430-03 

.3927004  .15420  08  .17720-01  .39110-03  .10340-01  .39110-03 

.4084004  .16680  08  .18000-01  .44280-03  .1012E-01  .42220-03 

.4241004  .17990  03  .18410-01  .47/50-03  .lOUE-Ol  .39610-03 

.4398004  .19340  08  .19270-01  .53400-03  .10620-01  .48180-03 

.4555004  .20750  08  .20090-01  .67490-03  .10970-01  .56910-03 

.4712004  .22210  08  .20630-01  .83950-03  .11090-01  .8395t-03 

,4869^04  .23710  08  .20840-01  .84270-03  .  Iu960-01  .842/0-03 

.5027004  .252/0  08  .22370-01  .43630-03  .11650-01  .43630-03 

.5184004  .268/0  08  .22240-01  .40360-04  .11180-01  .47450-03 

.5341004  .28520  08  .23420-01  .11700-02  .12090-01  .10610-02 

,5498004  .30230  08  .24100-01  .11840-02  .12110-01  .  1076E-02 

.5655004  .31980  08  .24760-01  .74090-03  .1236E-01  .6330E-03 
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.58l2b04 

.3378b 

08 

.2506t- 

-01 

.4610E- 

-03 

.12626- 

-01 

.46106- 

-03 

.3969604 

.3363C 

08 

.262bt- 

-01 

.6166E- 

-03 

.1297E- 

-01 

.5616b- 

-03 

.6126L04 

.3753E 

08 

.2671E- 

-01 

.a043b- 

-03 

.1323E- 

-01 

.a599E- 

-03 

,6233t04 

. 5948E 

08 

.2728b- 

-01 

.90186- 

-03 

.13526- 

-01 

.92436- 

-03 

.6440L04 

.4148E 

08 

.2710E- 

-01 

.93  7b£- 

-03 

.143/6- 

-01 

.98546- 

-03 

.6597t04 

.4352L 

08 

.2766E- 

-01 

.10246- 

-oz 

.1400F- 

-01 

.10596- 

-02 

.6754t04 

.4562E 

08 

.28406- 

-01 

.I12b6- 

-02 

.13316- 

-01 

.1125E- 

-02 

.69Ut04 

.4777t 

03 

.2939b- 

-01 

.11646- 

-02 

.13886- 

-01 

.1175C- 

-02 

.  ^0691; 04 

.4996C 

08 

.2944b- 

-01 

.132/6- 

-02 

.14296- 

-01 

.13166- 

-02 

./226L04 

.'>221C 

08 

.2932E- 

-01 

.135/6- 

-02 

.1450E- 

-01 

.1333E- 

-02 

./3B3L04 

.54!>0E 

08 

.29/26- 

-01 

.14006- 

-02 

.13956- 

-01 

.1400b- 

-02 

.7')40l04 

.56a6C 

08 

.3066E- 

-01 

.15136- 

-02 

.1313E- 

-01 

.1518C- 

-02 

. ITZZcO^ 

.b963E 

08 

.32^4E- 

-01 

.16706- 

-02 

.1436E- 

-01 

.16706- 

-02 

.7804hU4 

.6l68fc 

0  8 

.32106- 

-01 

.17006- 

-02 

.1593E- 

-01 

.16836- 

-02 

.B01lt04 

.6418C 

08 

.31396- 

-01 

.1720E- 

-02 

.1537E- 

-01 

.17326- 

-02 

.8168C04 

.6672b 

08 

.31/9E- 

-01 

.17806- 

-02 

.15iaE- 

-01 

.17926- 

-02 

.8325tU4 

.6931c 

U8 

.3269E- 

-01 

.179/6- 

-02 

.13996- 

-01 

.18086- 

-02 

.84a2E04 

.7i95E 

08 

.3447E- 

-01 

.20056- 

-02 

.1451E- 

-01 

.20166- 

-02 

.8639C04 

.7464C 

08 

.3405E- 

-01 

.20016- 

-02 

.1589E- 

-01 

.20136- 

-02 

.8796L04 

.773bE 

08 

.3461b- 

-01 

.2069E- 

-02 

.1^756- 

-01 

.20816- 

-02 

.8954fc04 

.8017b 

08 

.351UE- 

-01 

.2168E- 

-02 

.l54iE- 

-01 

.2l6db- 

-02 

.9111C04 

.a30GE 

08 

.3630E- 

-01 

.22016- 

-02 

.1534E- 

-01 

.22016- 

-02 
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